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During the years of our First ISLSCP Field Experiment (FIFE) mvestigahon 
"Stadfe 0/ L Net Surface Radiative Flux from Satellite Radmnces ^ f ^ 

(NAG5-9oV our efforts have focused on radiadve tranfer inodel 
calibration of VISSR and AVFIRR solar channels, development and refinemen 
of algorithms to estimate downward solar (shortwave) and terrestria < 
irradfances at the surface, including photosyntheticaliy available radiation (PAR) 
and surface albedo, verification of these algorithms using in situ measureinents, 
production of maps of shortwave irradiance, surface albedo, and related products 
[eg cloud albedo, fractional cloud cover), analysis of the temporal variability 
dfumal, seasonal) of shortwave irradiance over the FIFE site, development of a 
spectroscopy technique to estimate atmospheric total water vapor amount, and 
Ifudy of opLum toear combinations of visible and near-infrared refle^n es 
for Ktimathig the fraction of PAR absorbed by plants. Among these activities the 
last two and^the one concerning PAR estimation from space were not part of our 
original technical plan, but they appeared interesUng to ^dertake 
addressing FIFE issues. Most of our investigations objectives have been 
acTO^plished, but we have not been successful in achieving useful acimraaes on 
esHmates of downward longwave irradiance, and we have done too little wof ot 
[ he longwave radiation lost by the surface. The satellite sensors available 
unfortunately, did not appear sufficiently sensitive to ‘I;® ^ 
governing the variability of the longwave fluxes, making ill-conditioned and, 
fherefore arduous the inverse problem. The sciendfic results of the investigation 

are summarized below. 

1. Radiative Transfer Model Validation 

The radiative transfer models of Tanr4 et al. (1985) and Morcrette (1984), or 
simplified schemes based on these models, are used in our algorithms to estimate 
surface radiation fluxes from satellite radiances. Although ^ 

been verified against exact calculations or line-by-line models, they have not been 
validated against in situ measurements. During HFE, concomitant radiosonde 
observation! sky photographs, aerosol turbidity, and radiation measurements 
were made, offering the opportunity to compare model output to surface- 
measured fluxes. The comparisons were made in clear sky conditions since the 
available data did not permit adequate characterization of some imi^rtant cloud 
parameters. The results are described in Froum et al. (1990; see Appendix 
Figures 1 and 2, which display scatter plots of calculated versus measured fluxes 
summarize the results. The agreement is good between calculated and measure 
fluxes, with correlation coefficients above 0.98 and standard deviations of 
23 2 Wm-2 (2 7%) and 13.0 Wm-2 (3.7%) for shortwave and longwave irradiances, 
respectively. The model of Tanr^ et al. (1985) overestimated shortwave irradiance 
by 13 2 Wm-2 on average, whereas the model of Morcrette (1984) underestimated 
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1 • Viv 7 4Wm-2 If significant, these biases may compensate 
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sSubm^ orth?rdela of Tfnr^ et al. (1985) and Morcrette (1984) (or radiation 
budget studies, at least in clear sky conditions. 

2 Calibration of VISSR and AVHRR Solar Channels 

■ We have performed a sensitivity study of the net f "sL'v 

the surface to calibration (Gautier and Frouin, 1988; see Appendix 2). This study 
MS sto^ hat in overcast conditions, a 10% loss in sensor sensitivity (maease 
m el)“ates into uncertainties of up to -70Wm-2 on instantaneous values 
and of up to -15 Wm-l on monthly averages. In tropical regions, the mon y 
avtageS uncertainty is typically -14 Wm-2 (Figure 3). This uncertainty appear as 
tri gn opposite to that of the calibration bias; in other words a pos^e 

bias orcalilrafon'lLuces a negaUve bias on net shortwave irradiance (cloud 
albedo is overestimated). , . j- 

VISSR and AVHRR calibration was achieved by applying the in irec 
method described in Frouin and Gautier (1987). This method, which uses space 
and the White Sands Monument area in New Mexico as calibra ion targe , 
“mputesT rain measured at satellite altitude using the model of Tanr4 e 
al (1985) The relevant atmospheric characteristice are estimated from either 
SmSical data (aerosol size-frequency distribution and refractive index) o 
observatiOTS at nearby meteorological stations (water vapor amount and aerosol 
gure 4 sLws the GOES-6 VISSR calibration gain (relates counP 

square! to reflectance) obtained (open circles) and how “ 'ifaf 

reported by other investigators (the compilaflon was made by Whitloi* et al 
1990) Some erratic fluctuations are observed from one gain estimate to the nex , 
in some instances reaching 10%. Changes of such amplitude, however are at the 
limit of the uncertainty^level. In addition to the relatively high-frequency 
fluctuations, the gain exhibits low-frequen^ variations vn* 
in spring and early summer. From July to December 1987, the gain increased Dy 
approxiLtely 36%. Compared to the prelaunch gam, the gam at the end of 198 

reads higher by 64%. 


3. Algorithm Development, Refinement, and Validation 


3.1 Downward Shortwave Irradiance, Including PAR, and Surface Albedo 

We have applied the satellite method of Gautier et al. (1980) to GOES-6 
VISSR data acquired during HFE. The version used accounts for 
and absorption, which was accomplished based on the model of Tanr6 et ah 
(1985). Throbjective was to assess whether the method can provide quantitative 
inforLtion on the temporal variability of downward shortwave irradiate (or 
insolation) and PAR during the experiment. In the calculahons smface a 
each observation time during the day was determined for each IFC from the 


minimum brightness count over the entire IFC. We have indicated (Frouin and 
Gautier, 1990; see Appendix 3) how the Gautier et al. (1980) method can be 
modified to yield PAR estimates. In this case, cloud absorption vanishes in the 
model equations and clear sky PAR is computed using the coefficients given by 
Frouin et al. (1989). The computational procedure included checks of the satellite 
navigation and of the data quality, and calibration of the VISSR visible channel 
(see above). Appendix 4 summarizes the results obtained for the five IFCs. The 
satellite estimates of insolation and PAR were in good agreement with in situ 
measurements. During IFC-2, for instance, the correlation coefficients were above 
0.95 with the standard errors of estimate of 21.6 Wm'^ (9%) and 8.2 Wm'^ (6.5%) 
for daily insolation and PAR, respectively (Figure 5). These errors are acceptable 
for climate studies. The satellite estimates also described well the diurnal 
variability of PAR and insolation (Figure 6). Using the ratio of PAR and 
insolation determined in situ with the satellite estimate of insolation, we 
obtained daily PAR values with a slightly higher r.m.s error than when PAR was 
directly estimated from the satellite data. This indicates than the large scale 
satellite climatologies of insolation now produced within the framework of 
ISCCP can be simply converted into useful climatologies of PAR for global studies 
of photosynthetic activity, although a more direct use of the satellite data would 
be optimum. The various techniques so far proposed to estimate PAR from 
satellite data, including the one of Frouin and Gautier (1990), are discussed in 
Frouin (1993) (Appendix 5). 

3JZ Downward Longwave Irradiance 

The net longwave irradiance flux at the surface, owing to its small 
variability, is less important than the shortwave one. It cannot be neglected, 
however, particularly in the presence of clouds that decrease surface insolation 
while increasing the downward component of the longwave flux. Several 
methods have been proposed to estimate the downward longwave flux at the 
surface from satellite measurements (see, for instance, Darnell et al., 1983; Frouin 
et al., 1988; Gupta, 1989). Although some of them have been used to generate 
global fields of net longwave flux, few validation studies have been made. We 
therefore investigated the performance of two of them (Frouin et al., 1988; Gupta, 
1989) by comparing their outputs with carefully taken in situ measurements. 
Computations using the same atmospheric temperature and moisture profiles 
indicated that, for clear sky conditions, the two model outputs are correlated. The 
values obtained using Gupta's (1989) model are systematically higher by 15 Wm-2- 
Since Frouin et al. (1988) uses Morcrette's (1984) model, and since this model 
slightly underestimates downward longwave irradiance (see section 1), we may 
conclude that the performance of the two models is similar, Gupta's (1989) model 
overestimating downward longwave irradiance in somewhat warm and humid 
atmospheres. Computations using vertical sounder (TOVS) data as input did not 
show good agreement with in situ measurements, especially in cloudy conditions 
(Figure 7). This was expected since other uncertainties are added: that of the 
retrieved profile itself and the uncertainty resulting from the distance between 


the TOYS observations and the pyrgeometer measurements. Details of this study 
are given in Br6on et al. (1990; see Appendix 6). 


3.3 Atmospheric Water Vapor Amount 

We have shown (Frouin and Middleton, 1990; see Appendix 7) that 
viewing the earth's surface in two spectral channels, one narrow, the other wide, 
centered on the same wavelength at the water vapor absorption maximum near 
940 nm yields accurate atmospheric total water vapor amounts independently o 
the surface reflectance properties (Figure 8). The concept was verified iwing 
concomitant SE-590 spectrometer and radiosonde data acquired dunng F^ IFCs. 
Under varied geometries and atmospheric conditions, and for several surface 
targets, the relationship between radiance ratio in the narrow and specttal 
chLtels was stable, foUowing an expected exponential law (Figure 9). The SE-590 
instrument, however, did not view the surface from above the atmosphere or 
even above the boundary layer, where most of the water vapor is concentrated, 
and therefore did not permit verification of the technique in actual conditions of 
aircraft or satellite viewing. In such conditions, simple physics showed that the 
effect of photons backscattered by aerosols is to lower unacceptably the retrieved 
water vapor amount when the aerosols are located above the boundary layer. 
Airborne experiments, on the other hand, were carried out (Frouin al., 1990; 
Appendix 8) and have demonstrated the suitability of the method over both land 

and sea. 


3.4 Fraction of PAR Absorbed by Plants 

Using the SAIL model (Verhoef, 1984), we have simulated for varied, yet 
realistic soil and canopy parameters, namely leaf optical properties, soil 
reflectance. Leaf Area Index (LAI), and Leaf Inclination Distribution Functi^ 
(LIDF), above canopy visible and near-infrared reflectances and daily-averaged 
(weighted by incident radiation) fraction of PAR absorbed by plants, AFAR. The 
objective was to find linear combinations of visible and near-infrared reflectances 
that improve the accuracy of APAR estimates when compared to that obtained 
using Normalized Difference Vegetation Index (NDVI). Figures 10 and 11 display 
some results (see, also Podaire et al., 1991; Appendix 9). The calculations were 
performed for LAIs of 0.01, 0.05, 0.1, 0.2, 0.5, 1, 2, and 5, erectophile, planophile 
and spherophile canopies, soil reflectances of 0.1, 0.2, 0.3, and 0.4, and typical 1^ 
optical properties. The soil reflectance was assumed white spectrally, and the 
various LIDFs were considered separately in the regressions. We see in Figure 10 
(bottom) that the influence of the background on the daily-averaged APA 
versus instantaneous relationships is substantial, especially at moderate L^s, but 
is substantially reduced when using linear combinations of reflectances (Fig. 11, 
bottom). In this case, the points corresponding to a same LAI are generally more 
aligned with the best fit line. The result is a drastic improvement in the APAR 
residual error. For the geometries of Figures 10 (bottom) and 11 (bottom), the 
residual error is reduced from 0.064 to 0.029 (planophile case), 0,052 to 0.018 
(spherophile case), and 0.049 to 0.013 (erectophile case). When using NDVI, the 
minimum residual error is obtained for nadir viewing (Figure 10, top), but when 



using linear combinations it is generally preferable to view the target at 60° 
viewing zenith angle, although in the case of a planophile canopy any viewing 
gZ^^“”.^vldes®stailar results (Figure 11, top). Analysis of q“ 

raSometer measurements collected during FIFE have comborated that an 
opHmized linear channel combination scheme is a potentiaily better predictor for 

APAR (Figure 12). 

4, Effect of Broken Clouds on Shortwave Irradiance 
’ In operational satellite algorithms to estimate surface shortwave 
irradiance, top-of-atmosphere reflectances are used to infer pixel-averaged cloud 
reflKtance and, then, cloud transmittance (flux transmittance)^ ^ 

rad^rXted by the surface and clouds is generaliy assumed, and the c ouds 
are considered plane-parallel. This approach suffers a number of limitabons, 
basicaily linked to the anisotropy of the radiance field emerging from . 

B^n clouds, in particular, can affect drastically the spatial distribution of 
surface shortwave irradiance, as illustrated in Figure 13, which contrasts the 
situations on August 4, 1989 a clear day, and August 8, 1989, a parhaUy cloudy day, 

at site 16. « i i ■ 

On August 8, cumulus clouds started to develop around 11:30 local Ume, 
were fully developed by early afternoon, and disappeared around l&OO local hme, 
which mLifested in short periods of alternate high and low irradiance value^ 
Xe high values largely exceeded those observed at the san« hme on Ay s 4. by 
as much as 130 Wm'l between 13:00 and 14:00 local time. The cloud field at 13.32 
Sc” toe, as captured by Dr. Eloranta's lidar installed a few kilometers from the 

pyranoineter, is displayed in Fig. 14. 

Qualitatively, when there are no clouds (August 4), the pyranometer 
measures Ecos6 where E is the solar constant and 8 is the sun zenith angle (we 
neglect in this discussion scattering and absorption by molerales and aerosols)^ 
When there are clouds (August 8), two cases are possible: the clouds obscure or (to 
the sun. If to clouds obscure to sun, the pyranometer measures 

E'<Ecos6, and the shortwave flux is essentially diffuse. If the clouds do not 
obscure the sun, the pyranometer measures E">Ecose, 

because to to direct shortwave flux is superimposed a diffuse flux y 
clouds In other words we expect E"=Ecos9+E' if we assume that to diffuse y 
does not depend on the locaHon of the clouds Figme 16, to^ 

the situation is more complex on August 8 since between 12.00 and 1430 local 
time we estimate E'=270Wm-2 and F'-Ecos0=13OWm-2 or about E'/2. Thus t^ 
excess flux when the sun is not obscured is not equal to the diffuse flux measur 

in the shade of the clouds. 

To go further in this interpretation, consider spherical clouds of radius R 
and, for the sake of simplicity, neglect cloud-cloud interactions. In the ca^ o ^ 
isolated cloud illuminated by a directional beam of irradiance E, Monte C 
calculations show that the average backscattered radiance where x is 

scattering angle, is nearly isotropic for average optical thickness above 16. It 


follows that L=E/4 since the cloud is conservative. Now suppose that we have N 
such clouds covering an area of surface S. The fractional cloud coverage is 
n=N 7 iR^/S. For any sun position (not too low above the horizon to avoid shade 
effects) the N clouds collect a flux NtcR^E. Since they are Lambertian, they scatter 
half of this flux toward the ground. Thus the average diffuse reflectance at the 
surface is N7 cR 2E/2S, that is nE/2. Assume that this diffuse irradiance is 
homogeneous spatially. Depending on whether the sun is obscured or not by a 
cloud, the model predicts that the pyranometer measures nE/2 or nE/2+Ecos0. 
An interesting aspect of the model is that E' does not depend on q, as observed in 
Figure 16. However, taking nE/2=270 Wm-2, Ecos9=860 Wm'2, 0=25° qt 14:00 local 
time gives n=0.57, which is not unreasonable, but too high (see Figure 14). 

So far we have been able to explain qualitatively the measurements, but 
the model remains too simplistic. It appears that the clouds may not be thick 
enough to behave as Lambertian bodies. Therefore we shall write E"=Ecos0+Ed 
with Ed=130 Wm‘2 when the pyranometer is lit by the sun and E'=Ed+Ei when 
the pyranometer is in the shade of a cloud, where Ei is an excess flux that must be 
explained by the cloud transmittance (the clouds are not completely opaque). Let t 
denote the cloud optical thickness. Because of the diffraction peak of the phase 
function (where the light scattered is practically transinitted), we expect an 
effective transnuttance equal to exp(-t/2). Assume that the light scattered (order 2 
and higher orders) is isotropic. We have Ei=Ecosq exp(-t/2) and Ed=nE[l-exp(- 
t/2)]/2. Using Ei=140 Wm-2 and Ecos(25<>)=860 Wm-2, we obtain t=3.6, which is 
reasonable (if the scattering coefficient is 0.5 /km, the cloud geometrical thickness 
is 72m). For the fractional cloud coverage, we obtain n=0.33, which is in better 
agreement with the observations. 

A more specific calculation is necessary to justify the above explanation. 
Using sky photographs and lidar measurements, we can estimate the average 
dimension of the clouds, their average spacing, cloud base altitude and thickness. 
This information was fed into a Monte Carlo code that computed the spatial 
distribution of the cloud transmittance at the surface. Figure 15 gives the results 
obtained at 13:50 local time for cylindrical clouds of optical thickness 12. The cloud 
transmittance reaches values as high as 113% in the region not obscured by the 
clouds. Using a typical clear atmosphere and a cloud-level wind of 4ms' blowing 
from the North (as measured from Pr. Brutsaert's radiosondes), we were able to 
reproduce many features of the pyranometer trace (Figure 13), including the near- 
constant minimum values, the maximum values when the sun was not 
obscured, as well as the fine structure associated with the cloud edges and cloud 
shadowing (Figure 16). 

5. Variability of Downward Shortwave Irradiance 

High resolution (about 1.5 km) maps of hourly and daily surface solar 
irradiance over the FIFE site have been produced for all the IFCs (1 through 5). 
Maps of other useful parameters, namely equivalent cloudiness (one minus the 


ratio of actual and dear sky irradiances), cloud albedo, clear sky irradiance, and 
surface albedo, have also been produced. 

Since cloud-radiation interactions are important to weather and climatic 
prediction, we have examined the role of clouds in modulating insolation 
(Frouin and Gautier, 1990; see Appendix 3). The daily cloud J 

difference between clear and actual sky conditions, reached alxiut 75% of the clear 
sky insolation on some days (e.g., days 180 and 181). Minimum (or maximum) 
values occurred every 2 to 4 days, corresponding roughly to the time scales of 
mid-ladtude atmospheric disturbances. Compared to the effect of changes in 
water vapor amount and aerosol type and amount, the cloud forang 
dominates, by typically one order of magnitude. While developing satellite 
techniques to estimate insolation (or PAR), one should therefore concentrate on 
the role of clouds (determining accurately the governing cloud parameters) rather 
than on the influence of the clear atmosphere. 

Correlation between cloud radiative forcing and fractional cloud coverage 
was also studied using sky camera data and meteorological reports (Br^on et al., 
1990- see Appendix 10). Using the sky camera data, we found linear correlation 
coefficients of 0.62, 0.31, and 0.49 for shortwave, longwave, and net radiative 
forcings, respectively. Using the meteorological station data, we found lower 
correlation coefficients. These rather low values show that cloud amount is not 
the only parameter driving radiative forcing. A governing parametw for cloud 
forcing is liquid water content, but this parameter was not measured. The average 
diurnal cycle of the shortwave, longwave, and net cloud radiative forcings 
exhibited a definite diurnal cycle, with maximum values around local noon for 
the shortwave forcing, and minimum values at the same dme for the longwave 
and net forcings. These cycles, if confirmed by other observations, should ^ taken 
into account when evaluating radiation budgets from heliosynchronous 

satellites. 
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Bk I Performance of the Tanr4 et al. (1985) model to estimate 

solar (shortwave) irradiance at the surface under dear skies. Values are half- 

hourly averages. 



Fie 1 Performance of fhe Morcrette (1984) model to esdinate downward 
temestrial (longwave) irradiance al the surface under dear skies. Values are 

half-hourly averages. 
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Fig. Evolution of the calibration gain of GOES-6 VISSR solar channels. 
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Fie. 5' Scatter plots of calculated versus measured daily insolation (top) and 
PAR (bottom) during FIFE'S second IFC. The measured values are averages ot 
station 5 and 26 data (insolation) and station 5 and 21 data (PAR). 
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Fig 6 Calculated versus measured insolation (top) and PAR (bottom) for 
selected days during FIFE'S second IFC. Satellite estimates are instantaneous, 
whereas measured values are half-hourly averaged. 
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Fig. I ^ Surface solar irradiance on August 4 and 8, ^ Konza p 

Kansas (site 16). (Courtesy of Don Deering and Tom Eck.) 



Fie. Volume imaging Udar backscattering coefficient on August 8, 1%9 at 13:32 
local time, at a few kilometers from site 16 in the Konza prairie, Kansas. (Courtesy o 

Eloranta.) 
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1. INTRODUCTION 

Simple (highly parameterized, analytical) 
radiative transfer models are important for climate 
studies. While retaining the essential physics of the 
problem, they alleviate the computational burden of 
more sophisticated or 'exact' schemes, which is 
particularly useful in general circulation models. 
The simplified assumptions are justified not only 
because of computational efficiency arguments, but 
also because some of the governing parameters are 
not determined or known with sufficient accuracy to 
justify a detailed description of secondary processes. 

Among simplified schemes, the models of Tanr6 
et al., (1979) and Morcrette (1984) have been used 
extensively to compute downward solar (shortwave) 
and terrestrial (longwave) irradiance at the earth's 
surface, respectively (e.g., Gautier et al., 1988; Frouin 
et al., 1988; Gautier and Frouin, 1990). Although 
these models have been verified against exact 
calculations or line-by-line models, they have not 
been validated against in situ measurements. 

During the First ISLSCP Field Experiment (FIFE) 
concomitant radiosonde launches, sky photographs, 
aerosol turbidity and radiation measurements were 
made, offering the opportunity to compare model 
outputs to surface-measured fluxes. The 
comparisons are presented here, but only in clear sky 
conditions because the available data did not permit 
adequate characterization of some of the important 
cloud parameters (e.g., cloud vertical extension). 
Furthermore, the model of Tanr^ et al. (1979) is only 
valid for clear skies. 

First, we briefly describe the radiative transfer 
models and summarize previous comparisons with 
exact calculations of more detailed models; then we 
present the FIFE data used in the comparisons, and 
finally we report on the comparison statistics and 
discuss some implications for climate studies. 


2. RADIATIVE TRANSFER MODELS 

The Simulation of the Satellite Signal in the Solar 
Spectrum (5S) code (Tanr4 et al., 1985; 1986) computes 
the solar radiation backscattered to space by the earth- 
atmosphere system. Given a Lambertian ground 
target and a clear atmosphere, the code estimates the 
target's apparent reflectance by taking into account 
the effects of scattering by molecules and aerosols, 
gaseous absorption, and spatial in homogeneities in 
the surface reflectance. The input parameters, 
namely solar and viewing geometries, atmosphere 
model, surface reflectance, and spectral band, can be 
specified from standard or user-defined conditions. 
In addition to apparent reflectance, the code provides 
gaseous transmittance and irradiance at the surface, 
as well as the various components of the satellite 
signal. Furthermore, exact calculations at selected 
wavelengths allows the user to assess the code's 
accuracy. 

Based on Tanr^ et al., (1979), the top-of-atmosphere 
signal is expressed as a function of the successive 
orders of radiation interactions in the coupled 
surface-atmosphere system. If p is the reflectance of 
the target, and that of its environment, the 
apparent reflectance is written as: 

p* (p, p'. 4>) = tg (p, p') (pa(p, p'. (|>) 

where 9 = arc cos(p) and 6' = arc cos(pl are the solar 

and viewing zenith angles, respectively, ♦ is the 
relative azimuth angle between solar and viewing 
directions, ris the atmospheric optical thickness, tg is 
the gaseous transmittance, p^ is the intrinsic 
atmospheric transmittance, f,/ is the diffuse 



atmospheric transmittance^ and S is the spherical 
albedo of the atmosphere. The first term enclosed in 
the curly brackets represents the contribution of 
photons backscattered to space without surface 
reflection, whereas the second term characterizes 
photons that have sustained one or multiple surface 
reflections. Absorption by atmospheric gases is 
considered as a single multiplicative factor 
dependent on the direct paths sun-to-surface and 
surface-to-sensor. Decoupling absorption and 
scattering processes is justified since, on the one 
hand, ozone absorption is located at altitudes where 
molecules are rarified, and on the other, water vapor 
and carbon dioxide absorption occurs above 850 nm 
where molecular scattering is negligible, and first and 
second orders of aerosol scattering (predominantly 
forward) restitute almost all of the diffuse radiation. 
According to Tanr^ et al., (1986), the error introduced 
by separating the two processes is smaller than one 
percent, except for grazing incidence or observation 
directions (/i,/i'<0.1). 

The atmospheric functions and S are 
approximated by analytical formulas determined 
empirically from exact radiative transfer calculations 
performed for a wide range of model atmospheres. 
Table 1 shows the disparity between 5S and exact 
calculations of the total diffuse atmospheric 
transmittance, , for various zenith angles 

and a wavelength of 850 nm (near the equivalent 
wavelength for the entire solar spectrum). Results 
are presented for two atmospheres, clear and hazy. 
The differences are small, generally less than 1%, but 

reach 1.6% when 9 equals 60” and the atmosphere is 
hazy. 


Tabic 1. Comparison of 5S and cxaci calculations of ihc 
atmospheric loid diffuse transmittance at 850nm. 



e= 

15° 


60° 

Atmosphere 

5S 

Exact 

5S 

Exact 

Clear 

0.973 

0.973 

0.926 

0.924 

Hazy 

0.926 

0.927 

0,819 

0.806 


The gaseous transmittance, tgis computed from 
two exponential random band models, that of Goody 
(1964) for water vapor, and of Malkmus (1967) for 
oxygen, ozone, and carbon dioxide. The spectral 
resolution, 20 cm ^ is sufficient to apply the random 
band models confidently. Table 2 compares the 
gaseous transmittance in the total solar spectrum 
computed with 5S and a well-known code, 
LOWTRAN-6 (Kneizys et al., 1983). For all the 
atmospheres and zenith angles considered, the 
agreement is better than 1.5%, with 5S giving values 
systematically higher. This may be due to the fact 
that LOWTRAN-6, unlike 5S, neglects the influence 
of temperature on the molecular absorption 
coefficient, and approximates molecular line 
absorption by a one parameter band model. 


Table 2. Atmospheric gaseous transmittance in the spectral interval 
250-4000nm computed with 5S and LOWTRAN'6. 


Atmosphere 


15° 


S»60® 

5S 

UOWITIAN^ 

5S 

ljOWmAN6 

Tropical 

0.828 

0.820 

0.793 

0.782 

Midlaittudc Summer 

0.838 

0.830 

0.803 

0.794 

Midlaiitude Winter 

0.876 

0.86S 

0.844 

0.834 

Subarctic Summer 

0.850 

0.841 

0.815 

0.806 

Subarctic Winter 

0.893 

0.880 

0.862 

0.850 


2.2 The wide-band model of Morcrette (1984) 

This model (see also Morcrette et al., 1986) is based 
on a highly parameterized scheme developed for 
GCMs, In a clear atmosphere, upward and 
downward longwave irradiances at altitude z are 
computed from the radiative transfer equation 
developed into: 

(z) = 7C I dv ( [B, - B(0)] t(z, 0; r) +B(z) 

Jo 

■[ ( 2 ) 
Jo 


F (z) = Jt I dv { B(Z) t(z,Z) - B(z) 

Jo 

+ 1 dz'%z')dl{z')t(z'.z;r)) 

^ dz' ( 3 ) 

where z is altitude, Bfz^is the Plank function at 

temperature T(z), t(z, z';r) is the atmospheric 
transmittance between altitudes z and z' evaluated in 
the equivalent direction p=l/r(r is the diffusivily 
factor), B, characterizes the surface emission, 2 is the 
lop-of-atmosphere altitude, v is frequency, and 
superscripts + and - denote upward and downward 
fluxes, respectively. In (2), the effect of a surface 
emissivity e less than unity is taken into account by 
writing Bg as eS(Ts ) + (1-e) F (0) where Tg is the 
surface temperature. 

To evaluate the integral over z on the right-hand 
side of (2) and (3), the atmosphere is divided in N 
layers of arbitrary thickness. Temperature and 
absorber amounts are specified at each level 
separating these layers. In spectral regions where 
atmospheric absorption is strong, the radiative 
energy is exchanged over short distances. Integration 
over the layer adjacent to the surface is therefore 
performed using a 2-point Gaussian quadrature. For 
the other layers, a simple trapezoidal rule is applied. 
When compared with a more accurate 32-point 
Gaussian quadrature for all layers, the simpler 
Integration scheme gives a maximum 4 Wm'^ error 
in the downward flux at the surface (Morcrette and 
Fouquart, 1985), 


To perform the integration over v Rodgers' (1967) 
emissivity approach is followed. Four spectral 
intervals are considered (0-500 cm ’ + 1250-2820 cm 
500-800 cm ', 800-1250 cm ', and 970-1110 cm '), in 
which atmospheric absorption is due mainly to water 
vapor and carbon dioxide (first three intervals) or 
ozone (last interval). Absorption by minor gaseous 
constituents (e.g., methane, nitrous oxide) and 
aerosols is neglected. The contribution of each 
spectral interval to F*(z ) and F'(z ) is evaluated using 
normalized transmissivity functions. Absorption by 
water vapor is treated with the statistical band model 
of Goody (1952) whereas absorption by carbon dioxide 
and ozone is modeled according to Malkmus (1967). 
At high altitudes, these models are modified to 
account for the Voigt profile of the absorption lines. 
The strong line approximation is used for water 
vapor and carbon dioxide, and the weak line 
approximation for ozone. To account for 
temperature and pressure dependence on absorption, 
the Curtis-Godson approximation is applied and 
equivalent absorber amounts are defined accordingly. 
For more details, see Morcrette (1984) and Morcrette 
et al., (1986). 

The wide-band model of Morcrette (1984) has been 
compared with a more detailed, narrow-band model 
(Morcrette and Fouquart, 1985) and with Scott and 
Chedin's (1981) line-by-line model. In the 
comparison, the agreement was within 2-3%. This is 
illustrated in Table 3 from Brton et al., (1990), which 
shows the downward longwave irradiance at the 
surface F(0) obtained with the wide-band and line-by- 
line models for typical atmospheres. The differences 
are small, generally a few Wm ^ Morcrelte's (1984) 
model, however, appears to be more sensitive thart 
the line-by-line model to changes in atmospheric 
temperature and water vapor (the range of values is 
231.3 Wm'^ instead of 225.5 Wm‘^). 

Table 3. Downwclling longwave irradiance at the surface (in 

Wm *) computed using 4 A and Morcrette (1984) models. 

(After Bfdon et al.. 1990.) 


Radiauve Transfer Mode! 


Aunosplicre 

4A 

(line- by dine) 

Morcrette 

(Widc-Band) 

Tro|ilcal 

390.6 

395.1 

Midlatitudc Summer 

341.8 

346.2 

Midlaiitude Winter 

213.1 

2i4.9 

Subarctic Summer 

289.3 

295.7 

Subarctic Winter 

165.1 

163.8 


3. DATA 

The data used to verify the 5S model and the 
wide-band model of Morcrette (1984) comprise 
shortwave and longwave radiative fluxes from 
pyranometers and pyrgeometers, fractional cloud 
cover from a sky camera, vertical profiles of air 
temperature and water vapor mixing ratio from 


radiosondes, and aerosol optical thickness and 
Angstrom exponent from sunphotometers. The 
measurements were made during FIFE’S four 
Intensive Field Campaigns (IFCs), which took place 
in 1987 during the various phases of the vegetation 
seasonal cycle, namely green-up (26 May - 6 June: 
IFC-1), peak greenness (25 June - 15 July: lFC-2), dry 
down (10 August - 21 August: IFC-3), and senescence 
(5 October - 15 October: IFC-4). The experimental site, 
located on the Konza Prairie near Manhattan, 
Kansas, is an approximate 15x15 km square with 
various topological features including plateau, slope 
and creek. Fig. 1 shows the location of the site and, 
within the site, the measuring stations selected in the 
present study. 


''8 




II 

Oil 

^11 

II 

H 

II 


Siie 

iiiiiM 1-70 
= = = n-177 


.muMimo""""' 




A Sky Camera 
O Sunphotometer 
• Pyrenometer» Pyrgeometer 
X Radiosondes 



Fig. I. Map of il>c FIFE site, centered al 39° 04* N and %° 30* W, 
showing the location of the sky camera, radiosondes, and radiation 
stations selected in the study. 


The radiative flux measurements are those made 
at stations 2 and 38 (see Fig. 1). Table 4 compares, for 
downward shortwave and longwave irradiances the 
values obtained at the two stations under virtually 
clear skies (fractional cloud cover less then 5%). The 
r.m.s. differences are small at the scale considered 
(half-hour), representing only 2.6% and 1.7% of the 
average shortwave and longwave irradiances, 
respectively, and the biases are negligible. The 
agreement between the datasets further indicates that 
spatial variability of the clear atmosphere is small 
over the distance between the two stations; therefore, 
in the model computations one can confidently use 
concomitant atmospheric data at nearby (not 
necessarily the same) locations. For the comparisons, 
we averaged radiative flux measurements at the two 
stations. 


Surface albedo was deduced by ratioing upward 
and downward shortwave irradiances. The albedo 
values at the two stations exhibited significant 
changes (not shown here), especially at large solar 
zenith angles. Nevertheless, we averaged those 
values for use in the 5S model further assuming that 
spectral surface reflectance is constant in the entire 
solar spectrum and equal to the broadband surface 




Table 4. Comparison suuislics of downwelling irradiancc measured 
by sialions 2 and 38 during FIFE. Fractional cloud cover, as 
estimated from the sky photographs, is less lhan 5%. Average and 
range values correspond to station 2. 


Paiamelcr 

Nb. 

P«. 

Ave. 

(Wn.J) 

Range* 

(Wm^) 

Corr. 

Coef. 

r m s. Diff. 
(Wnr^) 

Bias 

(Wtn2) 

Shortwave 

Irradiance 

326 

476.4 

294.5 

0.999 

12.4 

2.1 

Longwave 

Irradiance 

339 

355.1 

44.1 

0.991 

6.1 

0.9 


* One standard deviation 


albedo. This procedure is questionable, indeed, but 
somewhat justified because (1) the comparison is 
limited to the domain of validity of the 5S code, that 
is solar zenith angles less than 60°, and (2) the effect 
of surface albedo on downward shortwave irradiance 
is small. 

To select clear sky conditions, we used fractional 
cloud cover data (hourly values) derived from 
photographs taken by an uplooking automatic 
camera. The camera was located within a few 
kilometers of stations 2 and 38 (Fig. 1). In deriving 
factional cloud cover, the percentage of cloudy areas 
in the photographs were corrected by the viewing 
angle of each area relative to the vertical. Since the 
camera operated in the visible, no cloud data were 
available at night. We therefore limited the 
longwave irradiance comparison to daytime. 

The vertical profiles of temperature and water 
vapor mixing ratio were obtained from radiosonde 
measurements of pressure, temperature, and wet 
bulb temperature. The radiosondes were launched 
nearby stations 2 and 38 (Fig. 1). More than 6 
radiosondes were launched on average per day, 
mostly during daytime. The calculated water vapor 
mixing ratio sometimes showed negative values, 
especially during IFC-4, This was the result of 
improper wet bulb temperature measurements below 
273K. We therefore removed all the water vapor 
information corresponding to wet bulb temperatures 
below that value. The data was completed using 
climatological profiles, noting that the water vapor 
mixing ratio cannot be larger than the saturation 
value. 

Aerosol optical thickness and Angstrom exponent 
are those inferred from total atmospheric optical 
thickness measurements at different wavelengths by 
several sunpholometers distributed within the FIFE 
site (Fig. 1). The measurements were made at 
frequent time intervals (typically every half-hour) 
when the sun was not obscured by clouds (the basis of 
the method). In general, the aerosol optical 
thicknesses were small, rarely exceeding 0.4; the 
largest ones occurred during IFC-2. Diurnal 
variability was significant (typically 50% of the 
average values), and so were the differences between 
nearby stations. The Angstrom exponent exhibited 
much more variability, mainly due to errors in its 
determination from already noisy (to about 20%) 
optical thicknesses. Since the aerosol optical 
thicknesses were small and, therefore, the aerosol 


Influence on surface solar irradiance also small, we 
decided to use in the 5S model the average of all the 
concomitant aerosol optical thicknesses and of all the 
Angstrom exponents. 

4. RESULTS AND DISCUSSION 

4.1 Performance of the 5S moild 

Table 5 summarizes the data, except ozone 
amount, used in the 5S model to compute downward 
shortwave irradiance at the surface. The cases 
selected correspond to fractional cloud cover less 
than 5% (sky non obscured) and sun zenith angle less 
than 60°. Furthermore, we only kept the data when 
all the parameters were measured within ± 30 
minutes of the shortwave irradiance measurements 
at stations 2 and 38. In the calculations, ozone 
amount was fixed at 0.34 atm. cm since the effect of 

_2 

ozone variations is small (a few Wm ). The 5% 
limit in cloud cover Is somewhat arbitrary, but 
justified to the extent that no significant changes in 
the comparison statistics were found when lowering 
the 5% limit to 0%. 

As seen in Table 5, the optical thicknesses are 
small (0.23 is typical of a clear atmosphere), the 
maximum values occurring during IFC-2 (0.29 on 
Julian day 192), The Angstrom exponents vary 
moderately around, the value of 1, which is 
characteristic of continental aerosols (70% of dust-like 
component, 29% of water-soluble component, and 
1% of soot component). During IFC-2, however, 
values as low as 0.52 were computed, indicating the 
possible presence of larger, perhaps less absorbing, 
particles on average. The values obtained on Julian 
day 157, 0.13 for optical depth and 1.00 for Angstrom 
exponent, may be compared to those reported by 
Wrigley et al. (1990), 0.12 and 1.17, respectively. The 
agreement is good, giving confidence to our values. 

From the Angstrom exponent, we have an idea of 
the aerosol typ>e. By mixing various percentages of 
dust-like, water soluble, and soot components, we 
found for each case a combination that gives the 
Angstrom exponent obtained from the 
measurements and a single scattering albedo of 0.89. 
This value for the single scattering albedo 
corresponds to the determination of Wrigley et al. 
(1990), who reported values ranging from 0.889 to 
0.891 on Julian day 284. Without other information, 
we used 0.89 for all cases. Our procedure is subjected 
to errors, indeed, all the more as the aerosol 
scattering phase function is another unknown 
(fortunately the forward peak, of interest here, is not 
as variable as the backward peak); but we have to 
keep in mind that the optical thicknesses are small, 
reducing to some extent the uncertainties associated 
with unknown single scattering albedo and phase 
function. 


Table 5. Parameters used in the SS model. 6 Is the sun zenith angle, As is surface albedo, N is 
fractional cloud coverage, Uw is vertically integrated with vapor amounts, t 55 () is aerosol xipttcal 
tJiickncss ^ 350nm, and a is Angstr^kn exponent All these parameters, except 9, were measured in 
situ during FIFE by pyranomeicrs (As ), sky camera (N), radiosondes (Uw). and sunphoiomcters 
{^ 550 - “ )• 


Dale 


e 

(“) 


N 

(%) 


Uw 

(gem 2) 




a 


155.688 

29.34 

0. 18 

155.708 

24.57 

0. 18 

155.729 

20.33 

0.17 

155.750 

17.48 

0.17 

155.771 

16.76 

0.17 

155.792 

18.42 

0 . 17 

155.813 

21.93 

0.17 

155.833 

26.31 

0. 17 

155 . 85 A 

31.52 

0 . 17 

155.875 

37.07 

0. 17 

155.896 

42.81 

0.18 

155.917 

48.65 

0. 18 

156.792 

18.30 

0.17 

156.813 

21.81 

0. 17 

156.833 

26.20 

0.17 

156.854 

31.42 

0.17 

156.875 

36.97 

0. 17 

156.896 

42.71 

0 . 18 

156.917 

48.55 

0. 18 

157.833 

26.09 

0.17 

187.688 

30.09 

0.18 

192.583 

59.35 

0.22 

192.604 

53.48 

0.21 

192.625 

47.62 

0.20 

192.646 

41.80 

0.19 

192.667 

36.08 

0.18 

192.688 

30.57 

0.18 

192.708 

25.68 

0. 18 

192.729 

21.21 

0.17 

192.750 

17.97 

0.17 

192.771 

16.70 

0.17 

229.771 

25.16 

0.17 

229.792 

26.02 

0. 17 

229 . 81.3 

28.45 

0. 17 

229.833 

31.91 

0.17 

229.854 

36.37 

0.17 

229.875 

41.39 

0 . 17 

229.896 

46.78 

0.18 

229.917 

52.41 

0. 18 

229.938 

58. 18 

0. 19 

280.771 

44.54 

0.17 

280.792 

45.66 

0.17 

280.813 

47.79 

0.17 

280.833 

50.64 

0. 17 

280.854 

54.34 

0. 17 

280.875 

58.63 

0. 17 

284.833 

52. 14 

0. 17 


0.0 

1.93 

0. 13 

1.11 

0.0 

2.04 

0.14 

1.07 

0.0 

2.12 

0.13 

1.00 

0.0 

2.21 

0.13 

0.99 

0.0 

2.30 

0.13 

0.99 

0.0 

2.30 

0.12 

1.03 

1.7 

2.29 

0.12 

1.07 

2.9 

2.27 

0.12 

1.07 

3.2 

2. 14 

0.11 

1.05 

1.9 

1.97 

0.11 

1.00 

1.9 

1.81 

0.11 

1.01 

1.5 

1.82 

0.11 

1.02 

3.6 

2. 10 

0.13 

1.30 

3.2 

2.08 

0.13 

1.29 

3.0 

2.07 

0.12 

1.24 

2.7 

2.06 

0.11 

1.12 

1.9 

2.04 

0.10 

1.09 

0.5 

2.03 

0.11 

1.19 

0.0 

2.02 

0.11 

1.35 

0.0 

1.89 

0.13 

1.00 

4.8 

5.05 

0.20 

1.03 

0.0 

4.04 

0.25 

0.56 

0.8 

4.17 

0.26 

0.55 

3.0 

4.37 

0.27 

0.62 

2.4 

4.57 

0.27 

0.62 

2.3 

4.64 

0.27 

0.58 

2.0 

4.59 

0.27 

0.56 

1.1 

4.55 

0.27 

0.61 

0.3 

4.67 

0.29 

0.53 

1.5 

4.79 

0.29 

0.52 

4.5 

4.91 

0.29 

0.62 

0.0 

2.82 

0.10 

0.94 

0.0 

3.01 

0.10 

0.96 

0.0 

3.23 

0.09 

0.99 

0.4 

3.41 

0.09 

1.01 

0.9 

3.47 

0.09 

1.03 

1.5 

3.53 

0.08 

1.05 

2.8 

3.60 

0.08 

1.07 

2.9 

3.68 

0.08 

1.09 

3.2 

3.77 

0.08 

1.12 

0.8 

1.66 

0.07 

0.74 

0.8 

1.72 

0.06 

0.81 

0.8 

1.77 

0.06 

0.88 

0.4 

1.78 

0.06 

0.95 

0.0 

1.80 

0.06 

1.02 

0.0 

1.81 

0.05 

1.10 

0.0 

1.33 

0.08 

1.37 


Figure 2 shows the results of the comparison. On 
a half-hourly time scale, the correlation coefficient is 
above 0.99, the standard deviation is 23.2 Wm'^ 
(about 3% of the average value), and the bias is 13.2 
Wm‘2 (overestimation by the 5S model). Examining 
each IFC separately reveals that overestimation 
occurs for all IFCs, except IFC-2. In the case of IFC-2, 
the 5S model gives lower solar irradiances. As 
indicated above, the Angstrom exponent during IFC- 
2 exhibits relatively smaller values which may be 
associated with less absorbing aerosols (single 
scattering albedo closer to 1). Therefore we may have 
overestimated aerosol absorption in the model, and, 
consequently, underestimated surface shortwave 
irradiance. In fact, because of uncertainties in the 

model input parameters, the 13.2 Wm“^ bias may not 
be significant. 



MEASURED SW IRRADIANCE (Wm“2) 

Fig. 2. Esii mated versus measured downward shortwave irradiance at the 
surface (half hourly averages). Sky conditions are virtually dear 
(fractional cloud cover less ilian 5%). 


4.2 Performance of the wide-band model of 
Morcrette (1984) 

We ran the model of Morcrette (1984) with the 
temperature and water vapor mixing ratio profiles 
determined from the radiosonde observations. 
Ozone mixing-ratio profiles were climatological, and 
the carbon dioxide mixing-ratio was fixed at 350 ppm. 
As for the shortwave irradiance comparison, we 
limited the computations to cases for which the 
fractional cloud coverage was less than 5%. 
Furthermore, we kept only the cases for which the 
radiosonde launches occurred within ± 30 minutes 
from the longwave radiation flux measurements. 
The number of cases selected, 146, is much higher 
than that used to verify the 5S model, 47. 

Figure 3 shows the result of the comparison. On a 
half-hourly time scale, the correlation coefficient, 

0.98, is high and the standard deviation, 13.0 Wm’^ 
is small, representing 3.7% of the average value. A 


-7.4 Wm"^ bias is computed, indicating slight 
underestimation by the model. The agreement is 

good for values above 370 Wm’ , which mostly 
occurred during IFC-2 and IFC-3 (warmer and more 
humid atmosphere). Below that value, the model 
lends to underestimate surface longwave irradiance, 

in some instances by as much as 40 Wm“^, We find, 
as in the comparison with the line-by-line model, 
that Morcretle's (1984) model is more sensitive to 
changes in atmospheric temperature and water 
vapor. The difference, however, is that the best 
results are now obtained for relatively warm and 
humid atmospheres. 



MEASURED LW IRRADIANCE (Wm“^) 

Fig. 3. Same as Fig. 2, but longwave irradiance. 

5. SUMMARY AND CONCLUSIONS 

Using surface and atmosphere data acquired 
during FIFE, we have compared the outputs of two 
radiative transfer models, the 5S model based on 
Tanr4 et al. (1979) and the wide-band model of 
Morcrette (1984), to surface radiation fluxes 
(downward shortwave and longwave, respectively). 
The comparisons were made on a half-hourly time 
scale in virtually clear sky conditions (fractional 
cloud cover less than 5%). The results show a good 
agreement between calculated and measured fluxes, 
namely correlation coefficients above 0.98 and 

standard deviations of 23.2 Wm”^(2.7%) and 13.0 

W m (3.7%) for shortwave and longwave 
irradiances, respectively. The 5S model 

overestimated shortwave irradiance by 13.2 Wm^^ on 
average, whereas the model of Morcrette (1984), 

underestimated longwave irradiance by 7.4 Wm“^ . If 
significant, these biases may compensate partially 
when computing the radiation budget at the surface. 
On longer time scales (e.g., daily or monthly), more 
suitable for climate studies, we expect the standard 
deviations to be reduced substantially. The model of 





MorcreUe (1984) performed belter in relatively warm 
and humid atmospheres and appeared slightly too 
sensitive to atmospheric changes. From our study, 
unfortunately, we cannot conclude on the aspects of 
the radiative transfer models that need 
improvements. This would have required a more 
extensive and complete set of measurements. Using 
the 5S model (or a derived parameterization) with 
typical continental aerosols, a visibility of 23 km, and 
climatological water vapor amounts instead of 
measured atmospheric properties leads to a similar 
standard deviation and an even reduced bias for the 
shortwave irradiance at the surface. This suggests 
that standard and/or climatological, eventually 
seasonally-varying, atmospheric parameters may be 
sufficient to estimate clear sky shortwave irradiance 
at the surface accurately. 

The comparisons made in this paper demonstrate 
the suitability of the 5S model and the model of 
MorcreUe (1984) for surface radiation budget studies, 
at least in clear sky conditions. However, clear skies 
do not represent general atmospheric conditions, and 
validation studies are necessary to assess model 
uncertainties in the presence of clouds. The problem 
is much more complicated, indeed, because clouds 
are highly variable and their properties are difficult 
to measure accurately. Programs such as the 
Atmospheric Radiation Measurement program of 
the Department of Energy, however, may provide the 
appropriate measurements. In any case, it is 
anticipated that the uncertainties will be larger in 
cloudy conditions. Clouds, in particular, are 
generally not plane-parallel, but plane-parallel theory 
is used in Morcrette's (1984) model. Experiments 
may be useful here, for instance, to determine 
suitable corrections. When estimating the surface 
radiation budget from space observations, plane 
parallel approximations are currently used to 
estimate the governing cloud parameters. It may be 
found, because of the uncertainties in deriving 
secondary cloud parameters, that improving some 
aspect of the radiative transfer models are 
unnecessary. 
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ABSTRACT 

We examine the effect of radiometric calibration uncer* 
tainties on satellite- derived net shortwave irradiance at the 
earth’s surface. Using Ref. 5*s satellite technique, we ex- 
press the net shortwave irradiance sensitivity to calibration 
as a function of two basic components that depend on sur- 
face and cloud albedo sensitivities, respectively. The analy- 
sis of these sensitivities for a wide range of atmospheric and 
surface conditions, as well as radiation geometries, shows 
that a 10% uncertainty in the calibration induces up to 
70 Wm“* errors in instantaneous net shortwave irradiance 
(negative when the calibration uncertainty is positive). The 
maximum relative errors are obtained in overcast conditions 
when cloud albedos are high. On a monthly time scale, the 
induced error becomes typically 13 Wm“^ in the tropics 
and 16 Wm”^ in higher latitude regions during summer. 
The error almost vanishes at high latitudes during winter. 
A 10% positive uncertainty in the calibration gives a net 
shortwave irradiance error similar to that induced by the 
3-hour sampling of the ISCCP Project. 

Keywords: Net shortwave irradiance, radiometric cali- 
bration, satellite, surface albedo, cloud albedo, sensitivity, 
climate. 

1. INTRODUCTION 

Although the earth’s climate is expected to undergo 
changes in response to radiative forcing induced by in- 
creased trace gases, it is not clear how these changes will 
take place or how they will affect atmospheric and oceanic 
circulation, precipitation patterns, and cloud cover distri- 
bution. Monitoring radiative fluxes at the earth’s surface 
and the top of the atmosphere over climatic time scales can 
give us some clues as to the nature of these changes. 

Long-term monitoring of radiative fluxes necessitates 
assimilation of observations, data analysis, and radiation 
transfer models for determining fluxes and their eventual 
changes. Since it is anticipated that radiative flux changes 
will be marginally above the noise level, at least over the 
next decade, it is crucial that the observations have long- 
term consistency so that small changes be detected early. 
This can only be achieved if careful calibration of observing 
instruments is performed. 


Successful attempts to estimate radiative fluxes at the 
earth’s surface (e.^., Refs. 5, 4, and 9) have involved high 
space- and time- resolution spectral (visible and infrared) 
radiances from imaging satellite radiometers. While in- 
struments measuring in the thermal infrared are generally 
checked in flight with blackbodies of known temperature, 
instruments measuring in the visible lack on-board calibra- 
tion. These instruments, nonetheless, can be calibrated by 
other means and have been used in validated methods to 
estimate net shortwave irradiance at the earth’s surface. 

It is therefore important to assess the effects of potential 
calibration errors on net shortwave irradiance estimates. 
This is accomplished in the present paper using Ref. 5’s 
model. After briefly reviewing net shortwave irtaiance ac- 
curacy requirements for climate studies, we examine how 
the various satellite-derived model parameters, and, there- 
fore, the net shortwave irradiance, are sensitive to calibra- 
tion. We then discuss the results in view of other potential 
sources of uncertainty, in particular temporal sampling. We 
conclude with a few recommendations for climate monitor- 
ing. 

2. REQUIREMENTS 

The accuracy to which satellite- derived net shortwave 
irradiance at the earth’s surface needs to be determined de- 
pends upon the application of interest and its time scale. 
For instance, when attempting to describe surface net heat 
flux variations associated with the El-Nino/Southem Oscil- 
lation (ENSO) phenomenon over the Pacific Ocean, which 
are on the order of 200 Wm“^, an accuracy of about 
10 Wm”^ over monthly time scales and spatial scales of 
2® of latitude by 10® of longitude is required (TOGA Scien- 
tific Plan). In other oceanic regions where climatic changes 
are smaller over similar time scales, a better accuracy is 
required. 

Over land, shortwave irradiance flux is important for 
validating General Circulation Models (GCM), provid- 
ing GCM boundary conditions, and studying surface- 
atmosphere interactions. The requirements for these stud- 
ies are more stringent than for those over the ocean because 
there Is a need to know the spectral composition and an- 
gxdar distribution dependence of the shortwave irradiance 
in addition to the total amount. If the requirement is lim- 
ited to the amount, the accuracy is also 10 Wm“^, but over 
smaller time (1-10 days) and space (50-100 km) scales. 


Proceedings of the 4th International Colloquium on Spectral Signatures of Objects in Remote Sensing held at Aussois, France, 18-22 January 
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3, NET SHORTWAVE IRRADUNCE 
SENSITIVITY TO CALIBRATION 

31- Derivation of Calibrat ion Eff€ct» 

Our investigation of the sensitivity of satellite-derived 
net shortwave irradiance at the surface to calibration is 
based on computations performed with the model of Ref. 5, 
later refined by Ref. 2. The net shortwave irradiance at the 
earth *8 surface, STVneti is expressed as 

5W„. = (1 - 4.)[(1 - N) + ^■SW^cIoud] (1) 

where Ag is the surface albedo, N is the cloud cover, and 
and SW^oud denote the downward shortwave ir- 
radiance in clear and cloudy conditions, respectively. To 
estimate A,, and STVcloud* following equations 

are used: 

A — Sq g tpi 

' 5o(l-a)(l-a0to.(^o,^) 

= •^o(l “ a)(l -I- A,a') ioti^o) <w(^o) (3) 


*^^cloud — •^o(I Q^)(l — Oc) toi(^o) ^w(^o)« t,u(do)b (^) 

where 5 q is the solar constant (instantaneous), is the 
GOES Visible Infrared Spin Scan Radiometer (VISSR) ra- 
diance in the visible, a and a! are the scattering coefficients 
for direct and diffuse radiation, respectively, is the solar 
zenith angle, $ is the viewing zexiith angle, Ac is the cloud 
albedo, Oc is the cloud absorption, toi is the ozone transmit- 
tance, iu is the water vapor transmittance, and subscripts 
“a” and “b” refer to the water vapor absorption above and 
below clouds. In Eq. (4), Ac is obtained by solving the 
following quadratic equation: 


sol^{6o,e) =“ + (1 - “)(i - + “'(1 - «)(i - 

+ (1 - q)( 1 - a')(l - Ae - Oc)’i4, (5) 

and Oc is fixed at 20% of Ac. To estimate V, a threshold 
technique is employed (for details, sec Ref. 5) or, alterna- 
tively, Ref. Us spatial coherence technique. 

In order to examine how the SW^ci estimates arc sen- 
sitive to the calibration of the GOES VISSR solar channel, 
the only wavelength channel used in Ref. 5’s technique, we 
may linearize Eq. (1) about a reference state: 

( 6 ) 

where g is the calibration gain that allows conversion of 
satellite digital counts into shortwave radiances and the 
overbar denotes the reference state. To simplify Eq. (6), we 
have assumed that the calibration is entirely defined by the 
gain g. For a particular reference state, dSWj^ct/dg thus 
represents the sensitivity to calibration of the net short- 
wave irradiance. Our objective is to estimate dSW^ct/^g 
for a wide range of reference stales (atmospheric and sur- 
face conditions) and assess the effect of time averaging on 
dSW^^Jdg. 

Applying the operator d/dg to both sides of Eq. (l) and 
using Eq. (2) yields the following expression for dSW^ct/dg: 


dSW^ci 

^9 


dAg 

dg 




+ V s^Vdoud 


(l + a'A.) 
dN,. ^ 


clear 


+ “ A,)(SWcloud - SWeXcu) 


dAc 7^(1 - Ag)SWclond 
dg (1 - Ac - Oc) 


This expression is convenient since dSWj^ct/dg appears 
as a linear combination of dAg/dg, dN/dg^ and dAc/dg^ 
the sensitivities to the various satellite- derived parameters 
involved in the 5TVn«t computation. 

Let us now examine dAg/dg^ dNjdg^ and BAd^g sep- 
arately. First, BAgjdg is difficult to estimate when N is 
obtained with the threshold technique because clouds are 
inhomogenous spatially. In this case, depending on the spa- 
tial distribution of vertically integrated liquid water and, 
hence, Ac within the study area of fractional coverage 77, 
BNIBg can take a wide range of values; tracking down the 
most probable values would require an involved statisical 
analysis of a large variety of cloud fields. This is certainly 
beyond the scope of our study. Furthermore, even if we try 
to simplify the problem by taking Ac as a constant, dNjdg 
then becomes a Dirac function of Xmi ~ being 

the threshold, which is obviously not realistic. When the 
spatial coherence technique is used, however, BN jdg can 
be readily obtained. We then have: 


-Tf ^ Mi 


( 8 ) 


where superscripts “s” and “c” refer to the satellite ra- 
diance in the clear and cloudy regions of the study area, 
respectively. Since 


( 9 ) 

where CNmi is the satellite digital count in the reference 
state and / is a function of the digital count only, the ratio 
on the right hand side of Eq. (8) does not depend on g. 
This leads to 


BN 

dg 


0 


(10) 


At first sight, Eq. (10) suggests that Ref. Us method 
should be employed to minimize the influence of calibration 
errors on SWnci estimates. This cannot yet be confirmed 
because BN/Bg ^ 0 might reduce the absolute value of 
dSWjgci/dg as a result of compensations between terms (see 
[7])* Still, in view of the aforementioned difficulty of 
estimating BN/Bg when the threshold technique is used and 
since BN/Bg = 0 when the spatial coherence technique is 
used, which can always be done, we shall ignore hereafter 
BN/Bg in Eq. (7). 

Differentiating Eqs. (2) and (5) 3 rields the following ex- 
pressions for BAg/Bg and &Ac/Bg\ 


dAg ^ g + (1 - g)(l - ^)Ag 
Bg ^(l-g)(l-g') 


= [a(1.2Ac)(l - 1.2Ac) + (1 “ g)(l - a')AcH- 
a^(l_a)(l-a')A?)/ 

(1 - g)(l - g')[l - 2.4A, + (2a' -f 2.88A,)Ac] 

( 12 ) 
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where we have used Eq. 5 and the corresponding equation 
for clear sky conditions. We see that dAgldg is a linear 
function of while dAddg is a more complex function 
of Ac as well as A,. Both sensitivities are influenced by 
atmospheric scattering, but not by osone absorption. We 
also note that QA,fdg / 0 when A, - 0, but dA^dg = 0 
when = 0. 

3,2. Resiilts and Discussion 

In all the results presented and discussed below we as- 
sume, for ease of interpretation, that the calibration uncer- 
tainty is 10% of the gain. This corresponds to the uncer- 
tainty achieved with present post-launch calibration tech- 
niques, be they vicarious using ground targets (Ref. 3) or 
direct using suitably equipped aircraft (Ref. 8). 

First, we examine the sensitivity to calibration of the 
surface albedo (Fig. 1) and cloud albedo (Fig. 2), as ex- 
pressed in Fqs. (11) and (12). Figure 1 indicates that 
dAtjdg depends only slightly on that a 10% uncer- 

tainty in the calibration yields about the same uncertainty 
on i4*. The same type of dependence on and instead 
of A. is obtained for dAc/dg (Fig. 2). When A, becomes 
large (> 0.3), however, dAcfdg increases rapidly at small 
Ac values. Thus, over land, where A, can be large, we ex- 
pect up to 100% SWiict uncertainties for clouds of small Ac 
(i.c., low liquid water content). This effect is quite impor- 
tant, even though 5Wnet is considerably reduced when A, 
becomes large. 



Fig. 1: Sensitivity of surface albedo to calibration gain for 
solar zenith angles of 0 and 60°. 



Jiff, 2: Sensitivity of cloud albedo to calibration gain for 
various surface albedos and solar zenith angles of 0 
and 60°. 


Next we look at the 5Wnet sensitivity to a 10% un- 
certainty in the calibration (Fig. 3). The results arc pre- 
sented as a function of the most important parameters in- 
volved, namely A, (Fig. 3a), N (Fig. 3b), Ac (Fig. 3c), 
and ^0 (Fig. 3d). As indicated in the previous section, the 
SWnti sensitivity is made up of two terms (since we ne- 
glect dNfdg) that depend on dAgldg and dAddg, respec- 
tively. Their contribution is also shown in Fig. 3. For the 
atmospheric and surface conditions selected, we find simi- 
lar sensitivity ranges for all parameters, with a 75 Wm“* 
mATifTniTTi when JV = 1. For instance, the sensitivity 

decreases with from about —50 Wm~^ for = 0° to 
-5 Wm-’ for = 80°. 

Lastly, we analyze the sensitivity of monthly-averaged 
to calibration during the year for various latitudes. 
The results are presented in Fig. 4. As could be expected 
from Fig. 3d, the sensitivity is larger in the tropics most 
of the year (except in the summer). The monthly SW^^t 
sensitivity is only about —13 Wm“* in the tropical region 
compared to an instantaneous SW^net sensitivity on the or- 
der of -70 Wm"^. These results arc to be compared with 
the accuracy obtained in estimating over the ocean 

in the tropics (e.^., Refs, 6 and 7) which is about 18 Wm""^ 
on a daily basis, when assessed by comparing satellite pre- 
dictions with surface measurements. 
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4* Sensitivity of monthly-averaged net shortwave ir- 
Since li the earth’s surface to cahbration g«n 
during the course of the yew for vjmuos 
The surface and atmospheric conditions are those 

of Fig. 3. 


4. 


other sources or uncertainty 


Uncertainties in determining large-scale net shortwave 
irradiance at the surface have other sources than calibra- 
tion. For example, they may oripnate from applying an im- 
perfect model to the satellite radiances and from inadequa 
cloud sampling. Most of these uncertainties have a rimdom 
component such that when averages are computed, the un- 
certainties are reduced. Calibration uncertainty however, 
generaUy introduces a bia. (calibration fluctuations have 
fong time periods) that does not diminish through averag- 
ing procedures. While it is not our intention ,n the present 
paper to make an exhaustive study of the uncertainties that 
might affect SW„., computations, it is interesting to com^ 
pare the uncertainty due to calibration to that resulting 
Lm temporal sampHng, for instance. In order to provi 
orders of magnitude, we write SW„.t as a f"“ction 
cloud parameter, CL, which represents the integrated ef- 
fects of clouds on the shortwave irradiance: 

SW.a = (1 - A.)5Wd.„(l - CL) (13) 

If Ah denotes the sampUng interval, the sensitivity of 
SWnet to Ah can be expressed as: 

95Wn,t _ - ij)(i -nL)-^^SWrU^{l- As) 

where we neglected dA./d{Ah). To compute 9S Wn,./aAh 
we have to assume a variation of CL with time. In he 
tropics, for instance, a typical diurnal cycle exists that can 
be parameterized as 

CL(h) = 0.5 -0.2 cos j (15) 



ria 5* Sensitivity of monthly-averaged net shortwave irra- 
iW at^he earth’s' surface to 
A typical diurnal variation is assumed for the cloud 




When Eq. 15 is used and SW„t a,veraged over one 
month, we obtain the results presented in Fig. 5. We see 
that a5W„rt/5Afc«aches -12 Wm when Ah = 5 and 
that the largest contribution comes form samphng SWd^. 
If we change Ah from 1 to 3 hours, for example, SVVnet will 
change by: 


ASW„et = 5W„t(Ah = 3) - 5W„t(Ah = 1) 




V dAh /Ah=5 

\ dAh ) 

-15Wm“* 



This bias is similar in magnitude to that introduced by 
a calibration uncertainty of 10%. It is interesting to note 
that for the conditions considered, the sampling adopted 
in the International Satellite Cloud Climatology Pr^ram 
(ISCCP, ».e., 3 hours) introduces a bias of -15 Wm . 

5. SUMMARY AND RECOMMENDATIONS 

We have performed a study of the sensitivity of the net 
shortwave irradiance at the earth’s surface to calibration. 
This study has showed that a uncertainty in the calibration 
leads to uncertainties up to -70 Wm"* for instantan^us 
values of net shortwave irradiance in overcast conditions 
and up to about -15 Wm"* for monthly averaged values. 
In the tropical regions, the monthly averaged uncertwnty 
is —13 Wm"*. This uncertainty appears as a bias with a 
sign opposite to that of the calibration bias, i.e., a positive 
bias in calibration induces a negative bias on net short- 
wave irradiance. We have compared the uncertainty due to 
calibration to 'that introduced by temporal sampling. We 
have found that the ISCCP sampling (3-hour) introduces a 
bias of typically -15 Wm"* in the tropics. This indicates 
that, if the calibration bias is negative, the biases intro- 
duced by the calibration could compensate that introduced 
by temporal sampling. The amount of compensation could 


where ho is local noon. 
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be partly assessed by comparing monthly- averaged estima- 
tions of net shortwave irradiance from satellite observations 
to those obtained from pyranometer measurements. 

Considering the restdit just summarized, it is very im- 
portant to improve the calibration of spectral radiances 
used to compute net shortwave irradiance at the earth’s 
surface, such as those from imaging radiometers. This cali- 
bration can be performed on-board future instruments, but 
for present instruments it must be done using a combina- 
tion of complementary techniques. Calibration of the visi- 
ble radiances is important for ISCCP data sets since they 
are, at present, only normalized to one particular visible 
sensor, which is itself poorly calibrated. We reco mm end 
that a calibration strategy be established with the objec- 
tives of providing both long-term calibration changes and 
the best possible absolute instantaneous calibration. Such 
a strategy could involve a continuous monitoring of cali- 
bration targets and a complete calibration experiment that 
provides the necessary measurements to apply all possible 
calibration methods. 
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VARIABILITY OF PHOTOSYNTHETICALLY AVAILABLE AND TOTAL SOLAR 
IRRADIANCE AT THE SURFACE DURING HFE: A SATELLITE DESCRIPTION 


Robert Frouin and Catherine Gautier 
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1. INTRODUCTION 

Sunlight drives the energy, water, and biochemical 
cycles of the eanh surface/atmosphere system. The flux of 
solar energy through the various spheres of the system, 
namely the stratosphere, troposphere, and hydrolithosphere, 
is understood generally, but specific rates and underlying 
processes are not well-known. This prevents our accurate 
prediction of the system's response to natural and 
anthropogenic changes, including changes inherent to the 
system's internal dynamics. 

Over land surfaces, this study's domain of interest, 
insolation largely determines the surface temperature and the 
rate of evapocranspiration, with imponant consequences on 
air-surface interactions and the global hydrologic cycle. 
Precipitation, for instance, is sensitive to changes in 
evapocranspiration, at least on a regional scale. On the other 
hand, solar radiation at the wavelengths hetw'een 0.35 and 
0.7 jjm is used by plants in the process so fundamental for 
life: oxygenic photosynthesis. Photosynthetically active 
radiation (PAR) governs primary production, the rate of 
carbon fixed by terrestrial vegetation, and knowing the 
geographical location of the fixed carbon, as well as its rates 
and forms of release, is important in studying the increase in 
atmospheric CO 2 and trace gas concentrations. 

Various methods, based on visible and near-infrared 
observations from meteorological satellites, have been 
proposed to estimate insolation at the eanh's surface (e.g., 
Tarpley, 1979; Gautier et ai., 1980; Moser and Rashke, 
1984; Pinker and Ewing, 1985; and Dedieu et al., 1987^. 
These methods offer the opportunity to map insolation 
accurately (within 10% on a daily rime scale) over large areas 
and, eventually, the entire globe (global scale satellite 
datasets are now becoming available). The perspective of a 
global, long-term climatology of surface insolation, which 
can only be established using satellites, should help in the 
investigation of the climate system’s sensitivity to surface 
processes by providing validation capabilities and boundary 
conditions. 

While numerous studies have been devoted to 
estimating insolation from satellite dafa, no satellite-based 
method has been proposed for PAR. This is because, for 
many, applications involving small space and time scales, 
PAR can be measured directly. Funhermore, it has often 
proved satisfactory to take PAR as a more or less constant 
fraction of insolation. The relationship between PAR and 
insolation, however, depends on atmospheric conditions and 
radiation geometry (e.g., Baker and Frouin, 1987). Clouds, 
in panicular, which do not absorb at P.AR wavelengths but 
do absorb substantially in the near- infrared, increase the ratio 
of PAR and insolation. Therefore, deducing PAR from a 
satellite estimate of insolation (also subject to uncertainty) is 
not optimum. Moreover, since meteorological satellites 


(except METEOSAT) carry instruments that measure in 
spectral channels resembling more the PAR wavelength 
range than the entire solar spectrum, it appears appropriate to 
use these data to esdmaie PAR directly. 

In the present study, we apply the satellite method of 
Gautier et al. (1980) to GOES-6 Visible and Infrared Spin 
Scan Radiometer (VISSR) data acquired during the First 
ISLSCP Field Experiment (FIFE). The objective is to 
assess whether the method can provide quantitative 
information on the temporal variability of PAR and 
insolation during the experiment. First, we describe the 
salient features of the Gautier et al. (1980) method and 
indicate how the method is modified to yield PAR. The 
problem, in this case, is simplified since the VISSR solar 
channel is mostly sensitive to photosynthesis wavelengths, 
and cloud absorption vanishes in the radiative transfer 
model. Subsequently, wc discuss the calibration of the 
VISSR solar channels, and then we compare the satellite^ 
estimates to the in situ measurements. The suitability of 
deriving PAR from insolation is evaluated, with implications 
for global studies of photosynihcdc activity, since large-scale 
satellite climatologies of insolation are now becoming 
available. Finally, the imponance of clouds in modulating 
insolation and PAR is analyzed, and we conclude with a 
perspective on future work. 


2. METHOD 

The method used to compute surface insolation, 
originally developed by Gautier et al. (1980), is based on 
simple, physical modeling of the most imponant radiative 
processes occurring within the atmosphere, namely 
scattering and absorption by molecules, clouds, and 
aerosols. Since variability of surface insolation results 
pnmanly from changes in solar zenith angle and cloudiness, 
the method focuses on determining the etfect ot clouds on 
surface insolation (solar zenith angle can be computed 
accurately from simple formulas). The method accomplishes 
this by computing cloud albedo, the governing cloud 
parameter, from GOES VISSR measurements in the visible. 
The repeat coverage of the GOES VISSR data (one 
observ'ation every 30 minutes) allows one to adequately 
sample the diurnal cloud variability, a strength ot the 
method. 

In the computational procedure, we first determine 
from a time series of satellite images (typically 15 days) the 
minimum brightness count at each observation time during 
the day. This minimum count defines a threshold (taken a 
few counts higher) that is used to classify each GOES 
VISSR pixel as clear or cloudy. The procedure, indeed. 


98 



does not allow one to determine whether the pixel is panially 
contaminated by clouds or not. We utilize full resolution 
data, however, which minimizes the error introduced by not 
modeling the resulting effect on insolation. Once the pixel’s 
nature (clear or cloudy) has been determined, we apply clear 
and cloudy sky radiative transfer models accordingly. 

In clear sky conditions, surface insolation is expressed 
as: 

-^0(77) ^ cos 6 exp ( - C^/cos 0)/{\ - C^A^) 

X exp [ - cos 0 } cos fl ^ 1 (1) 

where So is the solar constant, rfr^ is the ratio of actual to 
mean eanh-sun distance, 9 is solar zenith angle, Uo and 
are ozone and water vapor amounts, respectively, As is 
surface albedo, and bo. C j , and are 

coefficients (C/ and C 2 depend on the type and concentration 
of aerosols). The term I-C^As accounts for photons that 
have sustained multiple surface reflections. Equation (1) 
differs from that of the original model, but not in essence. 
Ozone and water vapor amounts are specified from 
climatology and As is obtained by solving the following 
equation. 

a + (l- a)(l- Gj)(l- ^2) 

where Asoj is the albedo measured at the satellite (the surface 
is assumed to reflect solar radiation isotropically), CNmm is 

the minimum brightness count, a and a/ are direct and 

diffuse reflection coefficients, respectively, and 
characterizes ozone absorption. Equation (2) simply states 
[hat Asai is the sum of an atmospheric component (photons 
reflected back to space without surface reflection), and the 
signal reflected by the surface and diffusely transmitted to 
space. 

In cloudy sky conditions, the clear sky formulation is 
modified to account for reflection and absorption by clouds 
which are assumed to occur in one layer. Cloudy sky 
insolation is therefore given by: 

^ c ^ ~ c ~~ (3) 

where is cloud albedo and is cloud absorption (taken as 
a constant fraction of Ac). In Eq. (3) we neglect the effect of 
multiple reflections between the cloud and the surface. This 
effect, of order fc Ac is generally small (a few Wm‘^). 
Cloud albedo is obtained by solving the following quadratic 
equation: 

= a + (1- a)(l - o^)(l- a^yAc 

+ A, - a^ni- a,)A, 

where Asai is the top-of-aimospherc albedo, assuming that 
clouds reflect solar radiation isotropically. This equation, in 
fact, gives Ac in the GOES VISSR solar channel (mostly 
wavelengths in the visible). We assume, however, that 
takes the same value in the total solar spectrum, which is 
reasonably well justified (sec, for instance, Welch ec al., 
1980). 

The above formalism, developed for insolation, only 
requires slight modifications to become applicable to PAR. 
In this case, we eliminate Qc in the cloudy sky equations 
(clouds do not absorb in the visible), replace So by the 
extraterrestrial solar irradiance in the PAR spectral interval, 
and modify the coefficients of the clear model so that they 
represent the PAR wavelengths. For these coefficients w^e 
take the values given by Frouin et al. (1989). Since Qc does 
not need to be parameterized as in the insolation case, and 
since the VISSR solar channel mostly captures radiation in 
the visible (no narrow-band to broad-band transformation is 
necessary), we expect the model to perform better in the 
PAR case. 


3. INSOLATION AND PAR COMPUTATIONS 

We focus, in this preliminary report, on FIFE'S second 
intensive field campaign (EFC), which took place at the end 
of June and the beginning of July 1987 (peak greenness of 
vegetation) on the Konza Prairie near Manhattan, Kansas. 
Figure I locates the experimental site and, within the site, the 
radiative flux stations selected for validadon purposes. 

To compute insolation and PAR over the FIFE area, 
GOES-6 VISSR visible and near-infrared (solar channel) 
data were acquired at full resolution (0.9 km at nadir) every 
half-hour during daytime. The data, 8-bit coded, navigated, 
but uncalibrated numerical counts, were made available by 
the Space Science and Engineering Center of the University 
of Wisconsin, Madison. In a preprocessing stage, we 
calibrated and checked the data for navigational errors and 
quality. 

Calibration was achieved by regularly applying (twice 
a month) the indirect calibration method described in Frouin 
and Gautier (1987). This method, which uses space and the 
White Sands Monument area in New Mexico as calibration 
targets, computes the radiance measured at satellite altitude 
using a fairly accurate radiative transfer model, the model of 
Tanre et al. (1986). The relevant atmospheric characteristics 
are estimated from cither climatological data (ozone amount, 
aerosol size-frequency distribution, and refractive index) or 
observations at nearby meteorological stations (water vapor 
amount and aerosol loading). The method’s theoretical 
accuracy is 8%, but it might even be better, as comparisons 
with a more direct method involving high-flying aircraft have 
suggested (Whitlock et al., 1987). Figure 2 shows the 
calibration gain (relates count-squared to reflectance) 
obtained during 1987. In the figure, each point corresponds 
to a daily average performed with three White Sands 
observations near local noon and one observation of space. 
The venical bars represent one standard deviation (around 
the daily average) of individual calibrations. Somewhat 
erratic fluctuations are obsen.'ed from one gain estimate to 
the next, in some instances reaching 10%. Changes of such 
amplitude, however, are at the limit of the uncenainty level. 
In addition to the relatively high-frequency fluctuations, the 
gain exhibits a low-frequency variation, with minimum 
values in April, May, June, and July. Prior to 1987, similar 
yearly variations have been detected by other indirect 
calibration methods (sec Whitlock et al. 1987), which 
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Fig. 1: Map of the FIFE site, centered at 39^ 04* N and 
96*^ 30’ W, showing the location of the insolation and PAR 
stations selected in the study. 
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Fig. 2: Calibration gain of the GOES-6 VISSR solar 

channel during FIFE. The prelaunch value is 0.154 x I0“* 
reflectance/count^. 


suggest that the observed feature is actual and not the result 
of a bias in the method. From July to December 1987, the 
gain increased by approximately 36%. Compared to the 
prelaunch gain, the gain at the end of 1987 reads higher by 
64%. 

Navigation was checked by comparing the computed 
earth coordinates of reference points on the surface, in 
panicular small lakes (that appear very black on the visible 
images), to the actual coordinates. The agreement was 
generally good, within a few kilometers. In extracting the 
FIFE area from the satellite images, we therefore shifted the 
pixel lines and samples accordingly. Note, however, that 
navigational errors of a few kilometers, although they 
change the actual radiation geometry, do not affect the 
insolation and PAR results significantly (errors in the' 
parameterizadon of atmospheric process largely dominate). 

After the preprocessing stage, we produced 
instantaneous and daily e<':mates of insolation and PAR by 
applying the procedure described in section 2. In the 
calculations, we used the continental aerosol model of the 
International Radiation Commission (WCP, 1983) and fixed 
the atmospheric visibility (characterizes aerosol turbidity) at 
23 km. Climatological values of 3.0 gcm-^ and 
0.3 1 atm-cm were taken as reprcsentarive of the water vapor 
and ozone amounts during the EFC, respectively. Since the 
surface albedo of vegetated surfaces differs substantially In 
the visible and total solar spectrum (vegetadon reflects more 
in the near-infrared), we increased the VISSR-derived 
surface albedo by a typical factor of 1.64 in the insolation 
calculation (Eq. 1). The daily averages were obtained by 
simple trapezoidal integration, which is sufficient since the 
satellite observations are well-sampled in time 
(approximately 24 observations at 30 minute intervals during 
daytime). Such a simple scheme, however, biases low the 
daily integrals when the sampling frequency is not as high 
(case of missing data), especially in clear or mostly dear sky 
conditions. 


4 , COMPARISON WITH IN SITU 
MEASUREMENTS 

For selected days during the IFC Fig. 3 shows the 
diurnal variation of measured insolation at stations 5 and 26 
and the corresponding satellite estimates at panicular times 
during the day. In the figure, the surface values are half- 
hourly averages and the satellite estimates are spatial 
averages over the FIFE area. Julian days 176, 177, and 178 
are mostly clear, whereas the others are mostly cloudy. 
(Note that the Julian days are depicted w-ith respect to GMT 
and not local time). We notice slight differences between the 
values measured by the two stations, but not a systematic 
bias. The discrepancy is more pronounced in cloudy 
conditions (e.g., day 188), which plausibly results from 
spatial cloud variability (stations 5 and 26 are 15 km apart). 
In general the satellite estimates are in agreement with the 
measured values; they describe well the actual diurnal cycle 
Even rapid insolation changes, for instance at the end of day 
177, are detected. During day 189, the satellite estimates did 
not see the early afternoon variations, but the tw'o stations 
measured a significandy different incoming solar flux. 

The same type of results are obtained for PAR 
(Fig. 4). Satellite estimates are in good agreement with 
surface values, and reproduce well the diurnal cycle. The 
surface PAR data, however, are noisier than the surt'ace 
insolation data. Station 5, for instance, which measured 




Fig. 3: Calculated and measured insolation for selected days 
during FIFE’S second IFC. Satellite estimates are 
Instantaneous, whereas measured values are half-hourly 
averaged. 
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Fig. 4: Same as Fig. 3 but for PAR. 
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both PAR and insolation, reports anomalously low PAR 
values on day 178. Examining Figs. 3 and 4, we note the 
similar diurnal variation of PAR and insolation, which 
indicates, as we know, that the two quantities are highly 
correlated. 

Figures 5 and 6 show, respectively, the temp<^l 
variation of daily insolation and PAR during the entire LTC 
(21 days total). The agreement is generally within 15- 
20 Wm‘^ for insolation and 5-10 Wm'^ for PAR. Day 
182, however, stands out. The method, in that case, did not 
detect much cloudiness. In fact, the clouds during day 182 
exhibited brightness count values slightly less than the 
dear/cloudy threshold and, therefore, w'ere not seen. This 
is, indeed, the disadvantage of using a threshold to separate 
clear and cloudy pixels; notwithstanding, with exception of 
that day, the clouds were distinguished properly. If we now 
examine the scatter plots of estimated versus measured daily 
insolation and PAR (Figs. 7 and 8), we do not notice a bias 
in the satellite estimates. In cloudy conditions (the lowest 
values on the plots), the PAR estimates are closer to the line 
of slope 1 than the insolation estimates, which indicates a 
better performance of the satellite method in the case of 
PAR. This is somewhat expected, as mentioned in the 
introduction, since the modeling of the cloud processes is 
simplified for PAR. The atmospheric conditions, however, 
were rather clear during FIFE'S second IFC, and more 
comparisons in the presence of clouds are necessary to 
quantify the improvement. 



JULIAN DAY (1987) 

Fig. 5: Calculated and measured daily insolation during 
FIFE'S second IFC. 



JULIAN DAY (1987) 
Fig. 6: Same as Fig. 5 but for PAR. 


Table 1 summarizes the statistics of half-hourly and 
daily comparisons for both insolation and PAR. The 
statistics are drastically improved when passing from a half- 
hourly to a daily time scale. Even though the correlation 
coefficient is similar, yet above 0.95, the standard error of 
estimate decreases from 86.7 to 21.6 Wm'^ in the case of 
insolation and from 41.9 to 8.2 Wm'^ in the case of PAR. 
The bias always remains small, not exceeding a few Wm'^ 
on a half-hourly time scale, and reducing to about 1 Wm’^ 
on a daily time scale. For PAR, however, the overall 
statistics arc better. The standard error of estimate for daily 
comparisons, in panicular, represents only 6.5% of the 
average measured value in the case of PAR, but reaches 
almost 9% in the case of insolation. Note finally, that the 
results obtained for insolation arc similar, in terms of 
comparison statistics, to those reported by Gautier et al. 
(1980), Diak and Gautier (1983). and Raphael and Hay 
(1984), who applied the satellite method over varied 
terrestrial sites. 
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Fig 7: Scatter plot of calculated versus measured daily 
insolation during FIFE'S second IFC, The measured values 
ai'e averages of station 5 and 26 data. 
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Fig. 8: Same as Fig. 7 but for PAR. The measured values 
are averages of station 5 and 2 1 data. 


Table 1. Comparison siatisrics of calculated and measured 
insolation and PAR. 



Insolation 


PAR 

Ptirameier 

Half- 

Hourly 

Daily 

Half- 

Hourly 

Daily 

Correlaijon Coefficient 

.961 

.964 

.950 

.968 

Standard Error (Wm‘^) 

86.7 

21.6 

41.9 

8.2 

Bias(Wm‘^) 

-5.3 

-L5 

4.0 

1.4 

Number of Points 

469 

21 

471 

21 


5. DISCUSSION 

We have shown in the previous section that the method 
of Gautier et al. (1980) for insolation and its modified 
version for PAR provide reasonably accurate Insolation 
estimates, especially on a daily time scale. On a half-hourly 
lime scale, the assessment of accuracy is inherently more 
difficult to make since the satellite estimates are quasi- 
instantaneous whereas the in situ measurements are time 
averages. Furthermore, it is not clear whether a spatial 
average over 15x15 km at a given time is represented 
adequately by a half-hour average at a specific location. The 
discrepancies, however, are substantially reduced on daily 
imeerals. The accuracy of the daily insolation estimates is 
sufficient, or at least useful, for large-scale surface- 
atmosphere interaction studies (10 Wm ^ is typically 
required for climate studies) and that of PAR adequate tor 
global primary production computations. In the latter case, 
The major uncertainty comes from esdmating the fraction ot 
PAR that is intercepted and absorbed by the vegetation 
canopy. 

Insolation and not PAR, however, is being produced 
from global satellite datasets, such as those compiled within 
the framework of the International Satellite Cloud 
Climatology Project (ISCCP; Schiffer and Rossow, 1983). 
Since PAR and insolation are highly correlated (see, for 
instance. Figs. 3 and 4), it is interesting to compare the 
accuracy with which one can deduce PAR from a satellite 
estimate of insolation to that obtained when PAR is directly 
estimated from the satellite data (the method described in 
section 2). Figure 9 shows the ratio of daily PAR and 
insolation at station 5. We did not use the data of stations 21 
(PAR) and 26 (insolation) because, even though these 
stations were close, they were not at the same location, and, 
therefore, spatial cloud variability may introduce errors in the 
ratio (in fact an anomalous value was obtained for day 185). 
We notice that the ratio is fairly constant around the average 
value of 0.448. We do not detect at this time scale, 
however, a tendency for the ratio to increase when average 
cloudiness increases, as theory predicts in cloudy 
conditions. 

Using the average value of 0.448 for the ratio of daily 
PAR and insolation with the satellite estimates of daily 
insolation presented in section 4, we deduced PAR values 
that w'c compared to in situ measurements. The comparison 
statistics revealed the same correlation coefficient as in 
Table 1 (0.969) and a slightly larger standard error of 
estimate (lO.O instead of 8.2 Wm ^) and bias (2.4 instead of 



Fig. 9: Ratio of daily PAR and insolation during FIFE s 
second EFC. 
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1.4 Wm'^). Although based on 21 points, these figures 
show convincingly that satellite estimates of [isolation yie 
useful PAR estimates on a daily time scale. Keep m mind, 
however, that the ratio of PAR to insolation may change 
significantly depending on the type of atmosphere, n a 
moist and humid clear environment, it may reach values 
above 0.50. Large variations in the ratio, however are not 
likelv to occur at a specific site during a few week period 
(see Fig. 9). Yet in calculating PAR from insolation we 
determined in situ the ratio of PAR and insolation, which 
minimizes the uncenainty in this ratio. On the other hand, 
adapting the insolation model to PAR (and this may apply to 
other satellite methods) is rather simple and eliminates not 
onlv the uncenaimv in the PAR/insolation ratio, but also the 
uncertainty in cloud absorption and narrow-band to broad- 
band conversion, which atfects insolation estimates. 


Finally, since cloud radiarion interacrions are imponant 
10 weather and climatic prediction, it is worthwhile 
e.xamining the role of clouds in modulating insolation. 
Following Charlock and Ramanathan (1985), we define the 
cloud forcing on insolation (PAR) as the difference between 
insolation in clear and actual sky conditions. Figure 10 
shows the daily cloud forcing during the IFC for both 
insolation and PAR. In the calculation we used the satellite 
estimates and computed clear sky irradiances from the 
sky models. The cloud forcing reaches about 75% of the 
clear sky insolation on days 180 and 181. Minimum (or 
maximum) values occur every 2 to 4 days, which 
correspond roughly to the time scale of mid-latitude 
atmospheric perturbances. Compared to the effect or 
changes in water vapor amount and aerosol type and amount 
(Fig. 1 1), the cloud forcing largely dominates (by typically 
one order of magnitude). This indicates, in other respects, 
ihat while developing satellite techniques to estimate 
insolation (or PAR), one should concentrate on the role ot 
clouds (and, therefore, determine accurately the governing 
cloud parameters) rather than on the influence of the clear 
atmosphere. 
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Fig. 1 1: Forcing on daily insolation due to changes in water 
vapor amount (a), and aerosol type and loading (b). 



Fig. 10: Cloud forcing on daily insolation and PAR during 
FIFE'S second IFC as deduced from the satellite estimates. 


5. SUMMARY AND CONCLUSIONS 

In the first pan of the study, we have indicated how ihe 
atellite technique of Gautier et al. (1980) can ^ modified to 
ield PAR as well as insolation estimates. In the case ot 
»AR, cloud absorption vanishes in the model equations and 
tear* sky PAR is computed using the coefficients giyeri by 
Touin et al. (1989). The technique w^then applied to 
30ES-6 VISSR data acquired during FIFE s second ILL. 
lie computational procedure included checks of the satellite 
lavigation and of the data quality, as well as calibration of 
he VISSR visible channel. The calibration, achieved using 
pace and White Sands, New Mexico as calibration targets, 
evcaled an imponant sensor degradation, reaching 64% m 
December 1987. Compared to in situ measurements, the 
latellite estimates of insolation and PAR exhibited small 
)iases, correlation coefficients above 0.95, and standard 

errors of estimate of 21.6 Wm'^ (9%) and 8.2 Wm‘ 

'6 5%), respectively. These errors are acceptable for climate 
Judies. The overall statistics were found to be better tor 
PAR which is expected since, in that case, cloud absorption 
ioes’not need to be parameterized and no na^ow-band to 
broad-band conversion is necessary. 

also described well the diurnal variability of PAR and 
insolation. Using the ratio of PAR and insolation 
determined in situ with the satellite estimates of insolation, 
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we obtained daily PAR values with a slightly better r.m.s. 
accuracy than when PAR was directly estimated from the 
satellite data. This indicates that the large-scale satellite 
climatologies of insdladon now being produced within the 
framework of ISCCP can be convened simply into useful 
climatologies of PAR for global studies of photosynthedc 
activity, although a more direct use of the satellite data would 
cenainly be optimum. Examining the cloud forcing on 
insolation and PAR, we found that over the time period 
considered, clouds strongly modulate insolation with a 2-4 
day time scale, variations in the characteristic of the 
atmosphere yielding irradiance changes of typically one 
order of magnitude less. This is not surprising, indeed, but 
the situadon might be different when considering insoladon 
variability over longer rime periods (e.g., one year). 

The above results, however, only concern a 3 week 
period. Three other IFCs were conducted during RFE in 
1987, and they correspond to different seasons and states of 
the vegetation (surface albedo). In addirion, only data from 
a few surface stations were analyzed in the present study, 
and many more stations provided radiation data. Our next 
step, therefore, is to extend the conclusions drawn in this 
study to the other IFCs, in pardcular regarding accuracy of 
insolation and PAR methods, and contrast the surface flux 
variability (including spatial variability within the FIFE site) 
during each IFC. More work will also be devoted to the 
surface albedo, a parameter that needs to be specified as a 
boundary condition in general circulation models. We did 
not elaborate on this parameter in the present study because 
of its small influence on insolation, but this is not the case, 
indeed, when the net radiation budget is considered. 
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Abstract 


Current satellite algorithms to estimate photosynthetically available radiation 
(PAR) at the earth's surface are reviewed. PAR is deduced either from an 
insolation estimate or obtained directly from top-of-atmosphere solar radiances. 
The characteristics of both approaches are contrasted and typical results are 
presented. The inaccuracies reported, about 10% and 6% on daily and monthly 
time scales, respectively, are useful to model oceanic and terrestrial primary 
productivity. At those time scales variability due to clouds in the ratio of PAR 
and insolation is reduced, making it possible to deduce PAR directly from 
insolation climatologies (satellite or Qther) that are currently available or being 
produced. Improvements, however, are needed in conditions of broken 
cloudiness and over ice /snow. If not addressed properly, calibration /validation 
issues may prevent quantitative use of the PAR estimates in studies of climatic 
change. The prospects are good for an accurate, long-term climatology of PAR 
over the globe. 


Introduction 


Solar radiation reaching the earth's surface in the wavelength range 0.35-0.7 pm 
is used by aquatic and terrestrial plants in photosynthesis. Called 
photosynthetically available radiation (PAR), it governs primary production, the 
rate of carbon fixed by the plants. Knowing the geographical location and 
temporal variability of the fixed carbon and its forms of release is important in 


assessing the climaUc impact of anthropogenic changes such as the destruction of 
major vegetation systems or the increase in atmospheric carbon dioxide. PAR is 

defined by 


0.7 

PAR(Wm"^)= jKWdX 

0.35 


( 1 ) 


where l(X) is the downward spectral irradiance at wavelength X. Since 
photosystem processes are quantum reactions, it is useful to consider the 
equation 


PARCquantam ^s ') = 


0.7 

— fXI(^)dX 
he •' 


0.35 


(2) 


where h is Plank's constant and c is the velocity of light in vacuum. Eqs. (1) and 
(2) indicate that PAR depends on the spectral interval considered which, for 
operational constraints, may sometimes differ from 0.35-0.7 pm. 


Fig. 1 shows how primary production varies as a function of PAR over land (Fig. 
la) and ocean (Fig. lb). The land case corresponds to typical, live, horizontal 
leaves (Sellers, 1985, Fig. 13a) and the ocean case to a 20°C, homogeneous water 
body (calculations were performed with the model of Morel, 1988). Over land, 
primary production increases rather linearly with PAR, the slope of variation 
depending on leaf area index (higher slope as leaf area index increases). The 
relationship, however, is affected little by leaf area index for leaf area indices 
above 4. Over the ocean, by contrast, the effect of PAR is highly non-linear in the 
range of PAR values generally encountered. As PAR increases, primary 
production becomes quickly insensitive to PAR. Saturation occurs at PAR values 
as low as 200 Wm'2 when phytoplankton concentration is as high as Imgm'^. 
Unlike over land, where primary production becomes independent of leaf area 
index at high values of the index, primary production over the ocean increases 
substantially even when phytoplankton concentration is high. 


Rg. , provides some insight about the accuracy 

llLriries in the relationship between primary production and PAK, the 

LTth.TpTimary'’ production can be considered directly 
r o n ino Wm-2 over land and 0-50 Wm-2 over the ocean), the 10% accu y 
primary production translates equally into a 10% accuracy on PAR, but 20% an 
35 % accuracies will be suHicient over land and ocean respectively, when PAR « 
above 300Wm-^. Thus, a better relative accuracy on PAR is required at low PAR 
values, which occur either at low solar zenith angles or in the presence , 

under those conditions, unfortunately, satellite algorithms are less accurate. In 
view of available primary production models, however, the accuracy 
requirements on PAR may be relaxed. The models generally incorporate the fac 
thm the growth rate of many plants is proportional to the rate 
energy absorption by chlorophyll pigments, but this rate (absorbed an 

efficSicy fairs (functions of plant type, environmental condidons) are difhcuU 
to estimate with accuracies comparable to those mentioned above for PAR. In 
other words, useful estimates of primary production may still be obtaine wi 

larger errors on PAR. 

If we are to understand truly the interactions between the biosphere and the 
atmosphere and their effects on climate, we need to know the S“S"P 
distribution and temporal variability of primary production and, thus, PAR ov 
the globe. Until recently, our information was based on surface pyranome er 
networks (essentially over land) and a few PAR sensors deployed during researc 
experiments. The networks are clearly insufficient for global change studies, the 
oceans and polar regions, in particular, are virtually not sampled, and long-term 
time series (from well-maintained, regularly-calibrated sensors) are only existent 
at a few locations. Furthermore, pyranometers measure insolation, or t e so ar 
radiation incident in the spectral range 0.4-4 pm, and the relationship between 
PAR and insolation depends on atmospheric conditions and radiation geome ry 
(e.g., Baker and Frouin, 1987; Pinker and Laszlo, 1992). Clouds, which do no 
absorb at PAR wavelengths but do absorb substantially m the near-in ' 

increase the ratio of PAR and insolation. Data from the pyranometer networks 
can be complemented by estimates based on empirical formulas and cloud 
observations made routinely at meteorological stations (e.g., cloud cover. 


type). The formulas, unfortunately, have been established locally and are 
therefore difficult to apply confidently over large areas. Moreover, the dataset is 
uneven and too often of questionable quality. Because of these limitations, we do 
not yet have a clear picture of PAR'S modes of variability over the globe. 
However the situation is being changed with existing earth-observing satellites, 
which provide regular coverage of the earth and observations of the basic cloud 
properties governing PAR variability. 


Satellite Algorithms 

While numerous studies have been devoted to estimating insolation from 
satellite data (e.g., Tarpley, 1979; Gautier et al., 1980; Moser and Rashke, 1984; 
Pinker and Ewing, 1985; Dedieu et al., 1987; Darnell et al., 1988), only a few 
satellite-based methods have been proposed for PAR, including the methods of 
Frouin and Gautier (1990), Eck and Dye (1991), and Pinker and Laszlo (1992). Part 
of the reason is that for many applications involving small space and time scales 
PAR can be measured directly. Furthermore, it has often proven satisfactory to 
take PAR as a more or less constant fraction of insolation. Deducing PAR from 
insolation, in fact, is the basis of Pinker and Laszlo's (1992) method, which can be 
qualified as indirect (requires an insolation estimate). Noting that meteorological 
satellites (except METEOSAT) carry instruments that measure in spectral 
channels resembling more the PAR wavelength range than the entire solar 
spectrum, Frouin and Gautier (1990) use the satellite radiances directly. 
Uncertainties in insolation are not propagated in that case, and the modeling of 
cloud effects is simplified (no narrow-band to broad-band transformation is 
necessary, and cloud absorption vanishes in the equations). This method, also 
used by Eck and Dye (1991), can be qualified as direct (does not require an 
insolation estimate). In what follows, we contrast the salient features of the 
indirect and direct methods, and we present typical results. 

a Indirect approach 

In Pinker and Laszlo's (1992) method, insolation (estimated using the model of 
Pinker and Ewing, 1985) is converted into PAR using a relationship established 
theoretically. This relationship depends on atmospheric conditions, which need 
to be specified. Under clear skies, the ratio of PAR to insolation varies little 


,„„nd 0 48 except at high solar zenith angles or extreme (low as well as htgh) 
“r vapor'amounts (Fig. 2), and the effect of aerosol turbidity is only signiftc^l 
when horizontal visibility is less than 10km. This suggests that the raho of PAR 
to insolation can be considered constant to a good degree of approximation unde 
deaTskies. The situation is quite different under cloudy skies. Cloud optical 
thickness substantially changes the ratio of PAR and insolation, which can vary 
by more than 50% at low solar zenith angles (Fig. 3). This variability m the PAR^ 
to-insolation ratio is corroborated by in-situ measurements (Fig. 4). Pinker 
Laszlo's (1992) procedure is to therefore apply a variab e conversion 
LLlat on estiltes. This factor depends on cloud optical thickness and 
fractional amount, parameters derived from the satellite measurement. 
Applying this method to hourly ISCCP Cl data at 250 km resolution. Pinker an 
Laszlo (1992) have produced the first global map of monthly PA , e ^ 

demonstrating that global satellite datasets produced within the tome of ISC 
will soon result in a global, long-term climatology of PAR. Owing to no - 
linearity, conversion factors are applied before averaging 
insolation estimates over daily and longer time scales. It may be possible to app y 
conversion factors to daily or monthly insolation estimates without signif.can 
loss of accuracy. Fig. 5, established from surface data collected during the First 
ISLSCP Field Experiment (FIFE), shows that the PAR fraction o ai y mso 
remains fairly constant regardless of cloud conditions. The same filing was 
reported by Howell et al. (1983) and Rao (1984) on a monthly time scale. At those 
time scales the PAR fraction variability due to clouds is reduced because it 
strongly depends on sun zenith angle (Fig. 4). It may, therefore, prove use u o 
deduce PAR directly from the various insolation climatologies (satellite or other) 
currently available or being produced al daily or longer time scales (ej;. B.stop 
and Rossow, 1991). Fig. 6 shows a typical example obtained with METEOS 

data. 


b Direct approach 

Frouin and Gautier's (1990) method is based on the formalism developed by 
Gautier et al. (1980) for insolation, that only requires slight modifications (m ac , 
simplifications) to be applicable to PAR. Cloud absorption vanishes in j:loudy 
sky model equations (clouds do not absorb at PAR wavelengths), and t e c ea 
sky model coefficients represent the PAR spectral interval instead of the total 


solar spectrum. Cloud albedo, the governing cloud parameter, is computed as in 
Gautier et al. (1980) from geostationary satellite observations in the visible and 
near-infrared. Since the solar channels of geostationary satellite instruments 
(except the METEOSAT radiometer) mostly capture radiation in the visible no 
narrow-band to broad-band conversion of cloud albedo is necessary. Because of 
these simplifications in the radiative transfer modeling, we expect, at least in 
principle, more accurate results for PAR than for insolation. Furthermore, by 
estimating PAR directly from the satellite radiances, uncertainties due to errors 
on insolation estimates and on the ratio of PAR and insolation, which are 
inherent to Pinker and Laszlo's (1992) method, are avoided. Fig. 7 shows, for 
selected days during the First ISLSCP Field Experiment (FIFE), the diurnal 
variation of measured PAR at the study site (Konza prairie, Kansas) and the 
corresponding satellite estimates at particular times during the day. In the figure, 
the in-situ values are half-hourly averages and the satellite estimates are spatial 
averages over the FIFE area (15xl5km). Julian days 222, 223, 226, and 227 are 
mostly clear, whereas days 224 and 225 are cloudy. In general, the satellite 
estimates compare well with the measured values; they describe the diurnal cycle 
properly. The larger discrepancy observed during days 225 and 226 may be linked 
to spatial cloud variability, which is not accounted for in the modeling (see in the 
next section the discussion about effects of cloud heterogeneity). For daily 
averages. Fig. 8 shows the temporal variation of PAR at the site during intensive 
field campaigns 2 and 3. Satellite estimates correspond to measurements to 
within 10-15 Wm-2 (about 10%), and more than 85% of the observed variance is 
explained. These comparisons, although performed for a single geographical 
location, are strongly indicative of the method's ability to quantify PAR 
variability on daily or longer time scales. 

Instead of using radiances in the visible and near-infrared, Eck and Dye (1991) use 
radiances (or, equivalently, reflectances) in the ultraviolet and test their method 
with Total Ozone Mapping Spectrometer (TOMS) data. Noting that cloud 
reflectivity is constant across ultraviolet and PAR wavelengths and that clouds 
do not absorb radiation at ultraviolet and PAR wavelengths, they parameterize 
the effect of clouds on PAR as a simple, linear function of TOMS ultraviolet 
reflectance. Cloud-screening is achieved by applying a threshold technique, and 
the authors argue that using data in the ultraviolet makes it easier to 
discriminate clouds from high-albedo background surfaces, except for ice and 


snow! The cloud-screening, however, may not be efficient because the TOMS 
data are in the form of monthly, 500x500km averages, and there is no way of 
assessing from the TOMS data alone whether the 500x500 km areas are partially 
contaminated by clouds or not. Furthermore, the radiative transfer modeling is 
rather crude (e.g., no correction is performed for molecular scattering above the 
clouds). Nevertheless, the effects do not appear significant on a monthly time 
scale (individual errors somewhat cancel out), as comparisons with surface 
measurements, which reveal less than 6% relative differences, demonstrate (Fig. 
9). 


Issues 


The satellite algorithms so far proposed to monitor the variability of PAR over 
the globe utilize data from instruments (e.g.. Advanced Very High Resolution 
Radiometer, Visible and Infrared Spin-Scan Radiometer) that are generally not 
calibrated after launch. These instruments have been shown to exhibit 
significant, even large changes in sensitivity. The resulting errors on PAR can be 
important, as Fig. 10 illustrates. For a cloud containing lOOgm'^ of liquid water at 
40°N, for instance, a 10% loss of sensitivity translates into errors of up to 50Wm'2 
on monthly averages. Degradation of that amplitude is quite common, as many 
studies have demonstrated (e.g., Frouin and Gautier, 1987; Staylor, 1990; 
Whitlock et al., 1990; Brest and Rossow, 1992). Therefore, unless a check-of- 
calibration is maintained on a regular schedule during the lifetime of the 
satellites, and instruments from various satellites cross-calibrated properly, it 
will be difficult to extract a meaningful signal for climate studies from observed 
changes. 

Another issue deals with cloud spatial heterogeneity. The satellite estimates are 
generally less accurate in conditions of partial (broken) cloudiness (see for 
instance the results for days 224 and 225 in Fig. 6). This is not surprising as clouds 
are considered plane-parallel in the modeling, and top-of-atmosphere radiance is 
often assumed to be isotropic. Drastic assumptions of that sort are necessary, 
however, to close the system of equations and reduce the problem to one of 
estimating PAR from a single top-of-atmoshere radiance measurement. The 
drawback is that large errors on the PAR estimates may be introduced for some 


situations. Broken clouds, in particular, can significantly affect the spatial 
distribution of PAR, as the Monte Carlo simulations of Fig. 11 illustrate. For the 
cloud field considered, namely a regular network of cylinder clouds characterized 
by a radius of 0.5 km, a geometrical thickness of 0.2 km, an optical thickness of 12, 
and a distance between clouds of 2.5 km (typical conditions observed during the 
HFE experiment on August 9, 1989), the cloud transmittance (flux transmittance) 
exhibits strong spatial variance, depending on whether the sun disk is obscured 
by the clouds or not, and reaches over 110% in areas directly illuminated by the 
sun. In other words, more sunlight that would be observed in clear sky 
conditions reaches the surface in those areas. This effect, observed by many 
investigators on pyranometer traces, cannot be reproduced by assuming plane- 
parallel clouds. Furthermore, depending on the cloud field, it may not cancel out 
on daily or monthly averages. 

To assess the accuracy of the PAR estimates, one needs to compare them to other 
data, particularly surface measurements. The networks of well-calibrated PAR 
sensors, unfortunately, are generally inadequate, even over the continents. In 
fact, the networks of surface radiation instruments have been designed to 
monitor insolation not PAR -and deducing PAR from insolation is subject to 
uncertainty (see above). Furthermore, the satellite estimates are instantaneous 
whereas the surface measurements are local, making it difficult to compare the 
two types of data. When the method utilizes coarse resolution pixels (see Eck and 
Dye, 1991), validation by surface measurements becomes very difficult. One 
alternative is to compare low resolution PAR estimates to estimates obtained 
from higher resolution data using a validated satellite method; but the procedure 
is far from optimum. It is clear, however, that without proper validation 
strategy, satellite PAR estimates will not find quantitative use in global change 
studies of the carbon cycle. 


Summary and Recommendations 

Developing methods for estimating PAR from satellites is a recent activity that 
has strongly benefited from the work performed on insolation by many 
investigators. Satellite estimates of insolation can be converted accurately into 
PAR, which makes it possible to exploit already existing datasets (satellite and 


other). From the radiative transfer point of view, the problem of deriving PAR 
from top-of-atmosphere radiances in the visible is simpler for PAR than 
insolation because narrow-band to broad-band transformation is not necessary, 
and cloud absorption does not need to be parameterized (clouds do not absorb in 
the visible). In situations of partial cloudiness for which plane-parallel theory 
does not apply, the problem is as complex as for insolation. Although limited 
comparisons have been made, an inaccuracy smaller than 10% on a monthly 
time scale appears feasible by the methods reviewed. In view of the existing 
models of primary productivity, which involve terms other than PAR more 
difficult to estimate, a 10% inaccuracy is more than sufficient and should allow a 
correct description of the month-to-month PAR variability and reveal large scale 
seasonal and interannual phenomena. 

Many of the recommendations of previous workshops on surface radiation 
budget (e.g., Suttles and Ohring, 1986; Sellers et al., 1990) are in order for PAR. 
Some effort particularly should be put to rigourously specifying the required 
accuracy on PAR. As suggested by Sellers et al. (1990), sensitivity studies are 
necessary, but it is unrealistic to expect that they will provide a complete, 
universal answer; the space and time scales of geophysical phenomena 
influenced by PAR are too varied. Whatever the phenomenon under study it 
will always be safe to define the required accuracy so that the variability of PAR 
over the phenomenon's characteristic space and time scales, generally 
observable, is described properly. 

Regarding the calibration issue, a lot of progress has been made during the last 2- 
3 years to monitor sensor degradation of meteorological satellites, those used for 
PAR, after launch (e.g., within the frame of ISCCP, NOAA and GOES pathfinder 
activities). Despite the numerous studies a consensus sometimes has been 
difficult to reach on the calibration coefficients to use for some sensors. This 
underscores the need for instruments that possess on-board calibration 
capabilities and for detailed, realistic calibration plans prior to launch. In view of 
the potential of radiometers carried by meteorological satellites for PAR 
monitoring, it appears in order to equip future versions with a proper calibrator 
for their solar channels. In the long run, the strategy might prove more 
economical and rewarding, since costly aircraft calibrations would be downsized, 


and scientists would be relieved from tedious, time-consuming calibration tasks 
they too often have to perform themselves at the expense of other work. 

Regarding validation activities, care should be exercised when satellite-derived 
estimates are compared with in-situ measurements. In general, the two 
quantities are not the same. On the one hand, satellite-derived values are 
instantaneous and averaged spatially; on the other hand, surface measurements 
are local and averaged temporally. The space and time scales at which the 
comparisons should be made need therefore to be selected rationally, and 
instrument networks designed accordingly. Using a single instrument is not 
optimum; dense networks are more appropriate. Such networks were installed 
during various ISLSCP experiments but covered a limited time period. They 
should be operated continuously at sites representing world-wide conditions and 
include measurements of other parameters (e.g., cloud properties) to test 
individual parameterizations in the models. PAR sensors, which are 
inexpensive, should also be deployed to complement the networks of 
pyranometers already in place, at least in representative areas of the globe. Effort 
should also be made to create a database of PAR measurements from various 
research experiments and make it available for validation studies. Comparisons 
of algorithms such as those for insolation should be made (e.g., Whitlock et al., 
1990 ), but with the purpose of understanding the advantages and drawbacks of 
each algorithm instead of selecting one. 

One of the major limitations of the methods is their inability to provide 
reasonable estimates when plane-parallel theory is not applicable (case of broken 
clouds, liquid water spatial heterogeneity). Efforts to improve the techniques 
should therefore focus on situations of cloud heterogeneity. One approach is to 
perform radiative transfer calculations for realistic cloud fields, determine the 
cloud parameters that govern departures to plane-parallel theory, and investigate 
relationships between the governing cloud parameters and observable cloud 
characteristics (texture, moments, etc.). If this approach proves suitable, current 
strategies to create long-term, large-scale satellite datasets might have to be 
reviewed to include those cloud characteristics. 

Two other aspects of the methods should also be addressed, namely the presence 
of snow or ice at the surface and diurnal sampling. Over snow and ice it is not 



easy to distinguish clouds, and the methods proposed would likely fail. Efforts 
should be made to improve the methods in those situations, all the more as the 
polar oceans cannot be neglected in studies of the global carbon cycle because of 
their high primary productivity. Regarding diurnal sampling, the success of the 
satellite methods generally resides in their ability to sample diurnal cloud 
variability properly. Polar-orbiting satellites do not provide adequate sampling at 
middle and low latitudes. Statistically obtained correction factors may be used, 
but they do not offer the solution. The problem may be obviated, however, by 
complementing data from polar-orbiting and geostationary satellites, as is 
currently being done to generate ISCCP datasets. 

The sensors adapted to PAR monitoring from space are not limited to those used 
in the algorithms so far proposed. Other instruments, scanners as well as wide- 
field-of-view radiometers, have not yet been exploited, in particular those of the 
Earth Radiation Budget Experiment. In fact the current algorithms can be easily 
modified to become applicable to those sensors. Furthermore, their longevity, 
careful calibration and characterization, as well as the continuity of the mission 
well beyond the end of the century (Clouds and Earth's Radiant Energy System, 
CERES, investigation), make them an ideal tool for studying PAR'S inter-annual 
modes of variability and related questions of climate change. Looking ahead, 
apart from the future versions of meteorological satellites and the CERES 
scanner a battery of instruments will be available for PAR monitoring during the 
Eos era, in particular the MODerate resolution Imaging Spectrometer (MODIS) 
and the MEdium Resolution Imaging Spectrometer (MERIS). Our prospects are 
good for an accurate, long-term climatology of PAR over the globe. 
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Figure Captions 


Fie. 1. Primary production as a function of photosynthetically available radiation, 
PAR (a) Case of a green canopy with horizontal leaves and a leaf area index 
ranging from 0.1 to 6 (after Sellers, 1985). (b) Case of a homogeneous, 2QoC ocean 
containing 0.3, 1, 3, and 10 mgm-3 of chlorophyll pigments. 

Fig. 2. Ratio of photosynthetically available radiation, PAR, and insolation as a 
function of water vapor amount (top), ozone amount (middle), and aerosol type 
and visibility (bottom). (After Baker and Frouin, 1987.) 

Fig. 3. Ratio of photosynthetically available radiation, PAR, and insolation as a 
function of cloud optical thickness and sun zenith angle. (After Pinker 

Laszlo, 1992.) 

Fig 4 Surface-measured ratio of half-hourly photosynthetically available 
radiation, PAR, and insolation as a function of satellite-derived instantaneous 
cloud liquid water content during the First ISLSCP Field Experiment. The ratio 
varies between 0,25 and 0.75, corroborating theoretical calculations. 

Fig. 5. Surface-measured ratio of daily photosynthetically available radiation, 
PAR, and insolation as a function of satellite-derived daily cloud cover during 
the First ISLSCP Field Experiment. At this time scale the PAR fraction variability 
is small, with values ranging between 0.43 and 0.52. 

Fig. 6. Monthly photosynthetically available radiation, PAR, derived from 
METEOSAT data for June 1990. Monthly insolation was first obtained using the 
method of Dedieu et al. (1987) and PAR was then deduced by taking the ratio of 
PAR and insolation equal to 0.45. 

Fig 7 Surface-measured and satellite-derived photosynthetically available 
radiation, PAR, for selected days during the First ISLSCP Field Experiment. 
Satellite estimates are instantaneous whereas measured values are half hourly 

averaged. 


Fig. 8. Surface-measured and satellite-derived daily photosynthetically available 
radiation, PAR, during the second and fourth intensive field campaigns of the 
First ISLSCP Field Experiment. 

Fig. 9. Satellite estimates of monthly photosynthetically available radiation, PAR, 
versus surface estimates from pyranometer measurements adjusted to PAR. 
(After Eck and Dye, 1991.) 

Fig. 10. Typical error on satellite-derived monthly photosynthetically available 
radiation, PAR, due to a 10% increase in the calibration gain, g, of the sensor's 
solar channel. Clouds contain 100gm*2of liquid water, and the clear atmosphere 
contains 0.3 atm-cm of ozone and aerosols of continental type and optical 
thickness of 0.22 at 550 nm. Latitude is 390N. As fractional cloud coverage, N, 
increases, the error increases in magnitude, reaching -50 Wm'2 in June and July. 

Fig. 11. Monte Carlo simulations of the spatial distribution of cloud 
transmittance (in percent) on August 8, 1989 at 13:30 local time over the Konza 
prairie, Kansas. The clouds are cylindrical of radius 500m, separated by 2,500m, 
and located between 2,000 and 2,200m (geometrical thickness of 200m). The cloud 
optical thickness is 12. When the sun disk is not obscured by clouds, cloud 
transmittance reaches 113%, indicating that the surface receives more 
photosynthetically available radiation than in clear sky conditions. 
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Fig. 1. Primary production as a function of photosynthetically available radiation, 
PAR. (a) Case of a green canopy with horizontal leaves and a leaf area index 
ranging from 0.1 to 6 (after Sellers, 1985). (b) Case of a homogeneous, 20°C ocean 
containing 0.3, 1, 3, and 10 mgm'^ of chlorophyll pigments. 





Ozone Amount (atm-cm) 


0.50 

1 

C 1 


maritime 


CONTINENTAL \ 

1 


1 

o 

1 


1 

S 0.48 


o 

1 

c/5 

C 

_ - 

^ 0.47 

^ ^ ” 



< 

-*• 

^ 0.46 

/ US 62 


0 =30*’ 


/ • 

0.45r . j 


0 10 20 50 40 50 

Visibility (km) 


Fig. 2. Ratio of photosynthetically available radiation, PAR, and insolation as a 
function of vv^ater vapor amount (top), ozone amount (middle), and aerosol type 
and visibility (bottom). (After Baker and Frouin, 1987.) 
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Fig. 3. Ratio of photosynthetically available radiation, PAR, and insolation as a 
function of cloud optical thickness and sun zenith angle. (After Pinker and 
Laszlo, 1992.) 
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Fig. 5- Surface-measured ratio of daily photosynthetically available radiatioi\, 
PAR, and insolation as a function of satellite-derived daily cloud cover during 
the First ISLSCP Field Experiment. At this time scale the PAR fraction variability 
is small, with values ranging between 0.43 and 0.52. 
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Fig. 8. Surface-measured and satellite-derived daily photosynthetically available 
radiation, PAR, during the second and fourth intensive field campaigns of the 
First ISLSCP Field Experiment. 
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Fig. 9. Satellite estimates of monthly photosynthetically available radiation, PAR 
versus surface estimates from pyranometer measurements adjusted to PAR 

(After Eck and Dye, 1991.) 
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Fig. 11. Monte Carlo simulations of the spatial distribution of cloud 
transmittance (in percent) on August 8, 1989 at 13:30 local time over the Konza 
prairie, Kansas. The clouds are cylindrical of radius 500m, separated by 2,500m, 
and located between 2,000 and 2,200m (geometrical thickness of 200m). The cloud 
optical thickness is 12. When the sun disk is not obscured by clouds, cloud 
transmittance reaches 113%, indicating that the surface receives more 
photosynthetically available radiation than in clear sky conditions. 
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1. INTRODUCTION 

During recent years public and scientific interest 
in the planet earth's climate system has grown 
remarkably. This is attributable, in part, not only to 
an increased concern about possible environmental 
problems (e.g., global warming due to the famous 
greenhouse effect), but also to the emergence of 
technologies enabling new forms of climate study. 
Meteorological satellites, for instance, now provide 
global views of the earth's atmosphere and surface, 
while increasingly powerful computers permit 
more realistic simulations of atmospheric and 
oceanic circulations. 

One of today's main climate research topics 
concerns ocean-land-atmosphere interactions on a 
global scale. This includes studying energy, water, 
and momentum exchanges at the media s 
interfaces. Characterizing the energy budget with a 
reasonable spatial and temporal accuracy is of major 
importance to various applications, in particular: 1) 
diagnostic studies of heat transport, cloud forcing, 
radiative heating; 2) specification of boundary 
conditions for global circulation models (GCM); 3) 
amelioration of sub-grid process parameterizations 
in such models; 4) validation of climate models, 
and 5) determination of long-term trends (NASA, 
1986). 

Moreover, since the land surface, unlike the 
oceans, has a small hea-t capacity, radiation and 
water budgets are strongly related: the net heat 
budget, which is the balance between radiation and 
the sensible and latent heat fluxes, is close to zero 
when averaged over time periods longer than a 
day. The latent heat flux is directly related to the 
evaporation rate and is consequently determined 
largely by the sum of the radiation and sensible heat 
budgets. 

The net longwave irradiance flux at the surface, 
owing to its small variability, is less important than 
the shortwave, or solar, one. It cannot be neglected, 
however, particularly in the presence of clouds 
which decrease surface insolation while increasing 


the downward component of the longwave flux. 
Several methods have been proposed to estimate 
the downward longwave flux at the surface from 
satellite measurements (see, for instance, 
al. 1983; Frouin et al., 1988; Gupta, 1989). Although 
some of them have been used to generate global 
fields of net longwave flux, few validation studies 
have been made (WMO-ICSU, 1984). We therefore 
investigated the performances of two of them, 
choosing among the most promising ones, by 
comparing their computations with carefully taken 


In the first part of this paper we describe the data 
used in this study. The data were collected during 
the First ISLSCP Field Experiment (FIFE). In 
analyzing the data quality, we emphasize the 
difficulty to acquire reliable longwave irradiance 
measurements and, therefore, validate radiative 
transfer models. We then present briefly the two 
methods selected to estimate downward longwave 
flux at the surface, namely those of Frouin et al. 
(1988) and Gupta (1989). Finally, the two methods 
results are compared and verified against in situ 
pyrgeometer measurements. 


2. DATA 

FIFE an international surface-atmosphere 
experiment, took place in 1987 at and around the 
Konza Prairie Long Term Ecological Research 
(LTER) site near Manhattan, Kansas (Sellers et al., 
1988). The FIFE objectives were to gather the 
necessary data to permit interpretation of satellite 
observations suitable to infer climatologically 
significsnt land surface parameters. 

The experimental area is a 15x15 km square 
with various topological features including burned 
and unburned plateau, slope, and creek 
Although the central portion of the area is studied 
continuously and additional parameters are 
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measured throughout the year, the experiment 
concentrated on four Intensive Field Campaigns 
(IFCs) corresponding to four different slates of the 
surface vegetation: 


First IFC : 
Second IFC : 
Third IFC : 
Fourth IFC : 


"green up" 

"peak greenness" 
"dry down" 
"senescence" 


May 26th to June 6th 
June 25lh to July 15th 
August 10th to August 21st 
October 5lh to October 15th 


Our study uses 4 different types of data acquired 
during the IFCs, namely surface measurements, 
atmospheric soundings, sky photographs, and 
satellite-derived parameters. 

Several surface stations collected radiative flux 
measurements. At the time of our study, only data 
from stations 2 and 38 (see Fig. 1) were available 
through the FIFE information system. These 
stations were operated by Eric Smith's group 
(University of Florida) during the four IFCs. A first 
check of the data was made by comparing hourly 
measurements from one station to those from the 
other (Fig. 2). Some of the station 2 measurements 
are obviously bad (values higher than 450 Wm'^). 
A plot of the time series of these measurements 
shows that the unrealistic values are limited in 
time to the first days of the first IFC. At other times, 
the two stations’ data are well correlated with each 
other, but exhibit a relative bias (higher 
measurements at station 2) which is rather small 
for low flux values (a few Wm‘^) but up to 20 Wm'^ 
for the highest fluxes. Since we know with 
confidence that station 2 provided inconsistent 
measurements at the beginning of the experiment, 
we have an objective reason to favor station 38. It 
is possible, however, that the problems experienced 
by station 2 did not propagate a bias along the four 



Fig. 1: General presentation of FIFE experimental field. 

Station 2 and 38 are the two longwave measurement sites used 
in this study. The black circle indicates the soundings' launch 
location. 
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Fig. 2; Comparison of coincident hourly downwelling longwave 
fluxes measured at station 2 and 38 during the four IFCs of FIFE. 

The fluxes are inWm . 

IFCs. Further discussion with E. Smith convinced 
us that more work has to be done to favor one 
station rather than the other. Before the 
experiment, the two pyrgeometers had been 
carefully intercalibrated and gave very close values. 
The reason for the discrepancy found during the 
IFCs is still unclear, but a different instrument 
orientation is suspected (E. Smith, 1989, personal 
communication). 

During the four IFCs there was a strong diurnal 
cycle of downwelling longwave flux at the surface. 
The cycle showed a maximum at about 5:00 PM 
local time and a minimum around sunrise. This 
does not exactly follow the temperature cycle as 
could have been expected. However, the 
development of convective clouds in the 
afternoon, leading to an increase of the flux, may be 
responsible for this delayed maximum. The 
amplitude of the cycle is about 30 Wm for IFCs 1 
and 3, close to 20 Wm'^ for IFC 2 and only 10 Wm'^ 
for IFC 4. The 4 IFCs are not long enough for those 
figures to be statistically significant; the smaller 
amplitude for IFC 4 may be explained by its later 
period in the year, leading to a smaller daily cycle 
temperature. 

Cloudiness was another parameter acquired 
from the surface. This was done by an uplooking 
automatic camera. Cloudiness was deduced from 
the percentage of cloudy areas in the photographs, 
taking into account the viewing angle of each area 
relative to the vertical. Since this camera was 
operating in the visible spectrum, no cloudiness 
information could be obtained at night. In the 
study that follows, we considered that the sky was 
clear when, according to the photographs, less than 
1 /8 of the sky was covered by clouds. 

The atmospheric measurements used in our 
study were obtained from radiosondes launched 
from the experimental site (see Fig. 1). These 




measured pressure, temperature, and wet bulb 
temperature, from which the water vapor mixing 
ratio was derived. We first had to clean the data: 
the pressure, which usually decreased smoothly, 
sometimes showed 50 mb jumps. This problem 
was resolved by removing corresponding 
measurements. More bothersome, we found that 
the water vapor mixing ratio occasionally exhibited 
negative values! This was the result of improper 
wet bulb temperature measurements when the air 
temperature was low; in such cases, the water 
surrounding the thermometer simply freezes and 
the wet bulb measurement loses its significance. 
We overcame this problem by removing all water 
vapor information for which the wet bulb 
temperature was negative. The water vapor profile 
was then extrapolated to higher levels using 
climatological profiles, taking into account the 
valid information about the lower layers and the 
actual temperature profile so that the mixing ratio 
could not larger than the saturation level. 

The satellite products we used are the pre- 
processed TIROS-N Operational Vertical Sounder 
(TOVS) sounding data. The TOVS instrument, a 
multi-frequency radiometer designed to remotely 
sense, among other variables, atmospheric 
temperature and water vapor mixing ratio profiles, 
has been flying for many years on board the NOAA 
series satellites (Kid well, 1981). The TOVS observa- 
tions were processed with NOAA/NESDIS's 
operational algorithm (Kidwell, 1981). The TOVS 
products include surface pressure, surface 
temperature, air temperature at 15 pressure levels, 
three layers of water-vapor content, cloud top 
pressure, and effective cloud amount. 


3. METHOD 

Two methods that determine downwelling 
longwave flux at the surface from a description of 
the atmosphere were selected for this study. These 
are the method of Frouin et al. (1988), which 
employs the fairly detailed radiative transfer model 
designed by Morcrette (Morcrette, 1984; Morcrette 
and Fouquart, 1985; Morcrette et al., 1986), and a 
simpler parameterization recently proposed by 
Gupta (1989). 

Morcrette's model uses as input temperature, 
water vapor and ozone mixing ratio, and cloud 
amount profiles at levels which can be chosen by 
the user. The model divides the infrared spectrum 
into six intervals and explicitly fakes into account 
absorption by the radiatively important atmos- 
pheric gases. 

To determine the model input parameters, 
Frouin et al. (1988) proposed various methods. In 
all the methods, temperature and water vapor 
mixing ratio are obtained from TOVS data. In the 
most refined method (method 1), cloud base 
altitude is deduced from cloud top altitude and 
liquid water path, assuming a liquid water 


distribution within the cloud. In the other 
methods, simplifying assumptions are introduced 
which include directly relating liquid water path to 
cloud geometrical thickness (method 2), fixing the 
cloud geometrical thickness to its climatological 
value (method 3), and, finally, parameterizing the 
cloud effects only as a function of fractional cloud 
coverage (method 4). In the calculation presented 
later in this paper, we only employed method 3 
using a thickness of 500 m. 

The Gupta (1989) scheme, a much simpler one, 
limits the computational burden. It is a simple 
parameterization developed using regression 
techniques applied to a set of atmospheric profiles 
and the corresponding longwave fluxes at the 
surface, as determined using a radiative transfer 
model similar to Morcrette's. Although the Gupta 
(1989) method can be applied, in principle, to any 
complete atmospheric description, it has been 
explicitly designed to use TOVS products; the input 
parameters are surface temperature, 1000-850 mb 
and 850-700 mb temperatures, cloud base tempera- 
ture, total water vapor content, water vapor content 
below the cloud, and cloud amount. Since cloud 
base information is not accessible from the TOVS 
instrument, the method uses the estimated cloud 
top pressure and assumes a mean cloud thickness 
of 50 mb; the cloud base temperature and water 
vapor below the cloud are then obtained from the 
corresponding profiles. The same assumptions 
have been used in the study that follows. 

One of the goals of this study is to accurately 
compare longwave flux estimates from the models 
with in situ measurements. The objective is to 
estimate how precise the models really are because 
although model intercomparisons have been 
performed, careful compariso'ns with in situ 
observations have not. When model estimates and 
measurements show large differences, the 
differences can be attributed not only to the 
uncertainty of the model itself, but also to the 
measurement or to the quality of the input 
atmospheric profile. To validate a model, we 
therefore need high quality and reliable longwave 
radiation measurements and atmospheric profiles. 
Temperature and water vapor profiles can be 
obtained with a reasonable accuracy from 
radiosonde data. This is not the case for the cloud 
parameters, particularly cloud base pressure and 
effective cloud emissivity. Therefore, in cloudy 
conditions, a discrepancy between estimated and 
measured fluxes will not be conclusive, and we 
decided to concentrate on clear periods. Clear sky 
radiosoundings were selected using the 
information given by the sky camera. 

Another goal of the study was to evaluate the 
Frouin et al. (1988) and Gupta (1989) methods in 
estimating the longwave flux from satellite 
observations over the FIFE area. 

This includes not only the algorithms 
themselves, but also the inversion process that 
retrieves geophysical parameters from TOVS data. 


In the evaluation, we included cloudy observations 
as well as clear ones. Only TOVS observations 
closer than 200 km from the center of the 
experimental area were used. In order to keep the 
more accurate TOVS products .'ve selected only 
"dass-1" soundings as defined in Guptas 
paper In other words, the radiative calculations 
were performed when all the parameters 
(temperature, water vapor, and cloud parameters) 
were present in the products. 


4. RESULTS AND DISCUSSION 

Figure 3 illustrates a comparison between the 
Morcrette et al. (1986) and Gupta (1989) model 
outputs when using clear sky radiosoundings 
acquired during FIFE. The longwave estimates 
range from 235 to 420 Wm’l • Morcrette's and 
Gupta's values are very close; the standard 
deviation of the differences is only 3.5 Wm' and 
the correlation coefficient is 0.999. However, a large 
bias of (15 Wm’^) between the two estimates is 
observed. This bias is larger than the accuracy 
required for estimating longwave flux at surface for 
most scientific applications; it is therefore 
important to evaluate which model is in fact most 
accurate. This can be accomplished by comparing a 
time series of the longwave flux measurements to 
the values found when applying the two models to 
the sounding data in clear conditions. Figure 4, 
which shows an example of such a cornparison, was 
obtained for the third IFC at station 38 using 
Morcrette's model. Similar comparisons for other 
IFCs, stations, or models are not presented here but 
give comparable results. These results show that, 
in general, the models accurately reproduce the 
variability of the clear sky flux. In particular, the 
increase in downwelling longwave flux associated 
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Fie. 3: Comparison between the Morcrette et al. (1986) and 
Gupta (1989) model clear flux estimates. The input profiles are 
the FIFE radiosoundings launched during clear periods of the 

four IFCs. Values arc displayed in Wm . 
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Fie. 4: Time scries of in situ-measured longwave fluxes and 
those estimated by the Morcrette el al. (1986) model using FIFE 
radiosoundings as an input. The estimates were done only 
during clear periods, as inferred by the sky-looking camera, 

Values arc in Wm"^. This figure is limited to the third IFC 
(August 10 to August 21). 

with the diurnal variations of air temperature is 
well-reproduced by the model. It is also interesting 
to note that the maximum values observed in late 
afternoon in clear conditions are almost as large as 
the overall maximum values observed. This 
suggests that, over land, the cloud effects are no 
more important in increasing the downwelling 
longwave flux than the atmospheric temperature 
(and moisture) effects are. As already mentioned, 
there is large uncertainty in the measurements 
themselves. The bias observed between the models 
is comparable to the one between the measure- 
ments Gupta's model reproduces well the station 2 
measurements, while Morcrette's models does the 
same for station 38. As a consequence, the results 
do not allow one to conclude which model is more 
precise. 

Figure 5 depicts the time series of downward 
longwave flux measurements at station 38 
and corresponding satellite estimates (open circles) 
for clear and cloudy conditions. The satellite 
estimates were obtained with the Frouin el ^ ^ 
(1988) method 3 (Fig. 5a) and the Gupta (1988) 
method (Fig. 5b). Even though Fig. 5 was obtained 
for the third IFC, the results are similar for the 
three other IFCs. Occasionally several satellite 
estimates are displayed for the same time, they 
correspond to different satellite soundings during 
the same satellite pass within the 200 km distance 
limit. This discrepancy, which can be larger than 
20 Wm'^, is mainly due to variable cloudiness 
parameters. These figures indicate that when using 
satellite cloud and atmospheric products instead ot 
in situ observations, the estimations of the 
downwelling longwave radiation degrade 
significantly. In comparison with the model results 
presented earlier, it suggests that the TOVS satellite 
products used are not sufficient for achieving an 
accuracy of 10 Wm'^. Both methods provide results 
that are smaller, in general, than the in situ 
measurements. The reason for this is not 
understood. 




Fig. 5: Time series of in silii measured loTigwave flux and thal 
cslimatcd by the Frouin cl al. (1988) method 3 (Fig. 5a) and the 
Gupta (1989) model (Fig. 5b) using TOVS retrievals as input. 
No distinction was made between clear and cloudy retrievals. 

Values are in Wm'^. This figure is limited to the third IFC 
(August 10 to August 21). 

Investigating the quality of the TOVS products 
by comparing them to their radiosonde equivalent, 
we did not detect any bias in the temperature or the 
water vapor content. We therefore suspect that the 
cloud parameters were responsible for the bias. The 
underestimation can possibly result from a lower 
than actual cloud effective emissivity, which is 
equivalent to an underestimation of cloudiness 
and/or cloud emissivity in the TOVS data, or from 
an overestimated cloud base height. The only 
cloud parameter we can compare to surface 
observations is the cloudiness. We attempted to 
perform such a comparison but the study was not 
conclusive because either: 

a) the surface-estimated cloudiness displayed a very 
high temporal variability; 

b) TOVS-based estimates displayed a rather large 
spatial variability: for the same satellite pass, 
TOVS retrievals within 200 km of the area often 
exhibited several octa differences in cloudiness; 
or, 

c) the two types of measurements were not exactly 
comparable: the camera estimate gives a 
fractional cloud cover whereas the satellite 
instrument gives an "effective cloudiness", a 
product of the cloud cover and the mean cloud 
emissivity. 


5* SUMMARY AND CONCLUSIONS 

In this paper we have attempted to compare two 
methods (i.e.. Frouin et al., 1988, and Gupta, 1989) 
for estimating the downwelling longwave flux at 
the surface and assess their respective accuracies. 
Computations with these two models using the 
same atmospheric temperatures and moisture 
profiles indicated, that for clear sky conditions, the 
two model computations are highly correlated. The 
values, however, present a systematic bias of 
15 Wm'^. Since Gupta's parameterization was 
derived from a radiative transfer model applied to a 
set of atmospheric profiles, it is likely that the same 
bias would also be found when this radiative 
transfer model is compared to that of Morcrette. 
Additional work needs be done to understand and 
explain the systematic difference. 

These model computations were then compared 
with in situ measurements at two nearby locations. 
The clear sky longwave flux estimates from both 
models were found to be very consistent with in 
situ measurements. In this case, differences were of 
the same order of magnitude as the uncertainty of 
the measurements themselves. One of our goals^ in 
comparing the two models with in situ 
measurements was to favor one of them, as far as 
the bias was concerned. Unfortunately, the Gupta 
(1989) model results were found to be comparable to 
the station 2 measurements, while the Morcrette 
model results were found to be comparable to those 
of station 38. Since we have no reason to have 
more trust in either dataset, it is impossible, in 
view of this study, to favor one model or the other. 
We hope that further work with the in situ 
measurements at Florida State University will help 
us understand the differences between the two 
stations so that we can achieve more solid 
conclusions regarding the superiority of one model 
versus the other. 

Computations were then performed with the 
two models under all conditions and using TOVS 
data as input. Comparisons with in situ 
measurements indicated that the model 
computations were much more scattered around 
the in situ measurements than were the clear sky 
model estimates with the radiosonde data. This 
was expected since, to the method s uncertainty, 
other uncertainties are added: that of the retrieval 
profile itself and the noise due to the separation 
between the TOVS sounding and the pyrgeometer 
sites. More bothersome, however, is the fact that, 
on average for this small dataset, the flux is 
underestimated by both satellite-based methods. 
No explanation is available for this at present. 

The general conclusion to the present study is 
that, once again, it is difficult to definitively 
quantify the accuracy of radiative transfer models 
because of the lack of reliable surface 
measurements. We have shown, however, that 
the clear sky flux estimates obtained by the two 



models we investigated were within the 
uncertainty of the measurements themselves, 
which is of the order of 5 %, The cloud effect 
computed by the models, unfortunately, cannot be 
validated against in situ measurements because of 
the large uncertainty in the input variables, 
particularly the cloud parameters. 
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1. INTRODUCTION 

It has long been observed that direct solar 
radiation is absorbed substantially by water vapor in 
a cloudless atmosphere, and that the phenomenon 
is even more pronounced when the observations 
are conducted in the infrared part of the solar 
spectrum. Fowle (1912, 1913) first exploited these 
observations for measuring vertically-integrated 
(total) atmospheric water vapor amounts. He 
developed a spectroscopy technique that consists of 
viewing the sun from the earth's surface through 
two spectral channels located near the center and 
outside a near-infrared water vapor band. The ratio 
of the two voltage outputs, as measured by a 
radiometer, gives a direct measure of the water 
vapor amount integrated along the sun-lo-surface 
path. Fowle's differential absorption technique has 
been subsequently verified and applied in many 
studies (e.g.. Hand, 1940; Foskett and Foster, 1943; 
Gates, 1956, Siversten and Solheim, 1975; Pitts et al., 
1977; and Reagan et al, 1987). 

A practical corollary to Fowle's water vapor 
measuring technique, and this paper's focus, is to 
view the earth's surface (i.e., measure the solar 
radiance reflected by the surface) instead of the sun. 
Although such an extension of the technique is 
non-trivial, impeded primarily by uncertainties in 
surface reflectance effects, it is feasible, and its 
successful implementation would be invaluable to 
global earth studies. That is, by viewing the surface 
in place of the sun, this modified technique could 
be conveniently adapted for aircraft and satellite 
measuring systems, thus offering the capability to 
map water vapor amounts over large areas. 

The problem in estimating water vapor 
amounts by viewing the surface in spectral 
channels centered at the different wavelengths, as 
in Fowle’s technique, is that the ratio of the 
channels' voltage outputs depends not only on the 
ratios of the source intensity (the sun) and of the 
atmospheric transmittance (mostly a function of 
water vapor abundance) in the two channels, but 
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also on the ratio of the surface reflectance in the 
two channels. This ratio is not constant; on the 
contrary, it varies strongly with surface type and 
radiation geometry (see, for instance, Bowker et al., 
1985). Consequently, applying the current band 
ratio techtuque to a radiometer viewing the surface 
will yield significantly different water vapor 
amounts depending on the nature and bi- 
directional properties of the underlying surface. 

Our study extends the sun-viewing 
spectroscopy technique to surface viewing for use 
with aircraft or satellite-based instruments. We 
provide appropriate channels in the 0.940 water 
vapor band to retrieve total water vapor amount 
independently of the surface reflectance properties 
and other atmospheric constituents. The 0.940 /jm 
band is selected because it is not saturated, thus 
sensitive to even large water vapor amounts; 
furthermore, the spectral reflectance of vegetation 
in this band is not influenced by leaf water content, 
and soil moisture content generally shifts the entire 
reflectance curve. In other words, liquid water 
absorption by vegetation and soil (also by rocks and 
minerals) is not expected to affect differentially the 
spectral signal reflected by the surface. 

In the following section, section 2, we describe 
our technique to minimize surface reflectance 
effects. We quantify, by performing radiative 
transfer simulations, the gain in accuracy when 
employing the proposed instead of the usual 
channel combination. Section 3 presents a 
verification of the technique using SE-590 
spectrometer measurements and concomitant 
radiosoundings acquired during FIFE. Section 4 
examines, in the case of an airborne or spaceborne 
instrument, the effect of photons backscattered 
directly by the atmosphere toward the sensor, 
which may be a problem over low reflective 
surfaces, such as the ocean in the presence of 
vertically-extended aerosols. Finally, Section 5 
summarizes the findings and concludes with a 
discussion of the technique's perspectives in the 
context of the Earth Observing System (Eos) and 
future international scientific projects. 


2. TECHNIQUE 

To eliminate the surface reflectance effect, one 
may choose two spectral channels, one narrow, the 
other wide, centered on the same wavelength at the 
water vapor absorption maximum in the 0.940 /im 
water vapor band. If the surface reflectance depends 
linearly on wavelength or varies symmetrically 
with respect to the central wavelength, then it 
becomes constant when averaged in the narrow 
and wide spectral channels. In this case, spectral 
variations in the surface reflectance do not affect 
estimates of water vapor amount. Indeed, the 
question arises: does the reflectance of natural and 
non-natural surfaces exhibit a sufficiently 
symmetric or linear behavior in the spectral range 
of the 0.940 water vapor band? 

Using the 5S radiative transfer code (Tanr^ et 
al., 1986), we performed simulations of the radiance 
measured by a radiometer viewing a wide range of 
surface targets (water, snow, ice, rocks, minerals, 
vegetation, wet and dry soil, etc.) whose reflectance 
properties were taken from Bowker et al. (1985). In 
the calculations, the atmospheric path from the 
surface to the sensor was considered negligible, 
corresponding to the radiometer close to the 
surface. The results, presented in Fig. 1, for a 
23 km-visibility US62 atmosphere containing 
continental aerosols, indicate that variations in the 
radiance ratio, due to differences in the surface 
reflectance spectral properties, are substantially 
reduced when choosing the 0.935-0.955 0.920- 

0.970 /im channel combination instead of the 0.945- 
0.955 ^m, 0865-0.875 combination (the typical 
combination when viewing the sun). 
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Fig. 1 Simulations for 140 surface targets, including water, 
snow, ice, bare soil, and dry and wet vegetation, of the radiance 
ratio in sp>ectral channels centered in the region of the 0.940 
water vapor band. Two channel pairs are considered: 0.945- 
0.955 fim, 0.865-0875 fim (crosses) and 0.935-0.955 /im, 0.920- 
0.970 /im (open circles). 


If we now examine the resulting errors in 
water vapor amount (Fig. 2), the standard deviation 
around the actual value (1.4 gcm*^ for the selected 
atmosphere) is reduced from 0.3 to 0.1 gcm"^ (by a 
factor of 3) when using the narrow and wide 
channel combination. The radiance ratio in these 
channels, however, is less sensitive to water vapor 
amount as Fig. 3 demonstrates (the entire depth of 
the band is not measured in that case); but this 
should not be a problem since the instrumental 
noise level in the radiance ratio can be as low as a 
few thousandths. Note, finally that we have only 
considered one model atmosphere. Varying the 
nature and concentration of atmospheric aerosols, 
however, does not change the radiance ratio 
significantly (see Table 1). 
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Fig. 2 Simulations of the retrieved water vapor amount for the 
140 surface targets, sun geometry, and atmospheric conditions of 
Fig. 1. The surface is viewed in channels: 0.945-0.955 and 
0,865-0.875 (a), and channels 0.935-0.955 /cm and 0.920- 

0.970 fjm (b). 
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Fig. 3 SensiUvily of the radiance ratio to water vapor amount 
for channel combinations 0.945-0.955 0.865-0.875 |xm and 

0.935-0.955 /im, 0.920-0.970 /im. 
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Fig. 4 Map of FIFE area showing SE-590 speciromeler and 
radiosonde locations. 


Tabic 1. Radiance ratio for various aerosol models. The values 
are for a sun at zenith and a Lambertian, p>erfeclly reflecting 
surface. 


R ( 0.945 - 0.955 R < 0.935 - 0.955 fm ) 

R ( 0.865 - 0 . 875 fimT R ( 0.920 - 0 . 970 ;im ) 

Atmosphere 


US62,Cont)ncnlaL V=23 km 

.4042 

.3119 

US62,ConlincntaL V=5 km 

.4064 

.3119 

US62.Urban, V=23 km 

.4051 

.3119 

US62, Urban, V=5 km 

.4097 

.3120 

US62,Marilime, V=23 km 

.4035 

.3119 

US62, Maritime, V=5 km 

.4040 

.3119 


Table 2. Summary of concomitant SE-590 spectrometer and 
radiosonde data collected during FIFE. The SE-590 
spectrometer measurements were performed at sites 2, 4, 6, 8, 
10, 16, 36, and 44 and the radiosondes were launched near site 
31 (see Fig. 4 for locations). 


Intensive Field Number of Number of Data Sets 


Campaign 

Field Days 

SE-590 

Soectfomcicf 

Radiosonde 

1 (May/Junc 1987) 

8 

53 

65 

2 (June/July 1987) 

10 

64 

48 

3 (August 1987) 

7 

104 

50 

4 (October 1987) 

7 

55 

48 


3. VERinCATION 

The concept of viewing the surface in narrow 
and wide spectral channels centered on the same 
wavelength to eliminate the effect of spectral 
variations in the surface reflectance has been 
verified using SE-590 spectrometer data and 
concomitant radiosonde observations acquired 
during the 1987 FIFE experiment. Figure 4 shows 
the location of the- SE-590 spectrometer 
measurements and radiosonde launches, and 
Table 2 summarizes the dataset collected during 
FIFE'S four intensive field campaigns (IFC's). 

The SE-590 spectrometer, which measures in 
252 bands from 0.4 to 1.1 /im , collected data at 8 sites. 
At each site the instrument head was installed 
approximately 4.5 m above the ground in a 
relatively uniform and representative area of the 
site. Many of the measurements were made with 
the instrument viewing at nadir, but bi-directional 
datasets were also acquired in the solar principal 
plane during the second, third, and fourth IFC's at 
10° viewing angle increments. To characterize 
diurnal variations in the surface reflectance, the 
data were collected at about every 10° change in 
solar zenith angle. 


A total of 276 SE-590 spectra were acquired 
over 32 field days, and more than 6 radiosondes 
were launched on average per day within 6 km of 
the sites. In the preliminary results presented 
below, however, we only used the data collected at 4 
sites (4, 6, 10, and 36) during the four "golden" days 
(June 6, July 11, August 15, and (October 11) of FIFE. 

Figure 5 shows the radiance ratio in channels 
0.938-0.953 ftm and 0.919-0.972 pm deduced from the 
SE-590 spectrometer measurements as a function of 
water vapor amount along the optical path 
calculated from the radiosonde data. ITie points fit 
very well the law given by: 

R = 0.326 exp[- 0.177(-^)’'^ 

where R is the radiance ratio, U is the total water 
vapor amount, and eis the solar zenith angle. This 
law, where the exponential exponent varies as 
(U/cosO)^^^, is expected since the absorption regime is 
strong in the 0.940 /im water vapor band (see, for 
instance. Goody, 1964). Using (1) yields U/cos 9 to 
±0.5 gem*^, which indicates the good stability of the 
relation between radiance ratio and water vapor 
amount for varied surfaces and a wide range of 





solar zenith angles (the spectral reflectance of the 
sites varies with radiation geometry). Let us 
emphasize, however, that only 4 days of data were 
considered in Fig. 5. More data remain to be 
analyzed (see Table 2) to corroborate statistically the 
above results. 
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Fig. 5 Radiance ratio in channels 0.938-0.953 /iw and 0.919- 
0.972 computed from SE-590 spectrometer data, as a 
function of water vapor amount deduced from concomitant 
radiosonde data. 


4. EFFECT OF AEROSOLS 

So far we have considered the case of an 
instrument viewing the surface from just above the 
ground. In the case of surface viewing from higher 
altitudes or above the atmosphere, the signal 
measured by the instrument contains not only 
photons that have been reflected by the surface, but 
also photons directly backscattered by the 
atmosphere. In the presence of aerosols located 
high in the atmosphere over a low reflecting 
surface such as the ocean, the contribution of these 
photons, mostly backscattered by the aerosol layer, 
may surpass significantly the signal reflected by the 
surface; it is this signal that contains the water 
vapor information (the bulk of the atmospheric 
water vapor is generally encountered within 2 km 
of the surface). The water vapor amount, in that 
case, would be underestimated. 

To assess qualitatively the underestimation 
introduced, let us assume that the aerosols are 
located above the Rayleigh atmosphere. For not- 
too-horizontal sightings, small aerosol optical 
thicknesses, and a Lambertian surface, the solar 
radiance reflected by the surface-atmosphere system 
can be written simply as: 



where and are the signals backscattered by 
the aerosols and molecules, respectively, pis the 
surface reflectance, lo is the extraterrestrial solar 
irradiance, Ta is the diffuse atmospheric 


transmittance, and Tg is the gaseous transmittance. 
In this formulation, we have neglected the effect of 
photons that have sustained multiple surface 
reflections, which is justified for surface targets of 
low reflectance. 

Although extremely simplified, the above 
modeling retains the essential physics of the 
problem. The gaseous transmittance does not affect 
the signal backscattered by the aerosols. When p is 

A 

small, may dominate the right-hand side of 
Eq. (2). As p increases, the relative contribution of 

to R decreases and R becomes more sensitive to 
water vapor amount. 

Fig. 6a shows the radiance ratio in the 0.935- 
0.955 pm and 0.920-0.970 pm channels, computed 
using Eq (2), as a function of surface reflectance. 
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Fig. 6 (a) Radiance ratio in channels 0.935-0.955 pm and 0.920- 
0.970 pm as a function of surface reflectance for 23 and 5 km 
visibility atmospheres. The dashed line corresponds to a pure 
Rayleigh atmosphere, (b) Resulting aerosol-induced error on 
the retrieved water vapor amount. In the calculations, the 
aerosols are assumed to be located above the Rayleigh 
atmosphere, which maximizes the error. 
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When the reflectance is low (p<0.1), the radiance 
ratio departs substantially from the 0.28 value 
corresponding to a pure Rayleigh atmosphere (no 
aerosols), all the more as the visibility decreases. In 
the 5 km visibility case, the difference reaches 0.09 
(32%) for a surface reflectance of 0.01. 

The resulting error in water vapor amount is 
displayed in Fig. 6b. Absolute values belovv 
0.3 gem -2 (21% of the actual value) are only obtained 
when pis above 0.2 and the atmosphere is clear 
(23 km visibility case). For both model 
atmospheres, the underestimation reaches over 
0.9 gcm-2 (63%) when p is less than 0.02. The errors, 
however, should be regarded as upper limits for the 
aerosol backscattering effect. Tropospheric aerosols 
are generally located within the boundary layer, 
whi^ reduces the effect considerably. Radiative 
transfer calculations for realistic atmospheres are 
necessary to quantify the error introduced as a 
function of the relative vertical distribution of 
aerosols and water vapor. 


5. SUMMARY AND CONCLUSIONS 

We have shown, from radiative transfer 
simulations, that viewing the earth's surface in two 
channels, one narrow, the other wide, centered on 
the same wavelength at the water vapor absorption 
maximum near 0.940 pm provides accurate 
atmospheric water vapor amounts independently 
of the surface reflectance properties. The concept 
was verified in situ using concomitant SE-590 
spectrometer and radiosonde data acquired during 
FIFE. Under varied radiation geometries and 
atmospheric conditions, and for several surface 
targets, the relationship between radiance ratio in 
the narrow and wide spectral channels and water 
vapor amount was found to be stable, following an 
expected exponential law. The SE-590 instrument, 
however, did not view the surface from above the 
atmosphere or even above the boundary layer, 
where most of the water vapor is concentrated, and 
therefore did not permit verification of the 
technique in actual conditions of aircraft or satellite 
viewing. In such conditions, simple physics 
showed that the effect of photons backscattered 
directly by the atmospheric aerosols toward the 
sensor is to lower unacceptably the retrieved water 
vapor amount when the aerosols are located at 
high altitudes. 

The technique, however, is only suitable in 
cloudless conditions or, at least, when the sensor 
views a surface target lit by the sun. Still, compared 
to satellite microwave techniques, which are 
applicable under most weather conditions, it has 
the advantage of simplicity and constitutes a 
promising alternative over land, where microwave 
radiometry is not appropriate. 

By providing the means of estimating 
atmospheric total water vapor amounts simply, yet 
accurately from space, our study is relevant to 


international programs, such as International 
Land Surface Climatology Project (ISLSCP), for 
which atmospheric corrections to retrieve land 
surface parameters require estimates of water vapor 
amount, and the planned Global Energy and Water 
Cycle Experiment (GEWEX), whose objectives 
include establishing an observational basis for 
predicting water transport in the atmosphere. In 
this regard, note that the differential absorption 
concept is applicable to Eos platform's High 
Resolution Imaging Spectrometer (HIRIS) and 
Moderate Resolution Imaging Spectrometer 
(MODIS), instruments that possess channels in the 
region of the 0.940 /im water vapor band. 
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abstract 

A new technique is proposed to estimate atmospheric total water vapor amounts from space. The technique 
consists of viewing the Earth’s surface in two spectral channels, one narrow, the other wide, centered on the 
same wavelength at the water vapor absorption maximum near 940 nm. With these characteristics, the ratio of 
the solar radiance measured in the two channels is independent of the surface reflectance and yields a direct 
estimate of the water vapor amount integrated along the optical path. To test the technique, we designed and 
built a two-channel radiometer based on the above concept. Airborne experiments carried out with the new 
device demonstrate the technique’s feasibility under clear sky conditions over both sea and land. Over the ocean 
and in the presence of thick aerosol layers, however, total water vapor amounts may be underestimated by as 
much as 20%. Compared to satellite microwave techniques, which are applicable under most weather conditions, 
the proposed technique has the advantage of simplicity and constitutes a promising alternative over land, where 
microwave radiometry is inappropriate. 


1. Introduction 

Water vapor is an important constituent of the at- 
mosphere. This is manifested in the ability of water 
vapor to change phase within atmospheric pressure and 
temperature ranges, producing clouds and hydrome- 
teors (e.g., rain, snow, hail). When significant amounts 
of water vapor condense, the latent heat release be- 
comes a source of energy for the maintenance of at- 
mospheric processes. Water vapor also affects atmo- 
spheric energetics through radiative interactions. 

Previous efforts to obtain water vapor data on a 
global scale have relied primarily upon radiosondes at 
hundreds of weather stations scattered around the 
world. Remote spectroscopy systems operating from 
space have more recently been deployed, and their ma- 
jor advantage is in obtaining continuous spatial and 
temporal data from not easily accessible regions (e.g., 
over the oceans, deserts, and poles). 

Microwave measurements near the peak of the 
22.235 GHz resonance line from the Nimbus series 
and Seasat have proved very suitable to derive the ver- 
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tically integrated (or total) water vapor amount over 
the oceans under most atmospheric conditions (e.g., 
Staelin et al. 1976; Chang and Wilheit 1979; Grody et 
al. 1980; Prabhakara et al. 1981). Typical accuracies 
of 0.1 to 0.5 g cm"^ have been reported. The results 
obtained over land, however, have not been satisfac- 
tory, mainly because the surface emissivity in the mi- 
crowave spectral region depends strongly on soil type 
and moisture. This strong and variable surface emission 
camouflages the water vapor information in the mea- 
surements. 

Infrared measurements in the 6.3 /xm rotation-vi- 
bration band from Nimbus-6, the NOAA series, GOES- 
5, and GOES-6 have been used to infer the vertical 
distribution of water vapor (e.g., Smith and Woolf 
1976; Smith 1983), but with degraded accuracy, both 
in cloudy conditions and near the surface. The inver- 
sion techniques employed by these authors have typ- 
ically yielded a 30% accuracy in the estimated water 
vapor mixing ratio. 

It has long been observed, however, that direct solar 
radiation is absorbed substantially by water vapor in a 
cloudless atmosphere and that the phenomenon is even 
more pronounced when observations are conducted in 
the infrared part of the solar spectrum. Fowle (1912, 
1913) was the first investigator to exploit these obser- 
vations for measuring atmospheric water vapor 
amounts. He produced laboratory graphs relating the 
opacity of near-infrared water vapor bands to water 
vapor amount. His differential absorption concept has 
been subsequently verified and applied in many studies 
(e.g., Hand 1940; Foskett and Foster 1943; Gates 1956; 
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Siversten and Solheim 1975; Pitts et al. 1977; Reagan 
etal. 1987). 

Fowle’s concept is given further development in the 
present paper. We intend to demonstrate that it can 
be extended to the sensing of total water vapor amounts 
from space. For this purpose, we designed and con- 
structed a radiometer that measures the intensity of 
the solar radiation reflected by the earth’s surface in 
two spectral channels, one narrow, the other wide, 
centered on the absorption peak of the 940 nm water 
vapor band. The ratio of the radiometric signals mea- 
sured by the two channels is independent of the surface 
reflectance properties and yields a direct estimate of 
the water vapor amount along the optical path. We 
report on the first few flights of this new device. 

2. Differential absorption technique 

The differential absorption technique consists of 
viewing a source of radiative energy at two (or more) 
wavelengths through the same atmospheric path; the 
wavelengths are chosen so that the absorption coeffi- 
cients of a given gas, the amount of which is to be 
measured, are different. In the aforementioned inves- 
tigations, the technique was applied to water vapor (the 
gas studied) by viewing the sun (the source) directly 
through the atmosphere. Instead of viewing the sun, 
however, one can view the earth’s surface from above 
the atmosphere to estimate the water vapor amount 
along the optical path. In this case, it is the solar energy 
reflected by the surface that is measured, and this is 
done through a double atmospheric path (sun-to-sur- 
face and surface-to-sensor); but, a priori, one must 
know the surface reflectance. Here we derive a tech- 
nique which requires no a priori knowledge of the sur- 
face reflectance. 

We first consider the case of a direct path between 
the sun and the sensor. Denoting the voltage outputs 
of the radiometer in channels 1 and 2 by Hi and V 2 , 
we have the following proportionality; 

% ( 1 ) 

V2 hh 


along the optical path, U*. Measuring F,/F 2 , therefore, 
gives access to U*. 

We now consider a surface target illuminated by the 
sun. In this case, we have to account for the solar energy 
reflected by the target, which yields 

( 3 ) 

V 2 hhh 

where r, and h are average target jeflectances. If the 
channels are selected such that ri/r 2 is a constant, V\ / 
V 2 remains a function of U* only {Iilh does not de- 
pend on the type of atmosphere encountered), and we 
can still obtain U* by measuring F,/F 2 . 

In the present study, the technique employs two 
channels centered on practically the same wavelength 
at the absorption peak of the 940 nm water vapor band. 
The channels have narrow and wide spectral band- 
widths, respectively. Jhese characteristics, while dif- 
ferentiating between t\ and tz for a fix^d water vapor 
path, allow one to eliminate the ratio rjrz in (3) (ri 
Rs Fz). We are aware that the region around 940 nm 
is influenced not only by water vapor absorption, but 
also, although to a lesser degree, by carbon dioxide, 
ozone, and aerosol absorption as well as molecular and 
aerosol scattering. It is assumed, however, that the 
properties of all the attenuators except water vapor do 
not vary significantly (or vary linearly) across the 
channels’ bandwidth and, thus, cancel in the ratio. This 
assumption is justified, as radiative transfer calculations 
performed with various standard atmospheres dem- 
onstrate (see section 3). 

Note, furthermore, that ( 3 ) neglects the signal back- 
scattered by the atmosphere toward the sensor, which 
fnay not be justified under certain atmospheric con- 
ditions (i.e., thick aerosol layers) when the surface re- 
flectance is small. An analysis of this effect will be pre- 
sented later in the paper, when examining the experi- 
mental results (section 5). 

Let us now express the average atmospheric trans- 
mittance 7 , . For a homogenous path, two basic random 
band models can be used (for more details, see Paltridge 
and Platt 1976): 


where I\ and I 2 are the source intensities and and 
I 2 are average transmission functions. If channel / 
(/ = 1, 2) is characterized -by the spectral response 
li is defined as 


R,{\)d\ 

Jo 

where X is wavelength and t,(X) is the spectral atmo- 
spheric transmittance. When the two channels are lo- 
cated in a spectral region where atmospheric absorption 
is essentially due to water vapor, l^/lz can be expressed 
as a function of an equivalent amount of water vapor 


1 ) the model of Goody ( 1952) 



2) the model of Malkmus ( 1967) 



where di is the average line spacing, a, is the average 
Lorentz half-width, and S/ is the average line intensity. 
In these expressions it is assumed that the spectral in- 
terval considered is wide compared to a,. Since 5, (7 
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> 7 ra, in the water vapor bands above 0.7 ;im, both 
(4) and (5) reduce to 


ti as exp 




( 6 ) 


can be easily taken into account using, for instance, 
Kasten’s ( 1966) approximation formula. 

Using (10), the ratio of atmospheric transmittances 
in the two channels takes the final form 


For a nonhomogenous path, the pressure and temper- 
ature variations along the path {S, depends on tem- 
perature and or, depends on temperature and pressure) 
can be taken into account by scaling U appropriately. 
A one-parameter scaling approximation (e.g.. Goody 
1964) is sufficient, since the absorption regime is strong 
and the effect of tropospheric temperature changes on 
Si is small. The scaled amount then takes the usual 
form: 


i^«exp[-^"'(/«*C/o)"^] 

h 


( 11 ) 


V' 


=r 


_p 

p* 


T* 

T 


dU 


(7) 


where T* and P* are the temperature and pressure of 
the equivalent homogenous path, respectively, and 
is the water vapor amount integrated along the path. 
In the major water vapor absorption bands of the 
shortwave solar spectrum, including the 940 nm band, 
m = 0.9 to 1 and n = 0.45 (Selby et al. 1978; Stephens 
1984). The procedure to calculate l, is therefore to re- 
place U with U* in (6) and evaluate S, at temperature 
T* and a, at temperature T* and pressure P*. 

Using the scaling approximation for a vertical at- 
mospheric path characterized by a temperature profile 
j = a pressure profile P = Poe and a 

water vapor density profile p = poC , and taking 
r* = To and P* = Po yields 


,c( 


I 1 + 




Hp 


Hr 


-I 


( 8 ) 


Typically, Hp ~ 30 km, Hp ^ 8 km, and Hw 2 
km, which gives: 

u* =« 0.8po//(c = 0.8(/o (9) 

where Uo is the vertically integrated water vapor 
amount. Thus, U* is relatively close to Uq. The factor 
relating U* to Uo, how'ever, is not constant and its 
changes result mainly from the variability of //n- (the 
eftect of variations in Hr is negligible). Still, the rela- 
tionship is stable enough to deduce accurately Uo from 
U*, as shown quantitatively in section 3. 

The average transmittance 7, over a slant path can 
therefore be expressed as 


: exp[-|d,(/n*Uo) 


l/2l 


( 10 ) 


where m* is the equivalent air mass (1/cosf? when 
viewing the sun at zenith angle d and 1 /cosd + \ fcosd' 
when viewing the surface a_t zenith angle d' with the 
sun at zenith angle 0) and /J, is an average absorption 
coefficient. Note that atmospheric refraction cannot 
be neglected at high solar or viewing zenith angles {0, 
d' >80°). The resulting increase in air mass, however. 


where /3' = |3i - By measuring 7, /t 2 and knowing 
and the radiation geometry, it is therefore possible, 
at least in principle, to derive Uo from (11). 

3. Radiometer 

A schematic description of the radiometer is given 
in Fig. 1, and Table 1 lists the sensor specifications. 
The radiance passing through the collector is focused 
on one extremity of an optical fiber. Two condenser 
lenses form an image of the other extremity of the fiber 
on the detectors. Using fiber optics was appropriate 
because the instrument was designed to be flown on 
an aircraft. The detectors are silicon photodiodes op- 
erating at ambient temperature. Two interference filters 
are mounted in front of the detectors on a rotating 
wheel driven by a synchronous motor. The spectral 
response of the filters is shown in Fig. 2. The center 
wavelength is 938 nm for both filters, and their band- 
widths are 13 and 46 nm at half-power points, respec- 
tively. These characteristics for the filters ensure that 
Fi/Fi is practically equal td 1 in all conditions (surface 
type, illumination, and viewing geometry). Also shown 
in Fig. 2 is the spectral response of the detectors, which 
is fairly constant over the wavelength range of interest. 
The temperature dependence of the detector response 
is small, not exceeding 0.1% per degree Kelvin in the 
range of temperatures encountered in the troposphere^ 
In addition to the filters, a reference surface is mounted 
to the wheel so that the detectors view an optically 
black target at each rotation of the wheel. A timing 
device permits control of the radiation signals generated 
by the optical system. The signals from the detectors 
are amplified and converted by a 12-bit analog-to-dig- 
ital processor. Two amplifier sensitivities can be se- 
lected, and are appropriate for measurements over land 
(5, ) and sea (5,), respectively. Once converted m dig- 
ital format, the signals are fed into a data acquisition 
unit. The integration time is 0.05 s for each filter and 
the optical zero. 

The differential technique requires only a relative 
calibration of the radiometric outputs. This calibration 
was carried out in the laboratory by directing the light 
collector toward a diffuse target illuminated by a solar 
simulator. The diffuse target was placed about 30 cm 
from the solar simulator and the light collector was 
installed just behind the diffuse target. Typical labo- 
ratory conditions were 292 K for air temperature and 
80% for humidity. With this experimental setup and 
these conditions, atmospheric absorption along the 
optical path was negligible. The procedure was repeated 
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Fig . 1 . Radiometer schematic description. 


before and after each in situ experiment. The ratio V, / 
V 2 was observed to vary little from one calibration to 
the next, by no more than a few tenths of a percent 
around the value 1/0.775. 

Figure 3 gives the ratio of the transmission functions, 
7, and 72, for the narrow and wide band channels, re- 
spectively, computed with the 5S code (Tanre et al. 
1985, 1986) for different air masses and atmosphenc 
conditions specified in the code. The salient features 


f the 5S code are given in the Appendix. In the cal- 
ulations, the spectral atmospheric transmittance was 
onvoluted with both the spectral response of the in- 
erference filters and the spectral response of the de- 
ectors. We see that h/t 2 fits fairly well with the law 
iven by ( 1 1 ) and the coefficient deduced from the 
imulation is 0.178 g'*'" cm^ For a fixed water vapor 
imount, the dependence of ti/h on atmosphere type 
s negligible. The curve in Fig. 3, which was obtained 


Parameter 


Table I . Radiometer characteristics. 


Value/description 


Design Parameters 


Wavelength al half-power points 
Instantaneous field-of-view (total) 

Collecting aperture diameter 

Rotating wheel (supporting the optical filter) 

Detector type 

Detector operating temperature 
Amplifier sensitivity 
Dynamic range 

Narrow-band channel 
Broad-band channel 
Integration time 
Cadence of measurements 
Signal quantizing levels 


927-944 nm (narrow), 
3 ° (instrument viewing 


914-959 nm (wide) 

the sun); 1 1 . 5 ® (instrument viewing the surface) 


3 cm 

3 positions (2 filters -I- a zero); 2 rps 
photodiodes EG&G (type: HUV 4000) 


ambient 
S,: 0.12 W m"^ 


sr“VV; S 3 : 0.46 W m ^ sr ‘/V 


0-0.49 W m ^ sr ' 

0-1.84 W m'^ sr‘ 

0.05 s 

I measurement (2 filters + zero) every 0.5 s 
4096 ( 12 -bit coding) 


Narrow'-band channel 
Broad-band channel 


Noise ecjuivalcnt radiance 

S,: 0.12 X 10'^ W 
S,: 0.46 X 10'^ W 


sr'’; S 3 : 0.42 X 10"^ W m~^ sr ' 
sr'‘‘;S 3 : 1.61 X 10'’ W m'^ sf‘ 


Weight 

Size 

Power (high/low) 


Physical characteristics 
4 kg 

20 X 20 X 24 cm 
70/25 W 
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Fig. 2. Spectral response of the interference filters and detectors. 


by least-square fitting \og{t\lti) as a function o 
departs significantly from the simulated 
points at water vapor amounts above 1 5 g cm .At 
m*Uo = 17 g cm-', for instance, the discrepancy m 
7,/72 reaches 0.014, which translates into a 10% error 
on m*Uo. Small solar and viewing zenith angles are 
therefore favored for greater accuracy in the water va- 
por retrievals. Note aho that the radiometer is suth- 
ciently sensitive io tjtj variations in the 0-20 8 crn 

range ofwater vapor amounts considered in Fig. 

Vj and, hence, are measured to within a few 

thousandths. 


distance between the radiosondes and the radiometer. 
Time differences between radiometer measurements 
and radiosonde launches were accounted for by linearly 
interpolating the radiosonde data. The above-described 
procedure is indeed subject to error, especially since 
some of the radiosonde sites were located several 
hundreds of kilometers from the expenment sites. Most 
of the radiometer data, however, were acquired in Nia- 
mey where radiosondes are launched on a regular basis. 
Processing the radiosonde data revealed that total water 
vapor amounts ranged from 1 to 2 g cm ' ^ This range 
of Uo values is small (in a tropical atmosphere, Uq can 
exceed 4 g cm"'), but according to theory, f, //2 is gov- 
erned by m*Uo and not atmosphere type (see Fig. J). 
Therefore, calibrating the instrument in the range 1- 

2 g cm"' for Uq (using m* from 1 to 5) should also 

be valid for C/o> 2 gem-'. , u * • oic 

Figure 4 shows the ratio of the intercalibrated signals 
generated in the narrow and wide spectral channels as 
a function of m*C/o, from which can be deduced by 
regression. This was done as for the simulations in Fig, 

3 Table 2 gives the values of obtained for each day 
of measurements. The overall mean value of J is 0 85 
± 0. 14 ( l<y) cm. The dispersion of /3 , 7.5% ( l<r) 
of the mean value, is due largely to uncertainties in the 
radiosonde total water vapor amounts [ 1 0%- 1 5% errore 
are frequently reported; for instance, see Richner and 
Phillips (1982)]. Also, the mean value of P corre- 
sponds fairly well to the value predicted by the 5S code 
(0. 178 g-'^' cm) within the accuracy of this code and 
experimental uncertainties. We conclude that the two- 
channel radiometer may derive total water vapor 
amount with a 1 5% ( 1 a) accuracy ( a 7.5% uncertainty 
on 0' yields at 15% uncertainty on m*L/o) using ( l U, 


4. Measurements 

a. Ground-based measurements 

The two-channel radiometer was used in a sun- 
viewing configuration from the ground to validate (11) 
and derive experimentally the coefficient . Measure- 
ments were made at several locations and vanous dates 
in France and Niger: Dinard (48.38 N, 2.03 W) on 
^0 June 1980; Lille (50.39°N, 3.05°E) on 13 and 1> 
May 1980; Roscoff(48.43°N, 

ber 1980; and Niamey (13.32 N, 2.05 E) on 19, 20, 
2 1 23, 24, 26, and 29 November and 1 , 2, 3, 4, 5, and 
7 December 1980. For each day, the data were collected 
at regular time intervals when il < 80° ( 1 < m < 5). 
Radiosonde observations were available for N^iamey 
(one launch per day at 1200 UTC), but not for the 
other locations. In order to estimate Vq at those loca- 
tions, we used 0000 and 1200 UTC radiosonde obser- 
vations at the nearest launch sites, namely Brest 
f48 45°N 4 4l°W), Trappes (48.76°N, 2.01 E),and 
Camborn; (50.21°N, 5,31°W) for Dinard and Ros- 
coff, andUccle (50.80°N, 4.35°E), Crawley (51.08 N 
0.2 1 °W), and Trappes for Lille, and weight-averaged 
the computed total water vapor amounts according to 



0 5 10 15 20 

WATER VAPOR AMOUNT m'U^ (gem ^) 

Fig. 3. 5S simulations of the transmission function i jh for dif- 
ferent water vapor amounts and atmosphere types ( tropica , ’ 

itude summer, midlatitude winter, subarctic summer, subarctic win- 

ter). 
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Fig 4 Ground-based measurements in Lille, Dinard, Roscoff. and Niamey (radiometer viewing the sun): 
calibrated ratio of voltage outputs K, and Kj as a function of water vapor amount m Uo- 


with empirically fitted as 0.185 cm. It is im- 
portant to emphasize, however, that the 1 5% accuracy 
indicated above corresponds to one standard deviation 
in /?' and is therefore valid at the 68% confidence level 
if we assume a normal distribution of fn*Uo values for 
any V 1 /V 2 measurement. 

b. Airborne measurements 

In order to demonstrate the validity of deriving total 
water vapor amounts from space by measuring the solar 
radiation reflected by the Earth’s surface, we installed 
the two-channel radiometer aboard a small aircraft with 
the light collector viewing the Earth’s surface at nadir. 
Two flights were made on 16 May 1979 and 22 May 
1980 over sea and land in the northern part of France. 
Figure 5 shows the aircraft flight pattern for each day 
as well as the nearest radiosonde launch sites (Uccle, 
Trappes, and Crawley), and Table 3 displays the total 
water vapor amounts at these sites before and after the 
flights. The aircraft flew over varied surfaces (forest, 
crops, bare soil, and ocean) whose altitudes did not 
exceed 100 m above mean sea level. 


Figures 6a, 6b, and 6c present the data acquired dur- 
ing the flight of 16 May 1979 at 2800 m altitude over 
land. The two channels’ voltage outputs are plotted as 
a function of time in Figs. 6a and 6b, respectively, and 
their calibrated ratio is shown in Fig. 6c. It is striking 
that each channel output varies rapidly with time, 
echoing changes in the surface reflectance, while the 
ratio of these outputs is very independent of the surface 
properties. This proves quite well that using two chan- 
nels centered on the same wavelength allows one to 
eliminate surface reflectance effects. The ratio in Fig. 
6c, however, exhibits a few abnormal values. Due to a 
malfunction of the timing device, the measurement in 
the narrow band channel was erratically skipped. The 
problem was corrected prior to the second flight on 22 
May 1980. 

Figures 7a, 7b, and 7c show the data record for 22 
May 1980 when flying at 900 m altitude over land and 
sea. Again the ratio of the two voltage outputs (Fig. 
6c) is quite independent of the surface properties, even 
when passing from land to sea, and the technique also 
worked when measuring the small radiation energy re- 
flected by the sea surface. 



454 


JOURNAL OF applied METEOROLOGY 


VOLUME 29 


Table 2. Best linear fit between calibrated Vi/Vj and exp[ ^ X 
for the iffeicnt ground-bas^ “Jf 
viewing the sun). The mean <f' value is 0.185 ± 0.014 g cm. _ 


Table 3. Radiosonde data used to compute total water vapor 
amounts along the aircraft flight paths. 


Location 

Dirard, France 
(48.38^N, 
2.03^W) 

Lille, France 
(50.39‘^N, 
3.05"E) 

Roscoff, France 
(48.43°N, 
3.59°W) 

Niamey, Niger 
(13.32“N, 
2.05‘^E) 


Date 

(d/mo/yr) 


0’ 

cm) 


Residual 

error 


Station 


Date 

(d/mo/yr) 


Time 

(UTC) 


Water vapor amount Uo 
Cg cm'^) 


22/6/80 

0.175 

13/5/80 

0.117 

15/5/80 

0.165 

17/9/80 

0.189 

19/11/80 

0.174 

20/11/80 

0.167 

21/11/80 

0.207 

23/11/80 

0.174 

24/11/80 

0.184 

26/1 1/80 

0.176 

29/11/80 

0.207 

01/12/80 

0.210 

02/12/80 

0.182 

03/12/80 

0.194 

04/12/80 

0.202 

05/12/80 

0.194 

07/12/80 

0.175 


0.008 


0.009 

0.009 

0.005 


0.015 
0.018 
0.020 
0.013 
0.012 
0.009 
0.006 
0.003 
O.OlO 
0.002 
0.002 
0.00 1 
0.003 


Crawley 


Trappes 


Ucclc 


16/5/79 

16/5/79 

22/5/80 

22/5/80 

16/5/79 

16/5/79 

22/5/80 

22/5/80 

16/5/79 

16/5/79 

22/5/80 

22/5/80 


00 

12 

00 

12 

00 

12 

00 

12 

00 

12 

00 

12 


2.01 

1.94 
1.78 

1.77 

1,96 

1.89 

2.27 

1.95 

1.87 

2.46 

2.05 

1.78 


- When flying at such relatively low altitudes, the water 
vapor amount along the optical path is not m*Uo, with 
m* defined in section 2, but 


U = Uo/cose + U'o/cosd' 


( 12 ) 



where U'o is the water vapor amount along the path 
from the surface to the flying altitude. Figures 8a and 
8b give f/, as a function of time for each flight, deduced 
from the radiometer measurements and computed 
from radiosonde observations. The radiosonde data 
were weight-averaged according to the distance between 
the launch sites and the aircraft location dunhg the 
flights, taking into account differences between the al- 
titude of the launch sites and the altitude of the surface 
viewed by the radiometer. The procedure was per- 
formed for times preceding and following the flights, 
and the resulting water vapor amounts were interpo- 
lated linearly with time. As seen in the figures the 
agreement between the two types of f^^^^^mates s 
fairly good. For the flight of 16 May 1979 (Fig. 8a) 
the average value of Ue obtained by the differential 
method is 5.2 g cm-^ and compares with 4_5 g cm 
from the radiosonde data. For the flight of 22_May 
1980 (Fig. 8b), the values are 2.2 and 2.6 g cm . re- 
spectively. These results, however, are not conclusive, 
they are based on only two flights. Furthermore un- 
certainties in the water vapor amounts derived from 
the radiosonde data may be largely responsible for the 
discrepancies. A definitive assessment of the method 
will require more measurements. Perhaps an optimuni 
means of verification would be to install the instrumen 
aboard a helicopter flying above the planetary bounda^ 
layer (where most of the water vapor is concentrated 
and over surface areas where concurrent high-quali y 
radiosonde observations are made as well as standard 
sunphotometer measurements. 

5. Aerosol contribution to measurement errors 

A possible limitation of the method over the ocean 
should be pointed out. At the very low level of the 
radiation signal reflected by the sea surface, the aerosol 
scattering contribution cannot be ignored and may leaa 
to an underestimation of the water vapor amount. Let 
us assume, for the sake of simplification, that the sur- 











TIME (GMT) 

Fig. 7. As in Fig. 5, but for the flight of 22 May 1980 
at 900 ra MSL altitude over land and sea. 
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(15) 


where U{z oo ) is the vertically integrated water 
amount between altitude z and the top of the atm ■ 

"’oTveloping the right-hand side of (14) to first order 
in WZ, ®) and taking Wz, <o) - Vo exp(-z/Hw) 
yields 


z 

2Hw 


(16) 


lOh 


z 6 

O 
2 

^ 4 

□£ 

O 
0 - 
< 

^ 2 

cn 

U1 
1 — 

< 

5 0 


5/22/80 


. radiometer 

o RADIOSONDE 


^ o-o-« e-0-« o~0« 0-0 ^0-0 

0-0 -e 0-0 


7,(z,oo)-l-/3t(«^*^o)'^^exp^- 

Although this approximation for7,(z, oo) is not 
Cdefittion, w" ter vapor absorption .s strong m 
940 nm band), it is sufficient to provide the order ot 
magnitude of the aerosol backscattenng effect on the 

differential method’s acc^acy. 

Substituting (15) and ( >6) ‘"to ( ‘3) one can easily 
evaluate the integral on the nght-hand side o ( ), 

which gives 

Pa'^ A li 

Li ^ <^QA 

IHw \ 

Ha + HIw) 

Since the ratio of the radiometer voltage outputs in the 
namow and wide spectral channels ,s proportional to 
Li/Lz, we have 


' 4 cost!' IT 
X^l -Mm*f7o)‘'^(- 


(17) 


TK70““"nr^8r‘^''90 12 .oT 

TIME (GMT) 

rrp^x'niro“^ 

creases^^exponentially with altitude (optical thickness 
r StSg phase function P., single scattenng at- 
j 1 h^iaht II ^ Let US further assume that 

ra^srt aCi^^^^^^ va^ sntoo.hlv with 

laverngth, so‘^ that they tan “7'>*S„gTe S 
the narrow and wide spectral channels. Using the hrst 

order scattering approximation, the 

scatters toward the sensor in channel i the radiance 
(e.g-, Deschamps et al. 1983): 

, r7(- (13) 


Vx 

Vi 



I - j3i m*Uoi 

\ h ^ 1 


(18) 


Thus, the total water vapor amount detected is 


/ iHw y 


(19) 


dz 


In this expression, 7,(z, co) and dT.,{z)ldz can be 
writteh: 

7 ,(z,co)«exp{-A[m*L/(z,co)]'^^} (14) 


Taking tvpical values of 0.5 km (e_g., 
ai™80) and 2 km (e.g.. Roll '^ 5 ) for ft and 
rcsnectivelv over the ocean, we obtain loa~'J- 
S inJK'ales iha, ihe relriev^ lo.^ water vajmr 

amount is underestimated by «• . q 

ff « Hw the error becomes unacceptable, une nas 
to Doint out however, that (19) gives a supenor error 
Um^r As soon as the surface reflectance reaches a few 
percent the contribution to L, of photons reflected by 
Erface generally surpasses that of photons uniquely 

backscattered by the atmosphere. .ugerved on 

In fact no noticeable decrease in U. is observed on 

22 May 1980 when passing from land ^ 

1 1 9 UTC (Fig 8b). On the contrary, Ve ‘r'creaseo 
rom 2 0 to 2 gcm-^ This change may be attnbu^d 
X amplifier sensitivity S 3 used over land which 
no adapted to measurement over the ocean. After 
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switching to the more appropriate sensitivity S [ , t/, 
quickly recovered the value of 2.0 g cm observed 
over land. 


6. Summary and concluding remarks 


Bloomfield for editing suggestions, and R. Markworth 
for typing the manuscript. One of us, R. Frouin, vvno 
piloted the plane during the field experiments, wishes 
to thank E. Maes for his assistance in pre- and posttlignt 
operations at the Lille-Bondues airport. 

appendix 


From the ground-based measurements, we have first 
shown that the differential technique can be used to 
derive total water vapor amounts by viewing the sun 
through the atmosphere in two channels of different 
spectral bandwidths centered on the same wavelength 
near 940 nm. The experiment was calibrated by com- 
paring the ratio of the radiometric outputs generated 
for the two channels with water vapor amounts de- 
duced from radiosonde observations. The relative ac- 
curacy of the water vapor retrievals is 15% ( U). This 
value, however, incorporates uncertainties in the ra- 
diosonde data (including space and time interpolation 
errors), so that the actual accuracy is likely to be better, 

probably around 5%. . 

The two-channel radiometer was then used in an 
airborne configuration, measuring the solar radiation 
reflected by the Earth’s surface. In this configuration, 
the instrument was able to yield a measure of the water 
vapor amount along the optical path independent of 
the surface nature, even for the lower radiation signals 
reflected by the sea surface. In the presence of thick 
aerosol layers, however, the water vapor amounts de- 
rived over the ocean may be underestimated by as 
much as 20%. One way to remove this limitation, at 
least partially, is to view the sea surface in the sun glint, 
which substantially increases the contribution of the 
signal reflected by the sea surface. Additionally, viewing 
the sea surface inside and outside the sunglint region 
would provide information on the ratio of aerosol and 
water vapor scale heights by giving access to Lq and 

(7oa [see Eq. ( 19)]. • , u 

We conclude by suggesting that the differential ab- 
sorption technique presented herein can be applied to 
yield accurate space observ ations of total water vapor 
amounts under clear sky conditions over land and sea. 
Compared to satellite microwave techniques, which are 
preferred over the ocean since they are applicable in 
almost all weather conditions, our technique has the 
advantage of simplicity and would complement the 
microwave techniques over land where they fail. In- 
terestingly, the Earth Observing System of the 1990s 
will carry the High Resolution Imaging Spectrometer 
(HIRIS), an instrument that possesses adequate chan- 
nels to exploit our differential absorption concept. 
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surements, C. Deroo for her programming support, B. 


The 5S Code 

The Simulation of the Satellite Signal m the Solar 
Spectrum (5S) code (Tanre et al. 1985, 1986) com- 
putes the solar radiation backscattered to space by the 
Earth-atmosphere system as it may be observed by a 
satellite sensor. Given a Lambertian ground target, the 
code estimates the target’s apparent reflectance by 
taking into account the effects of gaseous absorption, 
scattering by molecules and aerosols, and spatial in- 
homogeneities in the surface reflectance. The input 
parameters, namely, solar and viewing geome nes, 
atmosphere model, surface reflectance, and spectra 
band, can either be specified from standard conditions 
or user-defined. In addition to apparent reflectance, the 
code provides gaseous transmittance and irradiance at 
the surface, as well as the various components of the 
satellite signal. Complementary results are also avail- 
able; exact calculations at selected wavelengths, in par- 
ticular, allow one to assess the code accuracy. 

Based on Tanre et al. ( 1979), the satellite signal is 
expressed as a function of the successive orders of ra- 
diation interactions in the coupled surface-atmosphere 
system. If p is the reflectance of the target, and p. that 
of its environment, the apparent reflectance is wntten 

as 


6\ <j>) = tg{6, e')lpa(0, <l>) 


+ 


[g-WcosO ^ /,/P''|]| (Al) 


I pgS 


vhere 8 and O' are the sun and satellite zenith angles, 
•espectively, </> the relative azimuth between sun and 
atellite directions, r the atmospheric optical thickness, 
'the gaseous transmittance, /^the atmospheric diffuse 
transmittance, and S the spherical albedo of the at- 
mosphere. The first term enclosed by the curly brackets 
represents the contribution of photons backscattered 
to space without surface reflection, whereas the second 
term characterizes photons that have sustained one or 
multiple surface reflections. Absorption by atmosphenc 
gases is considered as a single multiplicative factm de- 
pendent on the direct paths sun-to-surface and surtace- 
to-sensor. Decoupling absorption and scattering pro- 
cesses are justified since, on the one hand, 
located at altitudes where molecules are rarefied, and 
on the other, water vapor and carbon dioxide absorp- 
tion occur above 850 nm where molecular scattenng 
is negligible, and first and second orders of aerosol 
scattering (predominantly forward) restitute almost all 
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of the diffuse radiation. According to Tanre et al. 
(1986), the error introduced by separating the two 
processes is smaller than one percent, except for grazing 
incidence or observation directions (cos0, cosO <c 0. 1 ). 

The atmospheric functions td and S are approxi- 
mated by analytical formulas determined empirically 
from exact radiative transfer computations performed 
for a wide range of model atmospheres. Table A I shows 
the disparity between 5S and exact calculations of 
the total atmospheric diffuse transmittance, e 
+ various solar zenith angles and wave- 

lengths. Calculations were made for two atmospheres, 
clear and hazy. The differences are small, generally less 
than 1%, but may reach over 2% (hazy atmosphere, 

0 = 60 °). 

The gaseous transmittance, tg, is computed from two 
exponential random band models, that of Goody 
( 1964) for water vapor, and of Malkmus ( 1967) for 
oxygen, ozone, and carbon dioxide. The spectral res- 
olution, 20 cm”*, is sufficient (contains enough spectral 
lines) to apply the random band models confidently. 
Figure Al compares the gaseous transmittance in the 
spectral region of the 940 nm water vapor band com- 
puted using the 5S code and with a well-known code, 
LOWTRAN-6 (Kneizys et al. 1983). For the tropical 
and midlatitude summer atmospheres considered, the 
agreement is good near the peak water of vapor ab- 
sorption, but notable differences exist in the wings of 
the band. We recall here that for computational effi- 
ciency, LOWTRAN-6, unlike 5S, neglects the influence 
of temperature on the molecular absorption coefficient, 
and approximates molecular line absorption by a one 
parameter band model. Nevertheless, integration over 
the wavelength range of the narrow and wide band- 
passes considered in the present study provides very 
similar results with both codes. 


Table A- 1. Comparison between 5S and exact calculations of the 
total atmospheric diffuse transmittance. 


0 = 15° 60° 


Wavelength (nm) 


450 

550 

650 

850 

1600 

2200 


450 

550 

650 

850 

1600 

2200 


5S Exact 


Clear atmosphere 

.877 

.873 

.930 

.928 

.954 

.953 

.973 

.973 

.988 

.989 

.992 

.993 

ffazy atmosphere 

.810 

.806 

.866 

.864 

.895 

.895 

.926 

.927 

.963 

.964 

.973 

,976 


5S Exact 


.755 .753 

,841 .839 

.885 .883 

.926 .924 

.966 .963 

.976 .975 


.648 .633 

.721 .705 

.766 .750 

.819 .806 

.899 .890 

.925 .921 




Fig a- 1 Gaseous transmittance along a vertical path sun-to-sur- 
face computed with 5S and LOWTRAN-6 for tropical and midlat- 
itude winter atmospheres. 
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ABSTRACT 

eSS=sf=?:.rH-3 

satellites The probie f^^etions of radiance or reHectance 

secon , ^ indices not proportional to primary 

reflectances at optimum (one or two) viewing zenith ang . 

Keywords; radiation, plants, primary production. 

1, introduction 

Oof n RnsicaUv the major portion of human food is 
pro^dU by plants groiving over 

also affects the environmental context m which man and soci^ies 
de^etop in addition to its key role 

Uls'nnc and structurine communities, land primary 

at^spheri^ carbon dioxide and impacts the water 
available to the atmosphere. Anthropogenic changes, such 
thosl hnked to the destruction of major vegetation systen^, have 

;,nd their effects on climate, we need to know the ^ 

distribution and temporal variability of land primary pro uc iv y 

over the globe. 

To achieve this goal, satellite observations are essentiah A 
promising technique for sensing 

least in the case of light-limited situations, incorporates the tact 

that the growth rate of many plants >« P'X'fXge (Ref. 

rate at which radiant solar energy is absorbed by the fohag (, 

2): 

PP=efp^PAR O) 

where PAR is the incident photosynthetically 

oractically the solar radiation reaching the canopy m the . 

K sl^tra? interval, /p,. is the fraction of PAR .ntercepted by 


the canopy,, and e is -il 

tX"1or vttsTan%ies -°Pf ‘ 

Se«nl.IandL4gCperMJofPAR (Ref.3). 

studies demonstrate (e.g., Ref. 9). 

nori^izLl difference, ND. These indices are defined by. 


SR = 




ND = 




R^’^Ry 


(2) 


(3) 


u V r,d Rvf are UD welling radiances in the visible and near- 
where Ry and Rn are upwemi ^ Advanced 

infrared (for instance radiance in channel 1 satellites). 

Very High Resolution Radiometer aboard NOAA sa^n ' 
r^lJctivfly. instead of radiap, 

SR andND are sensitive to results from 

^sm 

correlates with /p^j?- 

s,:. s. n^» .ho~»> -bouM “-““"'r 

, J- / 2 14 15, 16, 17) have also 

Experimental studies (^^^^^ did not resolve the 

provided disparate results d , 2) for instance, 

apparent theoretical controversy. fnr^uJar beet which 

observed that SR is linearly related to 

supported their theoretical analysis, ye e observations 

refo^ed an exponential n f and vegetation 

inHiVAtP that the relationships between fp^R 6 

“aS: di;... s.™.siy «. 'lVP.olc.~py .od ood.,.y,ng 

...rf,rp as nointed out by Choudhury (Ref. 13). 
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Thatlheinlerdependcnoeof/pAjt and ''®8e‘aHon md^s is 
non-linear is an imporianl issue because linearity ‘ 

«me-Wal of SRand ND should also be linearly related to 
primary production (e.g.. Refs. 11, 12, 18). =ven though senous 
doubtZexist about the linearity of the relationships, 
when considering various ecosystems and climates. Tucker et . 
(Ref. 19) and Coward et al. (Refs. 20, 21) have .^Sf^ement 

between ND time-integrals obtained from AVHRR data over the 
Senegalese Sahel and North and South American biomes, 
resp^Hvely, and published production rates. 

Using SR or ND for estimating /p^ , has several limitations. First, 
as mentioned above, the relationships are generally non-lmear, 
which makes SR and ND time-integrals not proportional to 
primary pixxiuction. Second, SR and ND are non-linear functions of 
radiance. Since vegetation is highly heterogeneous spatially, sub- 
pixel variability is likely to introduce uncertainties in SR and ND, 
particularly when the sensor spatial resolution is coarse (case of 
AVHRK). For such sensors, which have the advantage of frequent 
global coverage, applying relationships established for 
homogeneous canopies is not satisfactory. Third, satellite-denved 
SR and ND may reduce to some extent the effects of sensor 
calibration uncertainties and atmospheric interference, yet they 
depend on atmospheric composition, in particular aerosol and water 
vapor amounts, and viewing geometry (e.g.. Ref. 22). Even time 
series of maximum AVHRR vegetation indices over a several-day 
period, which correspond to minimum atmospheric contamination, 
remain relatively noisy, and it has not yet been possible to identify 
whether the noise is due to residual variations in the atmospheric 
contribution or to variable directional surface properties. We need 
to address these limitations if we are to remotely sense fp^j^ from 
space accurately . 


Our objective, therefore, and the purpose of this paper, is to define 
optimum combinations of visible and near-infrared reflectances 
that: a) relate linearly to b) can be used independently of the 
type of foliage and substrate; c) eliminate the effects of sub-pixel 
spatial heterogeneity; and d) improve the accuracy of fp^j( 
estimates when compared to SR and ND. 


2. METHODOLOGY 

Instead of using radiance ratios, we express fp^j^ as a linear 
combination of visible and near-infrared radiances or, 
equivalently, reflectances. This procedure, when applied to a 
coarse resolution sensor such as AVHRR, should eliminate or, at 
least, substantially reduce sub-pixel variability effects. Linear 
combinations of reflectances, known as "greenness’' transformations, 
have been used for many years to study vegetation parameters, in 
particular by Hatfield et at. Ref. 15) and Asrar et al. (Ref. 17). 
These authors found that greenness obtained by combining 
reflectances measured by a Barnes Modular Multispectral 
Radiometer (MMK) in two visible and two near-infrared bands is a 
much more linear predictor of fp^n than simple ratio and 
normalized difference. They did not favor greenness, hovNever, 
because of the smaller sensitivity of this index to 3nd the 

strong dependence of the relationship between greenness 

upon solar zenith angle and canopy geometry. If known (e.g., from 
theoretical calculations), the dependence upon solar angle or, more 
generally, radiation geometry, should not be a problem because 
solar and viewing angles can be determined exactly. The problem is 
to eliminate the effects of variable canopy geometry and soil 
reflectance in the relationships. Given a sun px^sition, this may be 
possible for specific viewing angles. 

Our approach, therefore, is to simulate for varied soil and canopy 
parameters, namely leaf optical properties, soil reflectance, leaf 
area index (LAI), and leaf inclination distribution function (LTDF), 
above-canopy visible and near-infrared reflectances as well as fp/^R 
and daily averaged (weighted by incident radiation) fp^R, 

The simulations arc performed which the SAIL canopy reflectance 
model (Ref. 23). Various radiation geometries (solar and viewing 
zenith angles, relative azimuth angle) are considered, as well as 
direct and diffuse fractions of incident solar radiation. The 
absorbed fractionjof VAKIpar • computed as a function of solar 
zenith angle and as a function of latitude and season. From the 
reflectances, simple ratio and normalized difference arc derived. 

We focus on/^^ since this parameter rather than fp^R is required in 


primary productivity models. In addition since the senso^ 
TOtentially useful to monitor land primary productivity from space 
are, or will be carried by heliosynchronous satellites (AVHR^ on 
the NOAA series, the POUrization and Directionality of the 
Earth Reflectance instrument, POLDER, 

MODerate resolution Imaging Spectrometer, MODIS, on EOS) and, 
therefore, offer the possibility of viewing the same target under 
one or several geometries during a several-day period, we attcrnpl 
to estimate from single or multi-angle combinations of visible 
and near-infrared reHecUnces. Indeed, the multi-angle approach 
is only suitable when the characteristics of the surface target do not 
change significantly over the several-day period. 

Thus, we regress Wat each latitude and month during the year 
against simple ratio, normalized difference, and visible and near- 
infrared reHectances. One and two viewing geometries are 
considered for the combinations of reflectances. The regression 
sUtistics, namely correlation coefficient, regression coefficients, 
and residual error of estimate are analyzed to determine the solar 
and viewing angles that minimize the effects of vanable LIDF and 
soil reflectance. The improvement in the predicting power of the 

,.,.n.V.;r.9fi,.nc ic alcrt 


3. RESLILTS 

To illustrate our theoretical approach and show the promise of 
linear combinations. Figs. 1, 2, and 3 display selected results 
obtained with the SAIL model. The calculations were performed 
for LAIs of 0.01, 0.05, 0.1, 0.2, 0.5, 1, 2, and 5, erectophile, 
spherophile, and planophile canopies, soil reflectances of 0.1, 0.2^^ 

0 3 and 0.4, and typical leaf optical properties (reflectance and 
transmittance of 0.1 and 0.001, respectively, in the visible, and 0.45 
and 0.4 in the near-infrared). The soil reflectance was assumed 
white spectrally, and the various LIDFs were considered conjointly 
in the regressions. A US 62 standard atmosphere (Ref. 24) 
containing continental aerosols (Ref- 25) was overlying the canopy. 
We in Fig. 1 (bottom) that the influence of the background on 
the fpAR versus ND relationships is substantial, especially at 
moderate LAIs, but is reduced considerably when using linear 
combinations of reflectances (Fig. 2, bottom). In this case, the points 
corresponding to a same LAI are generally more aligned with the 
best fit line. The result is a drastic improvement in theW residual 
error. For the solar and viewing geometries of Figs. 1 and 2 ,bottom, 
the residual error is reduced from 0-058 to 0.033. When using ND 
the minimum residual error is obtained for a nadir viewing (Fig. 1, 
top), but when using linear combinations it is preferable to view the 
canopy at a 45° zenith angle (Fig. 2, top). The fp^R residual error 
can be further reduced by combining linearly visible and near- 
infrared reflectances at two viewing zenith angles (Fig. 3). Using 
reflectances at nadir and 60^ from zenith, for instance, reduces the 
/p^ residual error to 0.026 (Fig. 3, bottom). Smaller residual errors 
can even be obtained when the second viewing zenith angle is as far 
' as possible from nadir ( Fig. 3, top). 

The above results, however, are only valid for a sun at 6CP of zenith 
in July and at 45° latitude. For a sun closer to zenith, the minimum 
/pax residual error is encountered at higher viewing zenith angles 
when using ND, for instance at 45® when the sun is at 30® from 
zenith. In the case of uni-angle linear combinations, the minimum 
at 45° (Fig. 3, top) moves to 60° when the sun zenith angle 
decreases to 30°. The picture is more complex with multi-angle 
combinations because of the many angular possibilities. In general, 
for a particular sun configuration^several viewing zenith angle 
pairs provide similar good results ifp^R residual error around 0 020). 
For a sun at 30° from zenith, for instance, viewing at nadir, 15°, or 
30° from zenith and at 75° from zenith gives residual errors ranging 
from 0.019 to 0.021. The regression coefficients, however, are quite 
sensitive to the viewing geometries selected. 


4. DISCUSSION 

The results presented in section 3, although encouraging, should be 
interpreted with caution. No hasty generalization can be made at 
this point. First, the SAIL model has often showed weaknesses 
when compared to measurements; it docs not predict a hot spot and 
is only appropriate for agricultural plants that form a layer-type 
canopy. More accurate canopy reflectance models may be used, al 
least to provide a reference. Second, the background reflectance 
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Fig. 1 Top: residual error as a function of viewing geometry foj^a 

<un at 60° from zenith in July at 45°N. In the calculations, f^g 
'estimated from normalized difference. ND, is compared_to fp^g 
obtained with the SAIL model. Bottom: scatter plot of fp^g versus 
ND for a nadir viewing. Variable soil reflectance results in points 
aligned rather perpendicularly^ to the best fit line WosM line), 
especially af moderate LAIs fp^g and ND v^ues around 0.5), which 
indicates that the relationship between fpj^g and ND is not only 
non-Iincor, fcwl also strongly depends an the type of background. 




0. 074-1. 748pv +1.526p|^ 

Fig. 2 Same as Fig. 1. but linear combinations of msible and near- 
infrared reflectances (single viewing geojnetry), Pv Pn' 

respectively. Compared to Fig.I. the fp,,g residua 

substantially reduced when viewing around 45° from 
corresponding to a same LAI. but different soil reflectances, are now 
more aligned along the best fit line. 


may vary with wavelength in the visible and 
the case with most soils (e.g.. Ref. 27) and leaf litter (Ref 28). The 
canopy may also be composed of living as well as dead leaves or, 
more generally, leaves of different optical properties. 
Calculations, therefore, should include more realistic situations. 
Still, our study strongly suggests that linear combinations at 
specific viewing angles may be rru. h more accurate in predicting 
than indices that are non-iinear functions of radiances or 
reflectances. 


The relationships established theoretically, the gain in fp^R 
residual error when using preferential viewing geometries, the 
ability of the linear combinations to reduce soil and LIDF 
dependence, etc., remain to be verified using in situ measurment^ 
Unfortunately, few data sets exist that contain concomitant Jpar and 
reflectance measurements at various viewing angles. The result^ 
therefore, may not be statistically significant. A 
experiment to establish and verify the validity of the SAlL-based 
data fits, therefore, appears necessary. 


One should further emphasize that surface reflectances observed 
from space are inherently subjected to instrument noise and are 


contaminated by the atmosphere. Consequently, it will be necessap^ 
in the comparisons of the various estimators^ performance to 
the effects of instrument noise and atmospheric interference, which 
act differentially on simple ratio, normalized difference, and 
linear combinations. This can be done by simulating the top of 
atmosphere reflectances corresponding to the surface reflectances, 
correct those renecUnces for atmospheric effects assuming typical 
atmospheric characteristics, and translate the effects o 
uncertainties in these characteristics into above-canopy reflectance 
uncertainties. The procedure is then to introduce the above canopy 
reflectance uncertainties in the regression datasets, as well as 
typical instrument noise. 


r investigation should be viewed in the 
cebome radiometers, in particular MODIS on EOS and POLDER 
nosphere reflectances for atmospheric interference, are 
■ticularly adapted to monitor /par and, therefore, primary 
iductivity from linear (uni- and multi-angle) combinations of 
lectances. This should lead, during the EOS era, to a better 
iraclerization of terrestrial primary productivity on a global 
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Fi^. 3 Same as Fig. 2, but two viewing geometries. Compared to Fig. 
2^ the residual error is further reduced when reflectances at 
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1. INTRODUCTION 

Spurred by the possible consequences of the so- 
called greenhouse effect, or the trapping of 
longwave radiation within the atmosphere by 
water vapor, carbon dioxide, and trace gases, 
public awareness and scientific interest in the 
planet Earth's climate system has grown 
considerably during recent years. 

General Circulation Models (GCMs) have been 
trying to identify pertinent climate changes. 
Predictions of global changes, however, vary 
widely. For a doubling of carbon dioxide 
concentration, for instance, the predicted mean 
temperature increase ranges from about TC to 
4°C. It is now believed that the main factor of 
uncertainty is cloud feedback: changes in the 
Earth's climate may drive variations in cloud 
parameters (height, liquid water content, 
brokenness, distribution) that could amplify or 
cancel out the greenhouse effect. 

Clouds affect climate in many ways; however, 
their main effect regarding climate research, of 
particular interest in this study, is on radiative 
exchanges occurring within the 
surface/atmosphere system. Except for some 
pathological cases of very high surface albedo 
(Kobayashi, 1989), clouds increase the Earth's 
system albedo which, by reducing the solar 
radiation absorbed by the surface, leads to a 
cooling of the Earth. On the other hand, clouds 
trap terrestrial radiation within the troposphere, 
reducing the amount of longwave flux escaping 
to space, and thereby causing a warming of the 
surface. The net cloud radiative forcing is 
therefore a difference between the two effects, 
which are of the same order of magnitude. Its 
sign and amplitude depend on cloud 
characteristics, in particular height and water 
content (Somerville, 1984). 

Current cloud forcing predictions are of poor 
quality for the following reasons: 1) we still 

cannot correctly predict the occurrence of clouds 
under given atmospheric conditions and 2) the 
radiative effect of a cloud, given its li(piid water 


content and height characteristics, is subject to 
uncertainties, even though the radiative 
transfer through a vertically inhomogeneous 

atmosphere is now well understood (Plateau 
and Stephens, 1988). 

Given this situation, it is important to conduct 
observations of cloud radiative forcing and 
understand what parameters drive the main 
variations. Parameterizing cloud forcing from 
variables that are well predicted by GCMs when 
compared to present climate is a key objective. 
The study presented in this paper evaluates a set 
of radiative flux surface measurements, 
concomitant with other meteorological 
observations, so as to derive the cloud radiative 
forcing at the surface. The forcing is then 
analysed in terms of variability and relation to 
other meteorological parameters. 

In the first part of this paper we describe the data 
used in this study, which were collected during 
the First ISLSCP Field Experiment (FIFE). We 
then present the methods that are used to 
process these data and characterize the cloud 
radiative forcing at the surface. This forcing's 
variability is finally examined as a function of 
cloudiness and atmospheric changes. 

2. DATA 

The First ISLSCP Field Experiment (FIFE), an 
international surface-atmosphere experiment, 
took place in 1987 at and around the Konza 
Prairie Long Term Ecological Research (LTER) 
site near Manhattan, Kansas (Sellers et al., 1988). 
The FIFE objectives were to gather the necessary 
data to permit interpretation of satellite 
observations suitable to infer climatologically- 
significant land surface parameters. 

The experimental area was a 15 x 15 km square 
with various topological features including 
burned and unburned plateau, slope, and creek 
(Fig. I). Although the central portion of the area 
is studied continuously and additional 
parameters were measured throughout the year, 



Ihe experiment concentrated on four Intensive 
Field Campaigns (IFCs) corresponding to four 
different states of surface vegetation: 

"green up" May 26th to June 6lh 

"peak greenness" June 25th to July 15th 
"dry down" Aug. 10th to Aug. 21st 

"senescence" Oct. 5th to Oct. I5th 

Our Study uses 4 different types of data acquired 
during the IFCs, namely surface radiative 
measurements, atmospheric soundings, sky 
photographs, and conventional cloud 
observations. 

Several surface stations collected radiative flux 
measurements. We selected stations 2 and 38 
(see Fig. 1) and acquired the radiation data from 
the FIFE information system. These stations 
were operated by Eric Smith's group (Florida 
State University) during the four IFCs. The 
surface measurements were compared with 
other stations and with radiative models 
outputs using atmospheric profiles (Br4on et al, 
1990a). A bias was found in the longwave 
measurements, but was later corrected (E. 
Smith, personal communication), so that we 
can expect an accuracy at least equal to the 
manufacturer's specifications (within a few 
percentile). 


First IFC : 
Second IFC ; 
Third IFC : 
Fourth IFC : 



“ *5 ^ 

Fig. 1: General presentation of the FIFE 

experiment field. Station 2 and 38 are the two 
surface radiation measurement sites used in this 
study. The black star indicates the locatioh of the 
soundings. 

Fractional cloud cover was another parameter 
observed from the surface. Within the FIFE 
area, observations were made with an 
uplooking automatic camera. Cloudiness was 
deduced from the percentage of cloudy areas in 
the photographs, taking into account the 
viewing angle of each area relative to the 
vertical. Since the camera was operating in the 


visible spectrum, no cloudiness information 
could be obtained at night. Another limitation 
of the camera is that no information on cloud 
thickness or height could be provided. 

In the immediate vicinity of the FIFE area, two 
stations (MHK and FRI) collected standard 
meteorological data. We acquired their cloud 
observations through the FIFE information 
system. These observations consisted of hourly 
estimates of three, or less, cloud layers giving 
their height and fractional cover in octas. In 
order to use the data in our particular study, we 
classified them in three atmospheric layers: low, 
mid-level, and high, with limits at 900 and 3000 
meters. For each layer and each day, a mean 
cloudiness was computed, assuming a random 
overlap if several cloud layers were found 
within the same atmospheric layer. Fig. 2 
presents a comparison of the mean cloudiness 
as found during the 4 IFCs by the two stations. 



Fig. 2: Comparisons between the daily cloud 

amount estimates from two nearby 
meteorological stations for the three cloud layers 
(less than 900 meters, between 900 and 3,000 
meters, more than 3,000 meters). Values are 
given in octas. 

The scatter is rather large, even though the two 
stations are close to one another (about 10 
kilometers) and are located in a similar surface 
environment. It is interesting to note, however, 
that the correlation between the two estimated 
cloudiness datasets, increases from low to high 
clouds. With the levels chosen, mean 
cloudiness results mostly from high clouds 
(altitude greater than 3000 m) that can produce 
overcast skies on a daily average. Most low 
cloudiness values are less than 2 octas but 
reached 5.5 octas on one occasion. 

The atmospheric measurements used in our 



study were obtained from radiosondes launched 
from the experimental site (see Fig. 1). These 
measured pressure, temperature, and wet bulb 
temperature, from which the water vapor 
mixing ratio was derived. Pre-processing and 
extrapolation of the profiles is described in 
Br^on et al. (1990a). 


determine the cloud effect. The longwave effect 
is simply the difference between the 
measurement and the ''clear" flux estimate. 
The shortwave estimate is defined differently 
because we have to compensate for the large 
diurnal cycle of the shortwave irradiance. We 
define the shortwave cloud effect by the 
formula: 


3. METHOD 

To derive the cloud radiative forcing from the 
surface observations, one must first estimate the 
downwelling irradiance in the absence of 
clouds. The two spectral intervals, shortwave 
and longwave, need different methods. 

The "clear" longwave flux was computed using 
Morcrette's radiative transfer model (Moraette, 
1984; Morcrette and Fouquart, 1985; Morcrette et 
al., 1986), with the FIFE soundings of 
temperature and water vapor mixing ratio as 
input. Ancillary data, namely ozone profile, 
carbon dioxide concentration, and temperature 
and water vapor profiles in the upper layers of 
the atmosphere, had to be set to their 
climatological value since no other information 
was available. Although necessary, these data 
have little effect on the results. After 
comparison with other radiative transfer 
models (WMO, 1984) and validation with in- 
situ observations (Breon et al., 1990b), it is 
believed that accuracy of models such as 
Morcrette's is better than 10 W m 2. During the 
FIFE experiment, there was an average of 6 
soundings per day, during day-time, which 
provided a time step sufficient to follow the 
longwave flux variability. "Clear" lonewave 
fluxes could therefore be interpolated from one 
sounding to the next. 

Main variations of the shortwave irradiance at 
the surface are driven by the solar zenith angle 
and cloud parameters (amount, optical 
thickness). Other parameters, such as the water 
vapor profile and the ozone content, have a 
second order effect and can be set to Hheir 
climatological value. Effects of aerosols, 
although small in many cases, can be large. 
During FIFE , aerosol turbidity measurements 
were too sporadic to establish their variability 
since they were taken only during clear days. 
We therefore estimated an ad-hoc atmospheric 
visibility for each of the IFCs by adjusting our 
shortwave irradiance estimates to the 
measurements obtained during clear periods, as 
seen by the clear-sky camera. Computations 
were done with an half hour time step, using a 
parameterization developed by Frouin et al. 
(1988) from the radiative transfer model of 
Tanr^ et al., (1979). 

From the "clear" flux values obtained by the two 
radiative transfer models, it was easy to 
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where SV/obs Is the instantaneous shortwave 
flux measured at the station, and SY/cUar and 
SW are the instantaneous and daily- 
integrated "clear" shortwave flux as estimated 
by the model, respectively. This formulation 
allows one to cancel out most of the diurnal 
variability for a constant cloudiness. It is 
obviously not usable for night period, when 
clear is equal to zero. When doing time 
averages, however, we did not simply average 
CLeff but gave each instantaneous value a 
weight proportional to the shortwave "clear" 
flux estimate. This method allows one, in case 
of variable fractional cover and non-continuous 
measurement, to obtain the statistically best 
possible estimate of a time average. 


Net cloud radiative forcing at the surface is the 
difference between the longwave and shortwave 
cloud effects. On a daily time scale, a negative 
value corresponds to a net cooling of the surface 
(shortwave forcing is larger than the longwave 
forcing); a positive value indicates that 
longwave heating overcomes shortwave 
cooling. 


4. RESULTS AND DISCUSSION 

Because of the notably variable cloudiness, as 
well as the uncertainties on each of the 
parameters, time series of the instantaneous 
cloud radiative forcings are extremely variable, 
especially close to the sunrise and sunset when 
SWc/ear values are small. Fortunately, when 
doing time averages, these values are given a 
low weight and the uncertainty is therefore 
reduced. Similarly, the uncertainty is 
statistically reduced when averaging several 
observations taken at the same local time. We 
therefore chose not to present instantaneous 
values, from which information is difficult to 
distinguish from the noise, but rather the daily 
cycles and the day-averaged values. 

Table 1 gives, for most days of the FIFE 
experiment, a daily average of the various 
measurements and estimates presented above. 
"Mod" refers to the model estimates of clear sky 
radiative fluxes, "obs" to the observations, and 
"cLeff" to the cloud effect. "Nobs'' gives the 
number of half-hourly values that were used in 


Tab. 1: Daily averages of radiative parameters and cloudiness during the four FIFE IFCs. 
SW refers to the shortwave irradiance and LW to the longwave irradiance. Clear sky 
estimates from a model are denoted "mod", observations "obs" and cloud effects 
"cl_eff". Radiative flux parameters are given in Wm‘2. "Cloudiness" is the cloud 
amount given in percent from the sky camera observations, Low_neb, Mid_neb and 
High_neg are the cloudiness, in octas, obtained from nearby meteorological stations for 
three cloud layers as defined in the text. 
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the averaging. "Cloudiness" is the fractional 
coverage from the sky camera (in per cent), and 
the three last columns give, in octas, the mean 
cloudiness of the three layers defined above, as 
observed by the meteorological stations. 

• 

Unexpectedly, for a few days, the longwave or 
the shortwave cloud effect took negative values. 
The negative values for the longwave cloud 
effect are -1.9 W m‘2 and -0.4 W m’2 (twice), 
which is negligible and could be attributed to 
either the measurement or the model estimate. 
For two of these cases, cloudiness is very low, as 
expected since no cloud effect is computed. In 
the third case, however, cloudiness is larger 
than 50% and, therefore, we would have 


expected a larger longwave cloud effect. One 
notices the high value of the estimated "clear" 
flux found for that particular day (May 14). A 
careful analysis of the soundings showed that 
May 14 was exceptionally humid, which in turn 
explains the high value of the longwave 
irradiance, and the absence of a discernable 
longwave cloud effect. 

Seven days (out of 44) showed a negative 
shortwave cloud effect. The values decrease to 
-16Wm'2 and are only found for very low mean 
cloudiness (less than 5%). These erroneous 
value can be explained by model and 
measurement uncertainties when the actual 
shortwave irradiance is very close to the 




































































































































































































































































































































































































































theoretical "clear flux" (clear days). One of the 
major uncertainties for the modeled "clear" flux 
is the effect of aerosols. It is possible that the 
actual visibility was higher than that modeled 
for some of these days, leading to a higher 
"clear" flux. 

The largest shortwave cloud effect is found on 
June 2 , with an absolute value of 264 Wm-2 (less 
than 24% of the clear sky radiation reaches the 
ground). Instantaneous observations showed a 
very variable insolation with values constantly 
low and reaching extremely low levels. During 
a half hour period centered on local noon, for 
instance, only 1% of the clear sky irradiance was 
measured. A detailed analysis of other 
meteorological observations showed a very high 
cloudiness during most of the day and gusty 
winds frequently changing directions. This 
strongly supports the presence of Cb type clouds. 
Surprisingly, despite the obvious presence of 
very thick clouds, no rain was reported for that 
day over the site. It is possible, however, that 
very localized rainfall occurred, but did not 
reach the measurement station. 

The largest longwave cloud effect (94W m-2) is 
more than twice the value of the second largest 
(42 W m-2). n is found on October 10, an 
overcast day with a very low boundary layer 
temperature: An average Of 4°C of surface air 
temperature during daytime was measured, to 
be compared to about 15°C on average during 
the IFC. This low temperature leads to a record- 
low "clear" flux and explains the magnitude of 
the cloud effect. The shortwave cloud effect for 
this day is also one of the largest reported. Both 
effects partially compensated each other leading 
to a -27 W m"2 net cloud forcing. 

Net cloud radiative forcing ranges are from 
- 227 W m-2 to 37 W m"2, most days showing 
negative values. The few positive values are 
obtained for small mean cloud fractions. The 
highest negative values result mainly from the 
shortwave component. 

Correlations of radiative forcing with cloudiness 
values were done using the sky camera and the 
meteorological reports (Fig. 3a,b,c). Using the 
sky camera, we found a correlation of 0.62, 0.31 * 
and 0.49 with the shortwave, longwave, and net 
cloud radiative forcing, respectively. Using the 
observations from the meteorological stations, 
we found an even lower correlation. This 
shows, that cloud amount is not the only ' 
parameter driving the radiative forcing. • 
Another parameter governing shortwave cloud 
forcing is liquid water content, but this 
parameter was not measured. On the other 
hand, we know that the longwave flux should 
be more sensitive to the low cloudiness than to 
total cloudiness. We therefore correlated 








Fig. 3: Cloud forcing (in Wm’2) as a function of 
cloud amount (in %) as seen by the sky camera. 
The cloud forcing is given for the shortwave (Fig, 
3a), the longwave (fig. 3b) and the net (fig. 3c), 

longwave radiative forcing with low cloudiness; 
but, once again, the results were rather 
disappointing. 

Fig. 4 shows the average diurnal cycle of the 
shortwave, longwave, and net cloud radiative 
forcing as a function of day fraction. No values 
are given at night since the shortwave cloud 
forcing has no meaning. Sunrise and sunset 
occur approximately at day fractions 0.5 and 1.1, 
respectively. It is clear that the three quantities, 
which are not independent, show a diurnal 
cycle. The longwave radiative forcing has an 
amplitude of about 15 W m"^, with a minimum 
around local noon. The cycle of the shortwave 
cloud radiative forcing is even larger: it is 

minimum at sunset and sunrise but remains 
rather constant during a large part of the day. 
Note that, owing to our definition of the 
shortwave cloud radiative forcing, the 
shortwave cloud effect in Fig. 3 is not directly 
related to the actual amount of incoming 
radiation. As mentioned above, the accuracy of 
the shortwave radiative forcing is lower for low 
solar zenith angles than during the rest of the 


Daily cycle 



day, but we believe that the large cycle is not 
entirely due to uncertainties. The net cloud 
radiative forcing exhibits an even larger cycle, 
enlarged by the longwave component. 


5. SUMMARY AND CONCLUSIONS 

During FIFE, a large set of meteorological 
measurements was taken. "Clear" longwave 
and shortwave fluxes at the surface were 
computed using atmospheric soundings and 
radiative transfer models. Coincident radiative 
flux measurements were used to derive the 
surface cloud radiative forcing for both 
components. 

Both shortwave and longwave cloud forcings 
show a very large temporal variability due, for 
the most part, to the changing cloudiness, and to 
measurements and computational 
uncertainties. 

Daily values of the cloud radiative forcings were 
analyzed and compared to cloud parameters 
obtained from a sky camera situated on the FIFE 
site and from conventional meteorological 
observations in the immediate vicinity. The 
advantage of the sky camera was its immediate 
proximity to the surface radiative observations, 
but the stations near the FIFE area also gave 
information on cloud height, another 
important factor governing longwave cloud 
radiative forcing variability. 

No dominant correlation between either 
components of the cloud radiative forcing and 
cloudiness were found using the sky camera 
data. This is explained by the fact that forcings 
are a function, not only of cloudiness, but also of 
cloud liquid water content (mainly for the 
shortwave) and cloud base height (for the 
longwave). We expected that more results 
could be drawn from the conventional cloud 
observations which enabled distinction between 
low cloudiness and cirrus; once again, however, 
the correlation was very low. One reason for 
this is the poor spatial coherency of the 
cloudiness, especially for low clouds, as was 
shown by comparing two similar observations 
about 10 km apart. 

Surface cloud radiative forcing showed a rather 
large diurnal cycle with a minimum at local 
noon for the longwave component and two 
pronounced minimum at sunrise and sunset 
for the shortwave component. These cycles, if 
confirmed by other studies, should be taken into 
account when evaluating radiation budgets 
from space using h^liosynchronous satellites. 

One main deficiency of this study was the 
limited number of days during which 
observations were made. To provide 
convincing results on the mean radiative effect 
of clouds, at least one year of observations is 
needed. The measurements should include 
radiative fluxes, cloudiness, cloud base height, 
cloud optical thickness, as well as temperature 
and water profiles, and should be done within a 



/ very limited area. Such measurements are 
planned within the Atmospheric RadiatVm 
Measurement Program of the Department Of 
Energy (IX)E, 1990), which, therefore, should be 
useful to cloud radiative forcing studies. 
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